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Abstract 



Diffractive electroproduction of p and (f) mesons is measured at HERA with the HI 
detector in the elastic and proton dissociative channels. The data correspond to an integrated 
luminosity of 51 pb~^. About 10500 p and 2000 4> events are analysed in the kinematic 
range of squared photon virtuality 2.5 < < 60 GeV^, photon-proton centre of mass 
energy 35 < < 180 GeV and squared four-momentum transfer to the proton \t\ < 
3 GeV^. The total, longitudinal and transverse cross sections are measured as a function 
of Q^, W and The measurements show a transition to a dominantly "hard" behaviour, 
typical of high gluon densities and small qq dipoles, for larger than 10 to 20 GeV^. 
They support flavour independence of the diffractive exchange, expressed in terms of the 
scaling variable {Q'^ + My)/A, and proton vertex factorisation. The spin density matrix 
elements are measured as a function of kinematic variables. The ratio of the longitudinal 
to transverse cross sections, the ratio of the helicity amplitudes and their relative phases 
are extracted. Several of these measurements have not been performed before and bring 
new information on the dynamics of diffraction in a QCD framework. The measurements 
are discussed in the context of models using generalised paiton distributions or universal 
dipole cross sections. 
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1 Introduction 



Diffractive scattering is characterised, in high energy hadron interactions, by final states con- 
sisting of two systems well separated in rapidity, which carry the quantum numbers of the 
initial state hadrons. The process is related through unitarity to inelastic scattering and gov- 
erns the high energy behaviour of total cross sections. It is described in Regge theory f\\ by 
the exchange of the vacuum singularity, called the "pomeron", and may be interpreted as the 
differential absorption of the various virtual components of the interacting systems [|2l. It is a 
challenge for Quantum Chromodynamics (QCD) to explain diffraction in terms of quark and 
gluon interactions. 

Most diffractive phenomena - which include elastic scattering - are governed by large dis- 
tance, "soft" processes, which in general are not accessible to perturbative QCD (pQCD) cal- 
culations. However, for short distance processes, the presence of a "hard" scale offers the pos- 
sibility to use perturbative techniques to calculate diffractive amplitudes. Alternatively, at high 
energy the interaction properties of colour fields are invoked in models which characterize the 
incident particles as a superposition of colour dipoles with various size to calculate diffractive 
and total cross sections. 



An important testing ground for calculations in diffraction is provided by the study of ex- 
clusive vector meson (VM) production e + p — e + V + Y . This process is illustrated in 
Fig.[T] the intermediate photon of four- momentum q converts into a diffractively scattered VM 
(p, CO, (p, J /ip, ...) of mass My, while the incoming proton is scattered into a system Y of mass 
My, which can be a proton ("elastic" scattering) or a diffractively excited system ("proton dis- 
sociation"). In VM production, a hard scale can be provided by the photon virtuality Q, with 
Q'^ = —q^, the four-momentum transfer ^/\t\ from the proton, or by the quark mass (for heavy 
VM production). The reaction energy is defined by the photon-proton centre of mass energy 
W , with W"^ ~ Q"^ /x, where x is the Bjorken scaling variable. The high energy electron-proton 
collider HERA offers access to all these scales, over a wide range of values. 

The present publication is devoted to the study of the diffractive electroproduction of p and 
mesons, both for elastic and proton dissociative scattering. The data were taken at HERA 
with the HI detector in the period from 1996 to 2000. A common analysis of the four channels 




Figure 1: Diffractive vector meson electroproduction. 
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is performed. Measurements of the production cross sections and of the spin density matrix 
elements, which give access to the helicity amplitudes, are presented as a function of the kine- 
matic variables Q'^ (including the Q'^ dependence of the polarised cross sections), W , t and, for 
p mesons, the dipion mass. 

The measurement of kinematic dependences and the comparison between different VMs 
provide tests of a large spectrum of predictions. The data cover the interesting transition from 
the low domain, dominated by soft diffraction, to the higher Q"^ domain where hard diffrac- 
tion is expected to be dominant. This offers the opportunity to test the relevance of soft physics 
features present in the photon and VM wave function, and to study the development of features 
predicted by pQCD calculations. Quantitative tests of pQCD and colour dipole calculations are 
provided by the comparison with the data of various model predictions. Two important aspects 
of diffraction are tested: the flavour independence of the diffractive process and the factorisa- 
tion of the process into a hard scattering contribution at the photon vertex and soft diffractive 
scattering at the proton vertex ("Regge factorisation"). In addition, valuable information is pro- 
vided by precise measurements of empirical parameters, in particular the Q'^ and t dependences 
of the cross sections and the ratio of the proton dissociative to elastic cross sections, as well as 
the contributions of various backgrounds. 

The present studies confirm with increased precision previous HI measurements on p [|3]-[6l 
and electroproduction, mainly in the elastic channel but also in proton dissociation [[5][6||. 
The samples analysed here include data taken in 1996 and 1997, and the present results super- 
sede those presented in |l4l|5l|7]|. Thanks to the larger statistics, the scope of the investigation is 
significantly extended. 

This analysis complements other HI measurements of exclusive diffractive processes: pro- 
duction of real photons, in photoproduction (Q^ ~ 0) at large \t\ ^ and in electroproduction 
at small \t\ (deeply virtual Compton scattering - DVCS) [9], production of p mesons in pho- 
toproduction at low [10] and large \t\ IfTTTl . of J /il) mesons in photo- and electroproduction at 
low lfT2]| and large \t\ fO], of V^(2s) |[T4l and T IfTSl mesons in photoproduction. 

The ZEUS collaboration at HERA has performed measurements of DVCS fW\, p lfT7ti20ll. 
uj 11211121, [|20l|23l|2l, J/iJ [|20ll25ti27|| and T [28] production. Resuks at lower energy have 
been published, in particular for p electroproduction, by the DESY-Glasgow Il29ll , CHIO [|30ll , 
NMC [I3l]|, E665 (31 and HERMES (31 collaborations. The experimental and theoretical 
status of diffractive VM production before the high energy fixed target and HERA experiments 
is presented in detail in the review [|34il . 

The paper is organised as follows. The theoretical context and the models which will be 
compared to the data are presented in section [2l The HI detector and the event selection cri- 
teria are summarised in section [3l where the kinematic and angular variables are defined. The 
various signal samples are defined in sectionHl which also contains a detailed discussion of the 
backgrounds, a description of the Monte Carlo (MC) simulations used for the analyses and a 
discussion of the systematic errors affecting the measurements. In section [51 the measurements 
of the VM line shapes and of the elastic and proton dissociative cross sections are presented, 
and VM universality and proton vertex factorisation are discussed. Section [6] is devoted to the 
polarisation characteristics of the reactions and their kinematic dependence. A summary of the 
results and conclusions are given in section |71 
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2 Theoretical Context 



Since the first observation of high inclusive diffraction f3F] and of VM production at HERA, 
a large number of theoretical studies has been published on diffractive VM production (see 
e.g. [[36] - [63l ). Reviews of theoretical predictions confronted by the data have been published 
recently [l64l[65]l . 

2.1 Cross section calculations 

Calculations are performed following two main approaches, sketched in Fig. |2l The approach 
based on coUinear factorisation, illustrated in Fig. [2ta), describes VM production using the 
parton content of the proton, in the presence of a hard scale. The colour dipole picture of 
Fig. [2b) provides a complementary way to describe high energy scattering. 




Figure 2: Representative diagrams for diffractive VM electroproduction: a) the coUinear fac- 
torisation, GPD approach; b) the high energy, low x colour dipole approach. 



CoUinear factorisation In a pQCD framework, a coUinear factorisation theorem [|36ll has 
been proven for the production of longitudinally polarised VMs in the kinematic domain with 
W'^ > My, > t^QCD ^iid |t| < ^QCD^ for leading powers of Q and for all values of x. The 
longitudinal amplitude, sketched in Fig. [2a), is given by 

77 ^^^^(x; t) = y y dx' fi/p{x',x'-x,t;n)-Hij(Q'^x'/x,Q'^,z;iJ,)-'^J{z;n), (1) 

where fi/p{x', x' — x, t; /i) is the generalised parton distribution function (GPD) for parton i in 
the proton and /x is the factorisation and renormalisation scale, of the order of Q. The GPDs (see 
e.g. [[661 ), which are an extension of ordinary parton distribution functions (PDF), include cor- 
relations between partons with longitudinal momenta x and x' and transverse momenta t; they 
describe the off-diagonal kinematics {x' ^ x) implied by the different squared four-momenta 
of the incoming photon and outgoing VM. The Hi j matrix elements describe the hard scatter- 
ing from the parton i in the proton to the parton j in the meson with wave function /u), 
where z is the fraction of the photon longitudinal momentum carried by one of the quarks. The 
scale evolution is calculated, in the HERA kinematic domain, using the DGLAP equations, 
and higher order corrections have been calculated [[571 [581 [631 . CoUinear factorisation holds for 
heavy VMs [|36j[, and its validity is extended to transverse amplitudes at sufficiently high Q"^ 
(see e.g. [[36l[45l[48l[65l). 
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Dipole approach At high energy (small x) and small \t\, VM production is conveniently 
studied in the proton rest frame, for all values of Q^. It is described as three factorising con- 
tributions, characterised by widely different time scales [l67l[68l, as illustrated in Fig.|2tb): the 
fluctuation of the virtual photon into a qq colour dipole, with a coupling depending only on the 
quark charge, the dipole-proton scattering (either elastic or proton dissociative scattering), and 
the qq recombination into the final state VM. The amplitude is 

T^*P^^P{x; t) = dz j dh ^^{z, r) ■ ^''^"-^(a;, r; t) ■ ^^{z, r), (2) 

where r is the transverse distance between the quark and the antiquark, and '^'^{z, r) and 
'^^{z,v) are the photon and the VM wave functions, respectively. The diffractive dipole- 
proton cross section a'^^^^{x, r; t) is expected to be flavour independent and to depend only 
on the dipole transverse size (the impact parameter between the dipole and the proton is inte- 
grated over). Photons with large virtuality and fluctuations into heavy quarks are dominated by 
dipoles with small transverse size. In this case, the two quarks tend at large distance to screen 
each other's colour ("colour transparency" [69]), which explains the small cross section. In 
several models [|38ll43l , the convolution of the VM wave function with the dipole is expected to 
play a significant role in VM production, by selecting specific dipoles. It can be noted that the 
Generalised Vector Meson Dominance model B4[|50ll is related to the dipole approach. 

Dipole-proton scattering is modeled at lowest order (LO) in pQCD through the exchange of 
a gluon pair in a colour singlet state [ITOl . and in the leading logarithm approximation (LL 1/x) 
as the exchange of a BFKL-type gluon ladder. In a (z, k^) representation, where represents 
the quark (or antiquark) momentum component transverse to the photon direction (i.e. the 
Fourier transform of the dipole transverse size), fcj-unintegrated gluon distribution functions are 
used. The contributions of gluons with small kt are of a non-perturbative nature, whereas at 
large kt they can be obtained from the logarithmic derivative of the usual, integrated, gluon 
distribution, G{x,Q'^). In the LO and LL 1/x approximations both gluons emitted from the 
proton carry the same fraction x of the proton longitudinal momentum and the cross section 
is proportional to the square of the gluon density [4Tl,'42l. Calculations beyond the LL 1/x 
approximation take into account the difference between the longitudinal momentum fractions 
carried by the two gluons ("skewing" effects) [146116511711 . 

At low X, VM production can be calculated [[52l [561 [591 [621 , in the absence of a hard scale, 
using universal dipole-proton cross sections obtained from deep inelastic scattering (DIS) mea- 
surements [[72l . This approach automatically incorporates soft, non-perturbative contributions. 
Such models often involve parton saturation effects, expected from the recombination of high 
density gluons [73,74J as inferred from the observation of geometric scaling |[62l[74l . DGLAP 
evolution can also be included, for instance in the model [.59,1 . 

2.2 Kinematic dependences and crz/aT 

The photon-proton cross section can be decomposed into a longitudinal and transverse part, 
and ax, respectively. At LO and for t = 0, the dependences ctl oc 1/Q^ and ctt oc 1/Q^ are 
predicted [[42l . and the ratio R = aL/crx is predicted to be i? = Q'^ /My. Modifications to these 
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dependences are expected (see e.g. (431), due to the dependence of the gluon density, the 
quark transverse movement (Fermi motion) and quark virtuality [|49ll . and the dependence 
of the strong coupling constant a^. 

In the dipole approach, the square of the scale /i of the interaction is of the order of 

/i^ ~ z{\ - z)Q^ + kl + ~ z{l - 2)(Q2 + M^), (3) 

ruq being the current quark mass. It is related to the inverse of the relevant "scanning radius" [f38l 
|43l[65]| in the dipole-proton interaction. 

For longitudinally polarised photons or for heavy quark production, the qq wave function 
^'~^{z, r) is concentrated around 2 ~ 1 — ^ ~ 1/2. This suggests that a universal hard scale, 
fi, following from the transverse size of the dominant dipoles, can be of the order of jj? ~ 
{Q^ + My) / 4. For transverse photons fluctuating into light quarks, in contrast, the wave function 
is non-zero at the end-points z ~ or 1. These contributions correspond to small kt values of the 
quarks forming the dipole, and hence to a large transverse distance between them. The scale ji 
is therefore damped to smaller values than for longitudinal photons with the same virtuality, soft 
contributions may be significant and formal divergences appear in pQCD calculations for z — t- 
0, 1 1143 1144 II65II . For moderate values, the z distribution of light quark pairs from longitudinal 
photons can present a non-negligible smearing around the value z = 1/2, which results in 
a contamination of soft, "finite size" effects [641 • It is estimated that the fully perturbative 
QCD regime is reached for light VM production by longitudinal photons for Q"^ above 20 to 
30 GeV^ iHH- 

The W dependence of VM production is governed by the x^^ evolution of the gluon dis- 
tribution, with A increasing from ^ 0.16 for = 2 GeV^ to ^ 0.26 for = 20 GeV^, as 
measured in the total DIS cross sections at HERA ll75l . For heavy VMs and for longitudinally 
polarised light VMs at sufficiently high Q^, a strong ("hard") W dependence of the cross section 
is thus expected, fixed for all VMs by the scale (Q^ + My)/4. In contrast, the W dependence 
of the transverse cross section is expected to be milder than for longitudinal photons, since the 
A parameter is taken at a smaller value of the effective scale. This may result in a 14^ depen- 
dence of the cross section ratio R = ai/ar. In the framework of Regge theory, the existence 
of two pomerons [|76ll is postulated to describe both the soft and hard behaviours of the cross 
section [,5 3 J . 

At low |t|(|t|<0.5 — 0.6 Ge for elastic scattering) , the t dependence of VM production is 
well described by an exponentially falling distribution with slope b, da/dt oc I*' (predictions 
for the \t\ dependence are also given e.g. in II49[|54|]62II ). In an optical model approach, the slope 
b is given by the sum of the transverse sizes of the scattered system Y, of the qq dipole and of 
the exchanged system, with possibly in addition a VM form factor. Neglecting the latter, the 
t slopes for heavy VMs and for light VM production by longitudinally polarised photons are 
expected to take universal values, depending only on ((5^ + My)/4, whereas the production of 
light VMs by transverse photons, which is dominated by dipoles with larger transverse size, is 
expected to exhibit steeper t distributions II431I56II . This may result in a t dependence of (Tl/ctt. 

2.3 Helicity amplitudes 

The helicity amplitudes Tx^^x.,, where Ay and A^ are the VM and photon helicities, respec- 
tively, have been calculated in perturbative QCD for the electroproduction of light VMs with 
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\t\ <^ Il48ll49l[62l[65l . In this domain, the dominant amplitude is the s-channel helicity 
conserving (SCHC) Too amplitude, which describes the transition from a longitudinal photon 
to a longitudinal VM. Other amplitudes are damped by powers of Q. Those leading to the 
production of a transverse VM, of which the SCHC Tn amplitude is largest, contain an addi- 
tional factor oc l/Q. SCHC violation implies for single helicity flip amplitudes an additional 
factor oc ^/\t\/Q, to be squared for the double flip T_ii amplitude. This leads, in the kinematic 
range studied here, to the following hierarchy of amplitude intensities (assuming natural parity 
exchange): |Too| > |Tn| > |Toi| > |Tio|, |T_n|. 

2.4 Comparison of models with the data 

Predictions for VM production are available from a large number of models. Quantitative cal- 
culations generally imply the choice of PDF or GPD parameterisations or, in colour dipole 
models, of dipole-proton cross section parameterisations. Model calculations also generally 
imply the choice of VM wave function parameterisations, often taken as following a Gaussian 
shape, with several variants [|39ll40ll44l[56l[65ir72| . In view of the large number of models, no 
attempt is made in this paper to provide exhaustive comparisons to the data. Instead, a few mod- 
els and parameterisations, representative of recent approaches, are compared to various choices 
of observables. Examples of the uncertainties on the predictions, due to the choice of parton 
distribution functions and wave function parameterisations, are given for two of the models. 

• The GPD model of Goloskokov and KroU (GK [61]) provides predictions within the 
handbag factorisation scheme for the longitudinal and transverse amplitudes in the SCHC 
approximation. Soft physics is described by a GPD parameterisation of the proton struc- 
ture, constructed from standard PDFs with adequate skewing features and t dependences. 
The end-point singularities are removed with the aid of a specific model for the VM wave 
function. Error bands are provided with the model predictions. 

• The model of Martin, Ryskin and Teubner (MRT (45]) for p meson production is based 
on parton-hadron duality. Open qq production is calculated in an appropriate spin-angular 
state and in a specific invariant mass interval, which is then assumed to saturate p pro- 
duction, thus neglecting any VM wave function effects. The dependence of the gluon 
density, described by the anomalous dimension 7 with G{x,Q'^) oc {Q'^/Ql)"', is used 
to calculate the longitudinal and transverse cross sections. Skewing effects are parame- 
terised [46. 71,] without explicit use of GPDs. Predictions using two alternative PDFs are 
compared with the present data: CTEQ6.5M (721 and MRST-2004-NLO [|78l. 

• The model presented in the review of Ivanov, Nikolaev and Savin (INS f65\) is framed 
in the /c^-factorisation dipole approach. The helicity amplitudes are calculated pertur- 
batively and then extended into the soft region by constructing parameterisations of the 
off-forward unintegrated gluon density. The and W dependences of the cross sections 
and the full set of spin density matrix elements are predicted. Two wave function models, 
"compact" and "large", are used for p mesons, corresponding to two extreme cases for 
describing the p — )■ e+e" decay width. 
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• The /cj-factorisation calculations of Ivanov and Kirschner (IK f48l) provide predictions 
for the full set of helicity amplitudes, including helicity flip transitions. Similar to the 
MRT approach, the relevance of pQCD for transverse amplitude calculations is justified 
by the scale behaviour oc (Q^/Qo)^ of '^he gluon distribution, which avoids divergences 
for z 0, 1. 

• The dipole approach of Kowalski, Motyka and Watt (KMW f59\) uses an impact param- 
eter dependent description of the dipole cross section in the non-forward direction 11791 , 
within the saturation models of Golec-Biemat and Wiisthoff (GW [|73l ') and of lancu et al. 
(Colour Glass Condensate - CGC [74]). The and W dependences of the SCHC lon- 
gitudinal and transverse amplitudes are predicted using the DGLAP evolution equations 
for \t\ < 0.5 GeV^. 

• The dipole approach of Marquet, Peschanski and Soyez (MPS [|62l ) proposes an exten- 
sion of the saturation model [1741 . geometric scaling being extended to non-forward am- 
plitudes with a linear t dependence of the saturation scale. The exponential t dependence 
at the proton vertex is parameterised with a universal slope obtained from previous VM 
measurements. 



3 Experimental Conditions and Variable Definitions 

The diffractive production and decay of p and mesons is identified using the following reac- 
tions: 

e + p e + V + Y, 

p n+ + n- (i37^~100%), 

^ K+ + K- (i37^ = 49.2 ± 0.6%). (4) 

The events are selected by requiring the detection of the scattered electron and of a pair 
of oppositely charged particles, and by requiring the absence of additional activity in the de- 
tector, except in the region close to the outgoing proton beam, where proton dissociation can 
contribute. 

The kinematic domain of the measurements is: 

< 60GeV^ 

< 180 GeV, 

< 3GeV^ 

< 5 GeV. (5) 

The large values of W'^ compared to Q"^, My, My and \t\ ensure that the process is diffrac- 
tive, i.e. due to pomeron exchange. The variable xp = (Q^ -|- My + \t\)/{W'^ + Q'^ — My), 
which corresponds to the proton energy loss, is always smaller than 10~^. 



2.5 < Q2 

35 < 

\t\ 
My 
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3.1 Data sets 



The data studied here were taken with 27.5 GeV energy electrons or positrons colliding with 
820 or 920 GeV protons (in the rest of this paper the term "electron" is used generically to refer 
to both electrons and positrons). The data sets are summarised in Table [IJ where a/s is the 
ep centre of mass energy and the lepton beam type is specified. The integrated luminosity of 
51 pb^^ corresponds to running periods with all relevant parts of the detector fully operational. 
The periods with high prescaling of the triggers relevant for the present analyses are discarded. 
The published results with 1 < < 2.5 GeV^ are also presented in Table[I]("Hl SV"). 
They were obtained in 1995 in a special run of 125 nb~^, with the ep interaction point shifted 
by 70 cm in the outgoing p beam direction. This data set is not re-analysed in the present 
publication. 



Data taking 


lepton 


proton energy 




luminosity 


year 


beam 


(GeV) 


(GeV) 


(pb-i) 


1995 (SV) 


e+ 


820 


300 


0.125 


1996 


e+ 


820 


300 


4.0 


1997 


e+ 


820 


300 


9.8 


1999 




920 


320 


4.8 


1999 


e+ 


920 


320 


4.6 


2000 


e+ 


920 


320 


28.1 



Table 1: Characteristics of the data taken in 1995 with a shifted vertex (SV) and of the data sets 
used in the present paper (1996-2000). 



3.2 The HI detector and triggers 

A detailed description of the HI detector can be found in fSO]. Only the components essential to 
the present analysis are described here. The origin of the HI coordinate system is the nominal ep 
interaction point, with the positive z-axis (forward direction) along the direction of the proton 
beam. The polar angles 6 and the particle transverse momenta are defined with respect to this 
axis, and the pseudorapidity hrj = — logtan(6'/2). 

A system of two large coaxial cylindrical drift chambers (CJC) of 2 m length and 0.85 m 
external radius, with wires parallel to the beam direction, is located in a 1.16 T uniform mag- 
netic field. This provides a measurement of the transverse momentum of charged particles with 
resolution /\pt/pt — 0.006 pt © 0.015 {pt measured in GeV), for particles emitted from the 
nominal interaction point with polar angle 20 < 6* < 160°. Drift chambers with wires perpen- 
dicular to the beam direction, located inside the inner CJC and between the two CJC chambers, 
provide measurements of z coordinates. Track measurements are improved by the use of the 
central silicon tracker [[8T1l (from 1997 onward). The interaction vertex is reconstructed from 
the tracks. 

The liquid argon (LAr) calorimeter, located inside the magnet and surrounding the central 
tracker, covers the angular range 4 < 6* < 154°. The backward electromagnetic calorimeter 
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Spacal (153 < 6 < 177.5°) is used to identify scattered electrons. In front of the Spacal, the 
backward drift chamber (BDC) provides a precise electron direction measurement. 

The "forward detectors" are sensitive to energy flow close to the outgoing proton beam 
direction. They consist of the proton remnant tagger (PRT), a set of scintillators placed 24 m 
downstream of the interaction point and covering the angles 0.06 < 9 < 0.17°, and the forward 
muon detector (FMD), a system of drift chambers covering the angular region 3 < 9 < 17°. The 
PRT and the three layers of the FMD situated closer to the main calorimeter detect secondary 
particles produced in interactions with the beam collimators or the beam pipe walls of elastically 
scattered protons at large \t\ and of decay products of diffractively excited systems Y with 
My > 1.6 GeV. 

For the data collected in 1996 and 1997, events with > 2.5 GcV^ were selected by 
inclusive triggers requesting an electromagnetic energy deposit in the Spacal. For the years 
1999 and 2000, diffractive VM events with > 5 GcV^ were registered using several inclusive 
triggers; in addition, a special trigger was dedicated to elastic production with Q"^ > 2 GeV^. 

To reduce the data recording rate to an acceptable level, data selected by certain triggers 
have been dowscaled. In the following, the accepted events are weighted accordingly. 

3.3 Event selection 

For the present analyses, the scattered electron candidate is identified as an electromagnetic 
cluster with energy larger than 17 GeV reconstructed in the Spacal calorimeter. This energy 
threshold reduces to a negligible level the background of photoproduction events with a wrongly 
identified electron candidate in the Spacal. The electron direction is calculated from the position 
of the measured interaction vertex and from the BDC signals, when their transverse distance to 
the cluster barycentre is less than 3 cm; if no such BDC signal is registered, the cluster centre 
is used. 

The VM candidate selection requires the reconstruction in the central tracking detector of 
the trajectories of two, and only two, oppositely charged particles. They must originate from 
a common vertex lying within 30 cm in z of the nominal ep interaction point, and must have 
transverse momenta larger than 0.15 GeV and polar angles within the interval 20 < ^ < 160°. 
This ensures a difference in pseudorapidity of at least two units between the most forward track 
and the most forward cell of the LAr calorimeter. The VM momentum is calculated as the 
vector sum of the two charged particle momenta. 

The existence of a gap in rapidity between the VM and the forward system Y is further 
ensured by two veto conditions: that there is in the central tracker no additional track, except if 
it is associated to the electron candidate, and that there is in the LAr calorimeter no cluster with 
energy above noise level, E > 400 MeV, unless it is associated to the VM candidate. These 
requirements reduce to negligible level the contamination from non-diffractive DIS interactions, 
which are characterised by the absence of a significant gap in rapidity in the fragmentation 
process. They imply that the mass of the diffractively excited proton system is restricted to 
My ^ 5 GeV. They also contribute to the suppression of backgrounds due to the diffractive 
production of systems subsequently decaying into a pair of charged particles and additional 
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neutral particles. Energy deposits unrelated to the VM event and noise in the calorimeter are 
monitored from randomly triggered readouts of the detector. The energy threshold of 400 MeV 
leads to an average loss of 13t^l% °f diffractive VM events. 

A cut is applied to the difference between the sum of energies and the sum of longitudinal 
momenta of the scattered electron and VM candidate, E^E — p^) > 50 GeV. For events where 
all particles except the forward going system Y are detected, this quantity is close to twice the 
incident electron beam energy, 55 GeV. The cut reduces the QED radiation and background 
contributions in which additional particles remain undetected. 

3.4 Kinematic and angular variables 

To optimise measurements in the selected domain, the kinematic variables are reconstructed 
from the measured quantities following the algorithms detailed in [4J. In addition to the nominal 
beam energies, they make use of well measured quantities in the HI detector: the electron and 
VM directions and the VM momentum. 

The variable is reconstructed from the polar angles of the electron and of the VM ("dou- 
ble angle" method [|82l ). The modulus of the variable t is to very good precision equal to the 
square of the transverse momentum of the scattered system Y, which is calculated as the vector 
sum pt^miss = —{pty +Pt,e) of the transverse momenta of the VM candidate and of the scattered 
electronlll. The electron transverse momentum, pt^^, is determined using the electron energy ob- 
tained from the "double angle" method. The variable W is reconstructed from the VM energy 
and longitudinal momentum [i83i . The electron energy measured in the Spacal is used only for 
the calculation of the variable T.{E — p^). 

Three angles characterise VM electroproduction and two-body decay (Fig. [3]). In the helicity 
frame used for the present measurements, they are chosen as follows. The azimuthal angle 
is defined in the hadronic centre of mass system as the angle between the electron scattering 
plane and the VM production plane, which is formed by the directions of the virtual photon and 
the VM. The two other angles, which describe VM decay, are chosen in the VM rest frame as 
the polar angle 9 and the azimuthal angle ip of the positively charged decay particle, h^, the 
quantization axis being opposite to the direction of the outgoing system Y. 



4 Data Analysis 

This section first defines the analysis samples. The backgrounds are then discussed, the Monte 
Carlo simulations used to extract the signals are introduced, and the predictions are compared 
to the distributions of the hadronic invariant mass and of other observables. Finally, systematic 
uncertainties are discussed. 

'More precisely, the quantity \pt,miss P is a measure oft' = \t\ — \t\min, where \t\min is the minimum value of 
|t| kinematically required for the VM and the system Y to be produced on shell through longitudinal momentum 
transfer At HERA energies and for the relevant values of My and My, \t\min is negligibly small compared to \t\. 
In the following the notations |i| is used for t'. 
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hadronic centre of mass 




electron scattering plane production plane 



VM direction in 
the hadronic centre 
of mass system 




VM rest frame 

Figure 3: Definition of the angles characterising diffractive VM production and decay in the 
helicity system. 



4.1 Analysis samples 

Four event samples, which correspond approximately to the four processes studied in this paper, 
are selected following the conditions summarised in Tables [2]and[3l These conditions are chosen 
to minimize background contributions. 



Vector meson 


mass range 


p sample 
(j) sample 


0.6 < m^^ < 1.1 GeV 
1.00 < niKK < 1-04 GeV 



Table 2: Sample definition for the two VM selection. 



The VM identification reUes on the invariant mass of the two particles with trajectories 
reconstructed in the central tracker; no decay particle identification is performed. For the p 
sample, the mass m^^r calculated under the pion mass hypothesis is required to lie in the range 
0.6 < m^TT < 1.1 GeV. For the (p sample, the range 1.00 < niKK < 1-04 GeV is selected, the 
invariant mass itikk being calculated under the kaon hypothesis. 



Diffractive process 


forward detector selection 


t range 


notag sample 
tag sample 


no signal above noise 
signal detected above noise 




t 
t 


< 0.5 GeV^ 

< 3.0 GeV^ 



Table 3: Sample definition for the two diffractive processes. 
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The events in the p and samples are further classified in two categories, "notag" and "tag", 
according to the absence or the presence of activity above noise levels in the forward detectors, 
respectively. Elastic production is studied in the notag sample with |t| < 0.5 GeV^ whereas the 
tag sample with |t| < 3 GeV^ is used for proton dissociative studies. 



Year 


VM 


Q 


2 range (GeV^) 


W range (GeV) 


1995 - SV 




1.0 


< 




< 


2.5 


40 


< 


W 


< 140 


1996-1997 




2.5 


< 


Q' 


< 


4.9 


35 


< 


w 


< 100 






4.9 


< 




< 


9.8 


40 


< 


w 


< 120 






9.8 


< 




< 


15.5 


50 


< 


w 


< 140 






15.5 


< 




< 


27.3 


50 


< 


w 


< 150 






27.3 


< 




< 


60.0 


60 


< 


w 


< 150 


1999-2000 


(p notag 


2.5 


< 




< 


4.9 


35 


< 


w 


< 100 






4.9 


< 




< 


9.8 


40 


< 


w 


< 120 






9.8 


< 




< 


15.5 


50 


< 


w 


< 140 






15.5 


< 




< 


27.3 


50 


< 


w 


< 160 






27.3 


< 




< 


60.0 


60 


< 


w 


< 180 



Table 4: Kinematic range of the measurements. 



The kinematic domain of the measurements is summarised in Table IH It is determined by 
the detector geometry, the beam energies and the triggers, with the requirement of a reasonably 
uniform acceptance. The accepted range depends on the data taking period; for the notag (p 
sample in 1999-2000 it extends to smaller values than for the tag </> sample and for the p samples, 
due to the special elastic (p trigger. For W, the regions with good acceptance are determined by 
the track requirement; the accepted W values increase with Q'^ and with y/s. 

The acceptance increases with Q^, mostly because of the non-uniform geometric acceptance 
of the electron trigger for Q'^ < 20 GeV^. Monte Carlo studies show that the total acceptance 
mcreases from 15% (18%) for p (0) elastic production at = 2.5 GeV^ to about 50% at 

= 8 (6) GeV^ and to more than 60% for = 12 (10) GeV^ and that they are essentially 
independent of W in the measurement domain. 

The raw numbers of events selected in the four samples defined by Tables [2l|4] are given in 
Table [51 together with the numbers weighted to account for the downscaling applied to certain 
triggers. 



Numbers of events 


p sample 


sample 


raw 


weighted 


raw 


weighted 


notag sample 
tag sample 


7793 
2760 


11775 
3824 


1574 
416 


1976 
495 



Table 5: Events in the different data samples: raw numbers and numbers weighted to account 
for the downscaling applied to certain triggers. 
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p notag 

\t\ < 0.5 GeV^ 


ptag 

\t\ < 3 GeV^ 


(j) notag 
\t\ < 0.5 GeV^ 


0tag 

\t\ < 3 GeV^ 


p. diss, events 
el. events 

0^ Stt 

p' 


10.7 ± 0.3% 

0.3 ± 0.1% 
0.6 ± 0.1% 
4.0 ± 0.2% 


13.1 ± 0.5% 

0.4 ± 0.1% 
0.7 ± 0.1% 
7.7 ± 0.4% 


9.7 ± 0.7% 

6.3 ± 0.5% 

1.7 ± 0.3% 
3.6 ± 0.4% 


11.8 ± 1.5% 

4.7 ± 0.9% 

2.8 ± 0.7% 
9.2 ± 1.3% 



Table 6: Background contributions to the four data samples defined in Tables [21111 The quoted 
errors are the statistical errors from the MC samples. 



4.2 Backgrounds 

Several background processes, which affect differently the four data samples and depend on 
the kinematic domain, are discussed in this section. Their contributions are summarised in 
Table[6l The non-resonant tttt contribution to the p signal, which contributes essentially through 
interference, is discussed separately in section U.2.1[ The e+e~ and backgrounds were 

found, using the GRAPE simulation [f85ll . to be completely negligible. 

4.2.1 Cross-contaminations between the elastic and proton dissociative processes 

The notag and tag samples correspond roughly to the elastic and proton dissociative processes, 
respectively. However, cross-contaminations occur, due to the limited acceptance and efficiency 
of the forward detectors and to the presence of noise. The response of these detectors is modeled 
using independent measurements, by comparing signals in the various PRT and FMD planes. 

The cross-contaminations are determined for each VM species without a priori assumptions 
on the relative production rates of elastic and inelastic events. In a first step, the contaminations 
are calculated from the numbers of tag and notag events and from the probabilities for elastic 
and proton dissociative events to deposit a signal in the forward detectors as obtained from the 
MC simulations. The crossed backgrounds are then determined in an iterative procedure from 
the simulations, after final tunings to the data. 

Proton dissociative backgrounds in the notag samples Proton dissociative events produce 
a background to the elastic signals in the notag samples when the mass of the excited baryonic 
system is too low to give a signal in the forward detectors (My ^1.6 GeV) or because of 
inefficiencies of these detectors. The background fraction increases strongly with because 
the proton dissociative cross sections have a shallower \t\ distribution than the elastic cross 
sections. In the notag samples with \t\ < 0.5 GeV^, the proton dissociative background amounts 
to 10.5%. 
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Elastic backgrounds in the tag samples Conversely, elastic background in the proton disso- 
ciative samples of tag events is due to unrelated signal or noise in the forward detectors. For 
|t| < 3 GeV^, it amounts to 12.5%, with larger contributions for small |t| values where the 
elastic to proton dissociative cross section ratio is larger. In addition, when \t\ is large enough 
for the scattered proton to hit the beam pipe walls or adjacent material (\t\^ 0.75 GeV^), elastic 
events may give signal in the forward detectors. 

4.2.2 Cross-contaminations between p and <p samples 

For p production, the contribution from the — i- K^K^ channel is removed by the require- 
ment m^^ > 0.6 GeV, which also suppresses the contribution of the — i- K^Kl channel 
(BTZ = 34%) with the Kg meson decaying into a pion pair close to the emission vertex and the 
Kl being undetected in the calorimeter. 

The largest background in the selected (p samples is due to the low mass tail of tt+tt^ pair 
production extending under the (p peak. It amounts to 6% and depends on Q^. The shape of the 
TT+TT^ distribution corresponding to small values of rnxK is discussed in section |431 

4.2.3 — )■ 3 TT and oj backgrounds 

A small contamination in the p samples is due to the channel — )■ tt+tt^tt^ (BTZ = 15%) 
when each photon from the vr*^ decay remains undetected because it is emitted outside the LAr 
calorimeter acceptance, because the energy deposit in the LAr calorimeter does not pass the 
400 MeV threshold, or because it is associated with one of the charged pions. This background 
contributes to the m^jr distribution mostly below the selected mass range; it amounts to 0.3% of 
the selected p notag sample with \t\ < 0.5 GeV^ and 0.4% of the tag sample with |t| < 3 GeV^. 
The background rate increases with \t\ because the non-detection of the vr" decay photons leads 
in general to an overestimate of the pt imbalance of the event, pt^miss, which mimics a large 
value. The — 3 vr contribution below the — )■ KK signal is negligible. 

Similarly, the diffractive production of u mesons decaying in the mode cu — )■ 7r"'"7r~7r° 
(BTZ = 89%) gives background contributions to the p and samples when the vr*^ decay pho- 
tons escape detection. In addition, the tu — )■ tt+tt' (BTZ = 1.7%) channel gives an irreducible 
background to the p signal. The background due to tu production contributes 0.6% to the elastic 
and 0.7% to the proton dissociative p samples, and 1.7 and 2.8% for the samples, respectively. 
The non-detection of photons leads to large reconstructed |t| values for these contributions. 
Note that for the cross sections quoted below, as for results in previous HERA papers, the cu — p 
interference is neglected: its contribution is small and cancels when integrated over the mass 
range. 

4.2.4 p' background 

The largest background to the p signal and the second largest background to the signal is due to 
diffractive p' production^. The p' mesons decay mostly into a p meson and a pion pair, leading to 

^The detailed mass structure f84l of the states described in the past as the /9'(1600) meson is not relevant for 
the present study. The name p' is used for all VM states with mass in the range 1.3 — 1.7 GeV. 
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final states with four charged pions (p' — i- p°7r+7r^) or with two charged and two neutral pions 
{p' — )■ p='=7r^7r°, — !• Ti^u:^). The 7r+7r^7r°7r° events can mimic large \t\ p or production 
when the photons from the 7r° decays escape detection, which induces a pt imbalance in the 
event and a distortion of the t distribution, similarly to the — )■ 3 tt and — )■ 3 tt backgrounds. 
At high \t\, this background affects mostly the notag samples. It is indeed distributed between 
the notag and tag samples following the elastic to proton dissociative production cross section 
ratio, whereas genuine high \ t\ p and </> mesons are essentially produced with proton dissociation 
and thus contribute mainly to the tag samples. 

No cross section measurement of diffractive p' production has been published in the relevant 
range. The p' contribution to the p signal is thus determined from the data themselves, using 
a method presented in the HI analysis of high \t\ p electroproduction [5J. The distribution 
of the variable which is the cosine of the angle between the transverse components of the 
p candidate momentum, pt^p, and of the event missing momentum, pt^miss, is sensitive to the 
relative amounts of p signal and p' background. The p' contribution gives a peak at ( = +1 and 
a negligible contribution at ^ = —1, since the p and the missing 7r°'s are all emitted roughly in 
the direction of the p'. In contrast, the p signal gives peaks at ^ = +1 and ^ = — 1. However, 
for genuine p production, is also correlated to the angle </> between the p production plane 
and the electron scattering plane, which is distributed according to the a priori unknown value 
of the combinations of spin density matrix elements Tqq + 2rli and Tqq + 2r{^ (Eq. (|4T1) of the 
Appendix). 

An iterative procedure is used to determine simultaneously the amounts of p' background in 
the notag and tag samples, the matrix element combinations Tqq + 2rf^ and Tqq + 2r{^ (assumed 
to be identical for elastic and proton dissociative scattering), and the \t\ distributions of p elastic 
and proton dissociative production. It is found to converge after a few steps. The results are 
also used to calculate the p' background to the signal. 

The p' background is estimated to contribute 4% to the notag samples with \t\ < 0.5 GeV^, 
and 8% to the tag samples with |t| < 3 GeV^. 

4.3 Monte Carlo simulations 

Monte Carlo simulations based on the DIFFVM program are used to describe p, to, and p' VM 
production and decay, detector response (acceptances, efficiencies and variable reconstruction) 
and radiative effects. 

The DIFFVM program lf86l is based on Regge theory and Vector Meson Dominance [|87ll . 
The My diffractive mass distribution for proton dissociative events contains an explicit sim- 
ulation of baryonic resonance production for My < 1.9 GeV and a dependence da/dMy oc 
1 /My^^ for larger masses If88l , with quark and diquark fragmentation simulated using the JET- 
SET programme (E9]. 

The p and MC samples are reweighted according to the measurements of the Q'^, W and 
\t\ differential cross sections and of the angular VM production and decay distributions: the 
angle 9 is distributed according to the measurements of the Tqq matrix element (Eq. (|39l)). the 
angle to those of the + 2rl^ and Tqq + 2rl^ combinations (Eq. (|4TI)). and the angle cp to 
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those of the cos S parameter, which in the SCHC approximation fixes the ip = (p — ip distribution 



For the u and p' backgrounds, the cross section dependences on the kinematic variables 
Q^, W and \t\ are taken to be the same as for p mesons at the same (Q'^ + My)/A value. 
For the two-body tu decay, the angular distributions are taken as for p mesons. For three- 
body tu and decays, the angular distributions are chosen to follow and cos^^ distribu- 
tions described by the same values of the matrix elements as for two-body decays. For p' 
decay^ the parameters Mi(OO) and Mi(lO) describe the angular distributions [|90l . The values 
|Mi(00)p = 0.5, |Mi(10)p = 0.5 are chosen for the present simulations. 

The ratio of proton dissociative to elastic cross sections is taken from the present p analysis 
and assumed to be the same for all VMs. All kinematic and angular distributions are taken to be 
identical for elastic and proton dissociative scattering, as supported by the present data, except 
for the \t\ dependence of the cross sections. 

The to p cross section ratio is set to that measured in this analysis. The w to p ratio is 
taken from ZEUS measurements ll2Tll22l . For p' production, a p' to p ratio of 1.12 is usecfl, as 
a result of the procedure described in the previous section. 

For p, and cu mesons, the particle mass, width and decay branching ratios are taken from 
the PDG compilation fi84J . The mass and width of the p' resonance are taken as 1450 MeV 
and 300 MeV, respectively. For p and meson decays into two pseudoscalar mesons, the mass 
distributions are described by a relativistic Breit-Wigner function BW{m) with momentum 
dependent width, as described in section |5l2| In addition, the p mass shape is skewed according 
to the parameterisation of Ross and Stodolsky [|91| , 



with the Q"^ dependent value of n measured in this analysis. 

The TTTT background in the mass region is taken from the skewed Breit-Wigner distribution 
for p mesons, modified for m^^ < 0.6 GeV according to the empirical form 



with masses expressed in GeV and the parameter k being taken to be 1.5. This parameterisation 
describes the low mass 171,,^ distribution well, as shown in Figs. HI [5] and |9l where the cut 
rriKK > 1-04 GeV is applied to suppress genuine production. 

^In the dominant p' pmr decay mode, the two pions do not form a p resonance and can be assumed 
to be in a spin state. The angular decay distribution thus includes the two possible polarisation states of the p 
meson, with the squared amplitude \Mi (00) p (| Afi (10) p) corresponding to the probability that it is longitudinally 
(transversely) polarised, giving in the SCHC approximation, with the notations of the Appendix: 

W{9,i;) = ^ — L_ {|Mi(00)|2 [^shi^9 + eRcos^0-§sin29cos^cosS + ^shi^9cos2'iJj] 

+ |Afi(10)p [ ^{l + cos^9 + eRsm^9+§sm29cosijcosS- § sin^ 6* cos 2V^ ]}, where A' = y^2eR{l + e). 

"^This number does not constitute a p' cross section measurement, but it is used as an empirical parameterisation 
for describing the p' background contribution under the p peak, for the p' mass and width chosen in the simulation; 
as a consequence, vaiying the latter values has negligible influence on the background subtraction. 



(Eq. (HT])). 




(6) 




(7) 
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Radiative effects are calculated using the HERACLES program [|92l . Corrections for these 
effects in the selected kinematic range with T.(E — p^) > 50 GeV are of the order of 1%. 

All generated events are processed through the full GEANT ||93l based simulation of the 
HI apparatus and are reconstructed using the same program chain as for the data. Of particular 
relevance to the present analysis is the description of the forward detector response; the activity 
in these detectors, not related to VM production, is obtained from data taken independently of 
physics triggers, and is superimposed on generated events in the MC simulations. 

4.4 Mass distributions 

The mTTTT and uikk mass distributions are shown in Figs. |4]to[6l separately for the notag and 
tag samples. The results of the Monte Carlo simulations, comprising signal and backgrounds, 
are also shown. They are reweighted and normalised to the data as described in the previous 
section. 

The mass spectra are presented from threshold to masses well above the actual measurement 
ranges defined in Table[2l The 171^,^^ spectra in Figs. |4] and [5] are presented in four bins in \t\, with 
the cut itlkk > 1-04 GeV. The itlkk spectra in Fig. [6] exhibit the reflection of p production and 
of backgrounds. 

The m^TT mass distributions are well described from the threshold at 2m7r up to 1.5 GeV. 
The backgrounds are small in the mass ranges selected for the physics analyses, shown as the 
shaded regions in the figures, but their contributions can be distinctly identified outside these 
domains. In the m^^ distributions, they are particularly visible at low mass and, as expected, 
they contribute mostly at large \t\, especially in the notag sample with \t\ > 0.5 GeV^ of Fig.|4l 
A decrease of the background with increasing Q"^ for the same ranges in \t\ is also observed 
(not shown here), which is explained by the larger transverse momentum of the virtual photon, 
resulting in larger pt values of the decay photons which thus pass the detection threshold and 
lead to the rejection of the events. 

In view of the small p' background in the final selected samples, an analysis of only the mass 
spectrum, performed in the restricted mass range 0.6 < mT,T, < 1.1 GeV, is not sufficient to 
constrain the p' contribution. Controlling this background is crucial for the measurements of the 
\t\ slope and of the Tqq matrix element. In the present analysis, the amount of p' background is 
obtained from the distribution of the variable C (defined in section |4! 2. 41 ). The value determined 
within the mass range 0.6 < m^rTr < 11 GeV also gives a good description of the mass range 
2m J, < m^TTr < 0.6 GeV, below the actual measurement. This demonstrates the reliability of 
the background estimate. 

The rnxK mass distribution shown in Fig.[6]is also very well described. The nir background 
under the (f) peak, which contains a p' contribution obtained from the p analysis, is smal|j. 

''For the notag sample with |<| < 0.5 GeV^, the background under the (f> peak amounts to 20.5% for ~ 

2.5 GeV^ (10% from the tt+tt^ low mass tail, 3% from uj and 7.5% from p' production), and to 5.5% for — 
13 GeV^ (2.5%, 0.5% and 2.5%, respectively). An empirical description of the background by ZEUS, using a 
simple power law shape, is in agreement with these detailed findings: it amounts to 18% for — 2.5 GeV^ and 
5%forQ2 ^ i3GeV^ lEl. 
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Figure 4: Distributions of the invariant mass m^^r (with the cut ttikk > 1-04 GeV to reject the 
— )■ KK signal) in four domains in \t\, for the notag sample. The dashed histograms show 
the MC predictions for the p' background, the dotted histograms the sum of the p', u and (p 
backgrounds, and the full histograms the p signal (including interference with irn non-resonant 
production) and the sum of all backgrounds. The mass and \t\ domain where the cross section 
measurements are performed is shaded. 
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Figure 5: Same as in Fig. |4l for the tag sample. 
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Figure 6: Distributions of the invariant mass rriKK' (upper plots) in the (j) mass region, for the 
notag and tag samples separately; (lower plots) over an extended mass range, showing the 
signal and the reflection of p production and the backgrounds. The dashed histograms show 
the sum of the p', cu and — > 3 vr backgrounds, the dotted histograms show in addition the p 
and non-resonant tttt backgrounds, and the full histograms the — )■ KK signal and the sum 
of all backgrounds. In (a) and (b), the mass domain where the cross section measurements are 
performed is shaded. 
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4.5 Kinematic and angular distributions 

Figures |7] and [8] present several kinematic and angular variable distributions for the samples 
selected as defined in Tables [IB They demonstrate that the simulations, taking into account 
the detector acceptance and response and the background contributions, correctly describe the 
data. 

Figure |7] shows kinematic variable distributions of the p and (j) notag samples. The structure 
observed in the electron polar angle distribution (a) results from the different kinematic range 
selections for the different years. The dip in the distribution (b) of the laboratory azimuthal 
angle 0e of the electron is due to an asymmetric electron trigger acceptance. The pt distributions 
of the decay mesons (e), (i) reflect the VM mass and the decay angular distributions. The good 
description of the difference between the azimuthal angles of the decay kaons in the sample (f) 
indicates that the reconstruction of pairs of tracks with small differences in azimuthal angles is 
under control. A description of similar quality is obtained for the tag samples. 

Figures [8];a)-(c) present distributions related to the spin density matrix elements. The C, 
distributions (d)-(e) are sensitive to the values of the matrix element combinations Tqq + 2r\^ 
and Tqq + 2r\^ and to the amount of p' background especially at high \t\ d& discussed in sec- 
tion |4]2]11 The \t\ distributions (f)-(i) are sensitive to the amount of diffractive backgrounds 
(proton dissociation for the notag sample, elastic scattering for the tag sample) and to the values 
assumed for the exponential t slopes. 

4.6 Systematic errors 

Uncertainties on the detector response and background contributions are listed in Table |71 They 
are estimated by variations in the MC simulations within the indicated limits, which are in most 
cases determined from the data. Global normalisation errors are given separately. 

The error on the electron polar angle Of., which affects the measurements and the accep- 
tance calculations, is due to the uncertainty on the absolute positioning of the BDC with respect 
to the CJC chambers, the uncertainty on the electron beam direction in the interaction region 
and the error on the z position of the interaction vertex. 

The uncertainty on the energy scale of the Spacal calorimeter affects the cross section mea- 
surements through the electron energy threshold of 17 GeV and the — Pz) cut. 

The uncertainty on losses due to the rejection of events affected by noise in the LAr calorime- 
ter or containing energy deposits unrelated to the diffractive event is estimated by varying the 
energy threshold, both in the data and in the simulation (where data taken from random triggers 
are directly superimposed to the simulated events). 

The uncertainties on the simulated cross section dependences on Q^, W and \t\ affect the 
bin-to-bin migrations and the extrapolations from the average value of the kinematic variables 
in a bin to the position where they are presented ("bin centre corrections"). 

An absolute error of ±0.10 is used for the ratio of the proton dissociative (with My < 
5 GeV) to elastic cross sections, which corresponds to about 20% relative error. It is estimated 
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Figure 7: Distributions of the polar angle 6e (a) and azimuthal angle (b) of the scattered elec- 
tron, of the Q"^ (c) and W (d) variables, and of the transverse momenta of the decay mesons (e), 
for the p notag sample; distributions of the difference between the azimuthal angles cj) of the 
decay kaons (f) and, in (g)-(i), of the same observables as in (c)-(e), for the (p notag sample. In 
panels (a)-(e), the dashed histograms present the MC predictions for the distributions of the p' 
background, the dotted histograms in addition for the to and backgrounds, and the full his- 
tograms for the p signal and the sum of all backgrounds; in panels (f)-(i), the dashed histograms 
describe the p' and to backgrounds, the dotted histograms in addition the nn background, and 
the full histograms the signal and the sum of all backgrounds. 
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Figure 8: Distributions of the VM production and decay angles (j) (a), cos 6 (b) and = (f)—ip (c) 
for the p notag sample; of the C variable for the p notag (d) and tag (e) samples; of the \t\ 
variable for the p (f) and </> (g) notag samples and for the tag samples (h)-(i). In panels (a)-(e), 
the dashed histograms present the MC predictions for the distributions of the p' background, 
the dotted histograms show in addition the to and backgrounds, and the full histograms the p 
signal and the sum of all backgrounds; in panels (f)-(i), the dotted histograms show the sum of 
the various VM backgrounds (p', cu, or p + vrvr), the dash-dotted histograms show in addition 
the diffractive background (proton dissociation in panels (f)-(g) and elastic production in panels 
(h)-(i)), and the full histograms the signal and the sum of all backgrounds. 
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Uncertainty source 


functional dependence 


VM 


variation 


Detector effects 


electron polar angle 6e 




all VM 


±1 mrad 


Spacal energy scale 




all VM 


±1% 


noise threshold in LAr 




all VM 


±100 MeV 


Cross section dependences 


da/dQ2 


(g2 + M2)-" 


all VM 


n±0.15 


da/dW 




all VM 


6 ± 25% 






P 

0, p', U 


el. : 6 ± 0.5 GeV"^ 
p. diss. : 6 ± 0.3 GeV"^ 

el. : 6± 1.0 GeV-2 
p. diss. : 6±0.7GeV"^ 


Backgrounds 


proton dissoc. / elastic 




all VM 


±0.10 ±20%) 


p shape skewing 




P 


n±0.15 


VM cross sections 




uj/p 


±0.02 (^ ±20%) 


p' decay 


Ml (00) and Ml (10) 


Hp 
p'/p 
P' 


±0.03 (^ ±15%) 
±0.40 (^ ±35%) 
see text 


p and (j) angular decay distributions 


'00 


fm 

f{\t\) 


p, 4> 

P, 


±15% 
±30% 


cos 5 







±0.05 


Global normalisation 


luminosity 




all VM 


±1.5% 


trigger efficiency 




all VM 


±1.0% 


track rec. eff. (per track) 




all VM 


±2% 


width of rel. B.-W. 


see text 


P 


±2% 


(j)^ KK Bn 


see [.84)1 





±1.2% 


TTTT under peak (k param.) 







±100% 


da/dM^ 


1/M|-" 


all p. diss. 


n±0.15 



Table 7: Variations in MC simulations for the estimation of systematic uncertainties. Numbers 
between parentheses indicate the relative variations. 
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by varying by ±0.15 the parameter n in the simulated dissociative mass distribution da/dMy oc 
1 /My', by varying the slope parameters of the exponential \t\ distributions of elastic and proton 
dissociative events within the experimental limits, and by calculating the cross section ratio 
using only the PRT or only the FMD. The latter covers uncertainties in the inefficiencies of 
these detectors. 

For the p cross section measurements, the error due to the extraction of the non-resonant tttt 
background is estimated through the variation of the dependent skewing parameter n of the 
Ross-Stodolsky parameterisation of Eq. 

The errors on the various cross section ratios are taken from the present analysis for the (p to 
p and p' to p ratios, and from the ZEUS measurements of the cu/p ratio [|22ll . 

The errors due to the uncertainty on the p' decay angular distribution are estimated by con- 
sidering the two extreme cases |Mi(00)p = 1, |Mi(10)p = and |Mi(00)p = 0, |Mi(10)p = 
1 of the pair of variables defined in [|90ll . 

The uncertainty on the angular distributions are described by varying the values of the matrix 
element Tqq (for the angle 9), of the combinations Tqq + 2r^^ and Tqq + 2rl^ (for the angle 0) 
and of the cos 5 parameter (for the angle ip = (p — (f). 

The uncertainty on the choice of the momentum dependent width of p mesons results in 
normalisation uncertainties of 2% (see section [5^ . 

For production, the uncertainty on the tttt background under the signal is estimated by 
varying the parameter k globally from to 3 (Eq. ^ in section |431 ), leading to a normalisation 
error of ±3% on the cross section measurements. 

For the proton dissociative cross sections, the error on the correction for the smearing 
through the experimental cut My < 5 GeV is estimated by varying the parameter n of the 
My distribution (da/dMy oc 1/My", with n ± 0.15), which leads to an additional normalisa- 
tion error of ±2.4% on the proton dissociative cross section measurement. 

The uncertainties on the luminosity measurement, on the triggers and on the track recon- 
struction efficiency are assumed to affect globally the normalisation only. 

Systematic errors due to limited MC statistics are negligible compared to the statistical 
precision of the measurements (the generated samples correspond to at least ten times the data 
integrated luminosity). 

All systematic errors on the measurements presented in the rest of this paper are calculated 
from separate quadratic sums of positive and negative effects of the variations listed in Table |7J 
In all figures, measurements are shown with statistical errors (inner error bars) and statistical and 
systematic errors added in quadrature (full error bars). In tables, the errors are given separately: 
first the statistical, second the systematic errors. Overall normalisation errors are not included 
in the error bars but are quoted in the relevant captions. 
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5 Cross Section Results 



In this section, measurements of the p and line shapes are presented first. The elastic and pro- 
ton dissociative cross sections are then measured as a function of (total and polarised cross 
sections), W and t (total cross sections); results for different VMs are compared. Finally, elastic 
and proton dissociative scatterings are compared, including tests of proton vertex factorisation. 
Model predictions are compared to the data. 



5.1 Measurement of cross sections 

The cross sections for p and cp production presented in this paper are extracted from the num- 
bers of events in the mass ranges 0.6 < 777,,^^ < 1.1 GeV and 1.00 < rnxK < 1-04 GeV, 
respectively. They are corrected for all backgrounds, including for p mesons the non-resonant 
dipion diffractive production (see section [5.2. II) . They include all corrections for detector ac- 
ceptance and response. When quoted at a fixed value of a kinematic variable, the cross sections 
are evolved from the average value in the bin using dependences measured in this analysis. 

The cross sections are quoted for the full resonance mass range from the two particle thresh- 
old up to the nominal mass plus five times the resonance width: 

2m^< m^^ <mp + hVp ~ 1501 MeV, 
2mK< rriKK <m^ + 5r^ ~ 1041 MeV. (8) 

For mesons, the cross sections take into account the branching ratio to the K^K^ channel. 

Elastic and proton dissociative cross sections are given at the Born level (i.e. they are cor- 
rected for QED radiation effects) in terms of 7*p cross sections (except for the mass shapes, 
which are given in terms of ep cross sections). The 7*p cross sections are extracted from the 
ep cross sections in the Weizs acker- Williams equivalent photon approximation [|94| using the 
definition 

1 dV(e + p-^e + l^ + y) 
f d^dQ^ 
where the flux F of virtual photons [|95l and the inelasticity y are given by 

^ ^ 1 -1/ + 1/V2 ^ p^ ^^^^ 

vr yQ"^ ' p-k' 

being the fine structure constant and p and k the four-momenta of the incident proton and 
electron, respectively. 



a(7* +p^V + Y) = t,. ^ (9) 



5.2 Vector meson line shapes 

The distribution of the invariant mass m of the VM decay particles is analysed assuming the 
relativistic Breit-Wigner distribution BW{m) with momentum dependent width T(m) [i96il : 

, , m My Tim) 

F m = Tv [ — ] , (12) 
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where My and Ty are the nominal VM resonance mass and width, q* is the momentum of the 
decay particles in the rest frame of the pair with mass m, and Qq is the value taken by q* when 
m = My. 

For p mesons, the mass extrapolation from the measurement domain 0.6 < m^^ < 1.1 GeV 
to the full range given by Eq. ([8]), including the correction for skewing effects, implies a correc- 
tion factor of 1.15 with a systematic error of 2% due to the theoretical uncertainty on the choice 
of the momentum dependent width [|96|. For cp production, a very small extrapolation outside 
the measurement domain is required, with negligible related error. 



5.2.1 p mesons 



Distributions of the m^T, mass in the range 2 m^, < m^^, < 1.5 GeV, with the cut rriKK > 
1.04 GeV applied to suppress the signal at low mass, are shown in Fig. |9] for elastically 
produced events in four ranges in Q'^, after subtraction of the proton dissociative, </>, cu and p' 
backgrounds and corrections for detector and QED radiation effects. The mass resolution in the 
p mass range, determined with the MC simulation, is about 10 MeV. 



Skewing The mass distributions are skewed towards small masses, especially at low Q^. Ac- 
cording to Soding's analysis [97J, this is due to the interference of the p meson with back- 
ground from p-wave Drell-type non-resonant vrvr pair production, with positive interference for 
m^TT < nip and negative interference for m^j, > nip. 

Following one of the forms of skewing proposed in IfTTll . the p mass shape is described as 



oc 



m^p - ml^ + i mpT^rriT, 



2 



(13) 



where resonant and non-resonant vrvr production are supposed to be in phase. The interference is 
proportional to //, which is taken to be independent of the m^jr rnass; the very small purely non- 
resonant contribution is given by /|/4. Figure [9] shows that the p mass shape is well described 
by Eqs. (|1 1H13I) over the full range 2m^ < m-,,-,, < 1.5 GeV, with the skewing parameters 
fitted in the range 0.6 < ttIt^^ < 1.1 GeV. No indication is found for significant additional 
backgrounds, also outside the mass domain used for the measurements. The p skewing effect is 
also often conveniently parameterised in the form proposed by Ross and Stodolsky [|9T1l , given 
by Eq. ©. 

For a fit over the whole range with the parameterisation of Ross and Stodolsky, the values 
of the resonance mass and width are 769 ±4 (stat.) MeV and 162 ±8 (stat.) MeV, respectively. 
The Soding parameterisation gives similar values, with larger errors. This is in agreement with 



the world average values as obtained in photoproduction II84I : rup = 768.5 ±1.1 MeV and 
Tp = 150.7 ± 2.9 MeV. Within errors, no difference is observed between the elastic and proton 
dissociative samples. 



30 




0.25 0.5 0.75 1 1.25 1.5 0.25 0.5 0.75 1 1.25 1.5 



m,, [GeV] m,, [GeV] 




0.25 0.5 0.75 1 1.25 1.5 0.25 0.5 0.75 1 1.25 1.5 



m,, [GeV] m,, [GeV] 

Figure 9: Distributions of the m^jr rnass for elastic p production with \t\ < 0.5 GeV^, expressed 
as ep cross sections, after experimental corrections and background subtraction, for four ranges 
in and in the W domains defined in Table HI The solid curves show the results of fits to the 
data in the mass range 0.6 < m^rTr < 1-1 GeV of the relativistic Breit-Wigner function with 
momentum dependent width defined in Eqs. (II im2l) . with skewing of the mass distribution 
following the Soding parameterisation given by Eq. (fT3l ); the dashed curves correspond to a 
non-skewed relativistic Breit-Wigner function and the dotted curves to the interference between 
resonant and non-resonant amplitudes. 
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Figure 10: dependence (a) of the Soding skewing parameter // defined in Eq. (fT3l) : (b) of tlie 
Ross-Stodolsky parameter n defined in Eq. for p elastic production. Measurements from 
HI [fTOl and ZEUS IITtI in photoproduction and E665 [[32ll in electroproduction are also shown. 
The present measurements are given in Table [T3l 



Figure \T0\ presents the dependence of the fitted values of the skewing parameters for 
elastic p production^, the mass and width of the resonance being fixed to the PDG values [[84l|. 
The skewing effects decrease with increasing Q"^, showing that the non-resonant amplitude 
decreases faster with than the resonant amplitude, as expected on theoretical grounds [|99ll . 
No significant dependence of the skewing parameters is observed as a function of or 

5.2.2 (p mesons 

The mass distribution for elastically produced kaon pairs is shown in Fig.dll after background 
subtraction and corrections for detector and QED radiation effects. It is described by the con- 
volution of the Breit-Wigner function defined by Eqs. (II 1H12I) with a Gaussian function of 
width a = 2 MeV describing the mass resolution, as evaluated using the MC simulation. The 
mass and width of the resonance, fitted over the interval 1.006 < uikk < 1.040 GeV, are 
1018.9 ± 0.2 (stat.) MeV and 3.1 ± 0.2 (stat.) MeV, respectively, reasonably close to the 
world average values of 1019.46 ± 0.02 MeV and 4.26 ± 0.04 MeV MM- Conversely, when the 
(p mass and width are fixed to the nominal values the fitted resolution, which is assumed to be 
Gaussian, is 1.0 ± 0.1 MeV. This value, which is slightly smaller than that obtained from simu- 
lations, is interpreted as to come from small systematic effects. As expected [98], no indication 
is found for skewing effects due to interference with non-resonant K^K^ production. 

5.3 dependence of the total cross sections 
5.3.1 Cross section measurements 

The measurements of the 7*p cross sections for p and meson elastic and proton dissociative 
production are presented in Fig. [121 as a function of the scaling variable {Q'^ + My). They are 

*The values of the parameters // and n slightly depend on the fit mass range. At low mass, this is related to 
the shape uncertainties reflected by the uncertainty in the parameterisation of Eq. O. For higher masses, the mass 
limit dependence may be due to additional interference of p mesons with heavier (p') resonances Ii98i . 
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Figure 11: Distribution of the ttikk rnass for elastic (f) production with \t\ < 0.5 GeV^, ex- 
pressed as ep cross section, after experimental corrections and background subtraction, for the 
and W domains defined in Tabled The solid curve shows the result of a fit to the data in the 
mass range 1.00 < rnxK < 104 GeV of a relativistic Breit-Wigner function with momentum 
dependent width defined in Eqs. (II 1H12I) . convoluted with the experimental resolution. 
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Figure 12: (Q'^ + My) dependence of the 7*p cross sections for W = 75 GeV: (a) p me- 
son production; (b) production. The upper points are for the elastic processes, the lower 
points for proton dissociative diffraction, divided by a factor 2 to improve the readability of 
the figures. Overall normalisation errors of 3.9% (4.6%) for elastic (proton dissociative) p pro- 
duction and 4.7% (5.3%) for production are not included in the error bars. ZEUS measure- 
ments lfT7] - [T9ll23ll24ll are also presented; when needed, they were translated to = 75 GeV 
using the measured W dependence. The superimposed curves are from the KMW model ||59l 
with GW saturation [73] (dash-dotted fines) and from the MRT model [45] with CTEQ6.5M 
PDFs [fTTll (dotted lines). The present measurements are given in Tables [T4lfT7l 
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quoted for W = 75 GeV using the W dependences parameterised as a function of following 
the measurements of section I5.5.2I Using the fits of the Q'^ dependence presented below, it 
is verified that the normalisations of the 1995 (SV) cross section measurement [4J and of the 
present measurement are in good agreement for p mesons (the ratio is 1.01 ± 0.10). For the </> 
data, the 1995 SV measurement is slightly lower than extrapolated from the present result (the 
ratio is 0.84 ± 0.11). This difference is attributed to the different treatments of the backgrounds. 
ZEUS measurements of p and electroproduction are also shown in Fig. [121 Whereas the p 
measurements agree well, (f) measurements of ZEUS are a factor 1.20 above the present data. 
When an improved estimation of the proton dissociative background, investigated for the latest 
ZEUS p production study |fT9ll , is used to subtract this background in their analysis, the 
cross section ratio of the two experiments is reduced to 1.06, which is within experimental 
errors ifTOOll . 
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Table 8: (Q'^ + My) dependence of the cross sections for p and (p elastic and proton dissociative 
production, parameterised in the form 1 / (Q"^ + My)"- , with (a) n constant and (b) n parame- 
terised as n = Cl ± C2 (Q^ + My). The 1995 (SV) measurements are normalised to those of 
1996-2000. 

The total cross sections roughly follow power laws of the type 1 / {Q^ ± My)" with values of 
n, fitted over the domain 1 < < 60 GeV2, given in Table [8t a). These values are compatible 
for elastic and proton dissociative scattering. They are also similar for p and mesons, which 
supports the relevance of the scaling variable (Q'^+My). 

The generally poor values of x^/d.ol. for fits with constant values of n confirms the ob- 
servation of [4J: compared to a simple power law, the cross section dependence is damped 
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for small values of (Q'^ + My) and steepens for larger values. An empirical parameterisation 
n = Ci + C2 (Q'^ + My) provides a significant improvement of the fit and a good description of 
the data (Table [8]^b)). It is interesting to note that the fitted values of the parameter ci are close 
to the value 2 expected in the Vector Meson Dominance model [87] for the exponent n when 



5.3.2 Comparison with models 

Predictions of the KMW dipole model [|59l with GW saturation [|73l are compared to the data 
in Fig. [121 The shape of the p elastic cross section measurement is well described. The nor- 
malisation of the prediction is low by 10%, while the overall normalisation error in the present 
measurement is of 4%. Predictions using CGC saturation [[741 (not shown) are nearly indistin- 
guishable, except for the highest bins in Q'^ where, however, the limited precision of the data 
does not allow to discriminate. The MRT model [|45l does not provide normalisation predic- 
tions, because of the uncertainty on the quark pair invariant mass window corresponding to the 
meson recombination. For this reason, the predictions for different PDF parameterisations are 
normalised to the data at = Q GeVl Both the CTEQ6.5M [1771 and the MRST 2004 NLO 
PDFs [[781 (not shown) lead to predictions which are compatible with the dependence of 
the data. It should however be noted that the normalisation factors required to fit the data are 
about 1.1 for CTEQ6.5M but larger than 2 for the MRST04 NLO PDF (see also UM)- This 
surprisingly large factor suggests that the gluon contribution in the MRST04 NLO PDFs is 
underestimated. 

For elastic (p production, the KMW predictions describe the shape of the distribution well, 
but are higher than the data by 25%. The MRT model gives a good description of the Q"^ 
dependence of the cross section, with normalisation factors similar to those for p mesons. 

5.3.3 Vector meson cross section ratios 

Figures [TSl a) and (b) present as a function of and (Q'^ + My), respectively, the ratio of the 
to p elastic cross sections, for which several uncertainties cancel, in particular those related to 
the subtraction of the proton dissociative background. The ratios are different because the same 
value of corresponds to different values of (Q'^ + My) for p and mesons, in view of the 
mass difference. A slight increase of the ratio with is observed for ^ 4 GeV^, whereas 
the ratio is consistent with being constant when plotted as a function of (Q'^ + My). Similar 
behaviours (not shown) are obtained for proton dissociative production. 

The cross section ratios, computed for the same domains in (Q'^ + My) for rho and phi 
mesons and for W = 75 GeV, are 



g2 ^0. 



(el.) 



0.191 ±0.007 (stat.) 



+0.008 
-0.006 



(syst.) ± 0.008 (norm.) 
(Q' + M^ > 2 GeV'), 



(p. diss.) 



0.178 ±0.015 (stat.) 



+0.007 
-0.010 



(syst.) ± 0.008 (norm.) 



(Q'±M2 > 3.5 GeV'), 



(14) 



35 



D 



0.3 



0.2 



0.1 



Elastic VM production 




0.3 



0.2 



0.1 



Elastic VIUI production 




10 

[GeV^] 



10 

Q^+Mv [GeV^] 



Figure 13: Ratio of the to p elastic production cross sections for W = 75 GeV: (a) as a 
function of Q^; (b) as a function of (Q'^+My). The overall normalisation errors on the ratios, 
which are not included in the error bars, are 4.0%. The present measurements are given in 
Table [H 



where the ratio of elastic cross sections includes the 1995 SV measurements {I < < 
2.5 GeV^). The measurements are close to the value expected from quark charge counting 
0/p = 2 : 9, but they tend to be slightly lower. 

Qualitatively, the behaviour of the ratio is consistent with the dipole model. At small Q^, 
the influence of the meson mass on the transverse size of the qq pair is larger, which implies 
that colour screening is expected to be larger for mesons than for p mesons. In contrast, for 
Q'^ :§> My, the transverse size of the dipole is given essentially by Q'^ and symmetry is expected 
to be restored. 

The dipole size effect also explains the strong increase with Q"^ of the J/tptop ratio, scaled 
according to the quark charge content J/V^ : p = 8 : 9, as presented in Fig.[T4Ta). and the fact 
that the ratio is nearly constant and close to unity when studied as a function of (Q'^ + My), as 
shown in Fig.[T4lfb) (note the different vertical scales). 

Although striking, the agreement with SU(4) universality is however only qualitative, with 
the scaled to p cross section ratios slightly below 1 and the scaled J/ip to p ratios slightly 
above 1. Scaling factors obtained from the VM decay widths into electrons [|43ll44l[65l are 
expected to encompass wave function and soft effects; the use of the factors given in ll65l 
modifies the scaled to p ratio very little and brings the scaled J/ipio p ratio slightly below 1. 



5.4 dependence of the polarised cross sections 

The separate study of the polarised (longitudinal and transverse) cross sections sheds light on 
the dynamics of the process and on the dependence of the total cross section. Soft physics 
contributions, related to large transverse dipoles, are predicted to play a significant role in trans- 
verse cross sections, whereas hard features should be significant in longitudinal amplitudes. At 
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relatively low values of the scale, {Q'^+My)/A < 3 GeV^, soft, "finite size" effects are however 
expected to also affect longitudinal cross sections. 

The extraction of the polarised cross sections presented in this section implies the use of 
the measurement of the cross section ratio R = ai/aT, which is performed using angular 
distributions and is discussed in section [63] 



5.4.1 Cross section measurements 

The total 7*p cross section can be expressed as the sum of the contributions of transversely and 
longitudinally polarised virtual photons: 

(Ttotir +p-^V + Y) = (TT + eaL = (TT{l + eR), (15) 

where e is the photon polarisation parameter, e ~ (1 —?/ + y^/2), with 0.91 <e < 1.00 
and (e) = 0.98 in the kinematic domain corresponding to the present measurement. 

The polarised cross sections, obtained from the measurements of the total cross sections 
and of R, with the value of e for the relevant Q"^, are presented in Fig. [15] for elastic p and 
production, as a function of (Q'^ + My). 

Results of power law fits with constant exponents are presented in Table |9] (the fit quality 
does not improve with a (Q'^ + My) dependent value of n). The fit values differ from the results 
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Figure 14: Ratios of u, (p and J/ip to p elastic production cross sections, scaled according to 
the quark charge contents, p : u : (j) : J/^/j = 9 : 1 : 2 : 8, plotted as a function of (a) Q'^; (b) 
(Q^ + My). The p cross section has been parameterised as described in Table [8tb). The ratios 
are determined for the HI (p (this analysis) and J/tJj [11211 measurements, and from the ZEUS 
p im, Lu Hal, </> and J/^Ij [|25ll26il studies. 



37 




Figure 15: (Q'^ + My) dependence of (a) longitudinal and (b) transverse 7*p cross sections for 
elastic p and meson production with W = 75 GeV. Overall normalisation errors of 3.9% for 
p and 4.6% for mesons are not included in the error bars. The superimposed curves are model 
predictions: GK 062 (shaded bands), MPS [162J (solid lines), INS with large wave function [[651 
(dashed lines), MRT [|45l with CTEQ6.5M PDFs [Ha and the same normalisation as in Fig.[l2] 
(dotted lines) and KMW [|59l with GW saturation [|73l (dash-dotted lines). The measurements 
are given in Tables [T9ll20l 
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Table 9: (Q'^ + My) dependence of the longitudinal and transverse cross sections for p and 
meson elastic production, parameterised in the form l/((5^ + M^)" with n constant. 



38 



n = 3 for the longitudinal and n = 4 for the transverse cross sections, obtained from a LO 
calculation of two gluon exchange P2ll . 

Model predictions for ctl and ctt are compared to the data in Fig. \T5\ The GPD predictions 
of the GK model [IMI are slightly too flat, both for and for ax, but the global normalisations 
are within the theoretical and experimental errors, which suggests that higher order effects, not 
included in the model, are weak. The KMW model [.591 describes well the shapes of the ai and 
o"t measurement and the absolute normalisation of ai, whereas the normalisation is too low for 
(Tt; this is the reflection of the good description of the shape for atot and of the prediction for 
R which is systematically too high (see Fig. [37] in section [63l) . The MRT |[45l predictions for 
the p polarised cross sections are reasonable, but for (p production they are too low for aL and 
too high for ar, which reflects the fact that the predictions for R are too low (Fig. [37l) . The 
INS fcf-unintegrated model with the compact wave function [65J gives predictions which are 
significantly too high both for aL and for (t^, and too steep for cjt (not shown); the predictions 
with the large wave function have better absolute predictions but are too steep for and for 
ax- The MPS dipole saturation model [|621 describes the data rather well. 




Figure 16: (Q'^ + My) dependences of the cross sections for p and (p elastic production 
with W = 75 GeV: (a) longitudinal cross sections, multiplied by the scaling factor (Q^ + 
MyY/Q'^; (b) transverse cross sections, multiplied by (Q^ + MyY/My. The superimposed 
model predictions are the same as in Fig. [151 



The same data and model predictions are presented in Fig. [TH where the longitudinal cross 
sections are divided by and the transverse cross sections by My, all being in addition mul- 
tiplied by the scaling factors {Q'^ + MyY to remove trivial kinematic dependences HlOlll . The 
breaking of the formal expectations (n = 3, n = 4) for the l/((5^ + My)" dependence of the 
longitudinal and transverse cross sections is manifest in this presentation. This is expected from 
the fast increase with of the gluon density at small x. Note that the cross sections in Fig. [16] 
are given for a fixed value of W and thus correspond to different values of x. The increase 
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with of the scaled longitudinal cross section is slower than that of the scaled transverse cross 
section. This is reflected in the Q'^ dependence of the cross section ratio R = ai/cT, which is 
slower than Q"^ /My (see section [631 Figs.[37]and[38l). 



5.4.2 Vector meson polarised cross section ratios 



Figure [T7] shows the to p and J/;/^ to p polarised cross section ratios, scaled according to 
the quark charge content of the VM {J /if) longitudinal cross sections are affected by very large 
errors due to the measurement errors on R and are not shown). 
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Figure 17: Polarised cross sections for the elastic production of (present measurements) and 
J li^ [11211 mesons, divided by the parameterisations of the p elastic polarised cross sections and 
scaled according to the quark charge contents, p : : J/^ = 9 : 2 : 8; in (a) longitudinal; (b) 
transverse cross sections. 



The ratios of the to p polarised production cross sections are within uncertainties indepen- 
dent of {Q'^ + My) and close to the ratio of the total cross sections (Fig. [141) . suggesting little 
effect of the wave functions. In contrast, the ratios of the J /^jJ to p transverse cross sections 
are very different from 1. This is because the polarisation states for p and mesons on the one 
hand and for J/ %jj mesons on the other hand are very different for the same iQ^ + My) value, in 
view of the dependence of R. The fact that the cross section ratios are consistent with being 
independent of (Q'^ + My) thus indicates that, within the present errors, no large difference is 
found between the small dipoles involved in transverse J/^ production and the dipoles involved 
in transverse p production, for (Q'^ + My) > 10 GeV^. 
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Figure 18: W dependence of the 7*]? cross sections for elastic (a)-(c) and proton dissocia- 
tive (b)-(d) production for several values of Q^: (a)-(b) p meson production; (c)-(d) produc- 
tion. The overall normalisation errors, not included in the error bars, are the same as in Fig. [I2l 
The lines are the results of power law fits. The present measurements are given in Tables [2T1I241 
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5.5 W dependences 



5.5.1 Cross section measurements 

Figure [T8] displays the W dependence of the 7*p cross sections for the production of p and 
mesons, for several values of Q'^. For the first time, measurements are performed for both the 
elastic and the proton dissociative channels. 

The W dependence of the cross sections is well described by power laws of the form 

a{Y +P + Y) (xW\ (16) 

represented by the straight lines in Fig. [TSl This parameterisation is inspired by the Regge 
description of hadron interactions at high energy, with 

5{t) = 4(ajp(t)-l), (17) 
ajp(t) = ajp(O) + a' ■ t. (18) 

In hadron interactions, typical values for the intercept and the slope of the pomeron trajectory 
are ajp(O) = 1.08 to 1.11 lUOlll and a' = 0.25 GeV^^ MM, respectively. 

5.5.2 Hardening of tlie W distributions witli 

The W dependence of the cross sections is presented in Fig. [19] in the form of the intercept of 
the effective pomeron trajectory, Q!jp(0), to allow comparison between different channels with 
different t dependences. The values of a;jp(0) are calculated for the present p and (p meson 
production from the W dependences following Eqs. (|16H18I) . using the measured values of (t) 
and the measurements of a' for p production given in Table \TOi the latter are derived from the 
evolution with t of the W dependence of the cross section. The measurements of a;jp(0) are 
presented as a function of the scale = (Q^ + My)/4 for p, </> and J/^p production, and as a 
function of p^ = for DVCS, as expected for the LO process. 

Up to (Q'^ + My)/A values of the order of 3 GeV^, the W dependence of the elastic cross 
section for both p and (p production is slightly harder than the soft behaviour characteristic of 
hadron interactions and photoproduction (Fig.[T9ta)). For the higher (Q'^ + My) /4 range, higher 
values of a;jp(0) are reached, of the order of 1.2 to 1.3, compatible with J/ ^/^ measurements. This 
evolution is related to the hardening of the gluon distribution with the scale of the interaction. 
Consistent results are obtained in the proton dissociative channel, but with larger uncertainties 
(Fig.Iiab)). 

5.5.3 Comparison with models 

In principle, the W dependence of VM production can put constraints on gluon distributions, in- 
cluding effects like saturation at very low x and large W values. All models predict a hardening 
of the W distribution with increasing Q^, following from the steepening of the gluon distri- 
butions. As examples, predictions are given in Fig. [20] for the GK GPD model I16I1I . the INS 
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Figure 19: Evolution with the scale /i^ = ((5^ + My-)/4 of the intercept of the effective pomeron 
trajectory, ajp(O), for p and production: (a) elastic production; (b) proton dissociation. HI 
measurements of DVCS [[H and J/ijj production [|T2| and ZEUS measurements of p [|T7l[T9l 
(for the low points, the value of a' in (T7] is used ), ^24] and J/ip production Il25ll26l are 
also shown. For DVCS, the scale is taken as p? = Q"^. The values 1.08 and 1.11 [102l|, typical 
for soft diffraction, are indicated by the dotted lines. The present measurements are given in 
Table |25l 



/cj-unintegrated model with the large wave function [65] and the KMW dipole with GW 
saturation [17311 . The MPS saturation model [62] (not shown) gives predictions for p production 
nearly identical to those of KMW. In general, relatively small differences are found between 
the model predictions for the W dependence, and the present data do not provide significant 
discrimination. Differences in normalisation between models in Fig. [20] reflect differences in 
the predicted Q"^ dependence of the cross sections. 

5.6 t dependences 

5.6.1 Cross section measurements 

The differential cross sections as a function of |t| for p and </) elastic and proton dissociative 
production are presented in Fig. [21] for different ranges in Q"^. They are well described by 
empirical exponential laws of the type da/dt oc e^^ 1*1 

The slope parameters h extracted from exponential fits in the range \t\ < 0.5 GeV^ for 
elastic scattering and |t| < 3 GeV^ for proton dissociation are presented in Fig. [22] as a function 
of the scale p"^ = (Q^ + My) j 4. The measurements of the proton dissociative slopes are the first 
precise determination at HERA for light VM in electroproduction; they constitute an important 
ingredient for the extraction of the elastic h slopes. In Fig.[22l p and measurements by ZEUS 
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Figure 20: Comparison with models of the W dependences of the 7*p cross sections given in 
Figs. fTSl a) and (c), for the elastic production of (a) p mesons; (b) mesons. The superimposed 
curves are model predictions: GK [|6T|| (shaded bands), INS with large wave function |[65l 
(dashed lines) and KMW |I59J with GW saturation [El (dash-dotted lines). 



and J/tp measurements are also presented as a function of (Q^ + My)/4, together with DVCS 
measurements (with fx^ = Q^). 

The present measurements of the b slopes for (Q^ + My)/4 < 5 GeV^ are higher than those 
of ZEUS L19J and also than those of a previous HI measurement JH. Two sources of systematic 
experimental differences are identified. The first is related to the estimation of the proton dis- 
sociative background, both in size and in shape. The subtraction of a smaller amount of proton 
dissociative background and the use of a steeper proton dissociative slope lead to shallower |t| 
distributions of the elastic cross section and to smaller b slope measurements. The use of a 
central value of 2.5 GeV^^ for the proton dissociative slope, as assumed in [jU, compared to 
the values measured here (Fig.[22tb)), leads to a decrease of the elastic slope determination by 
0.1 GeV^^, and a variation by ±20% of the amount of proton dissociative background induces 
a change of the elastic slope measurement by ±0.2 GeV~^ for Q'^ = 5 GeV^ and ±0.1 GeV~^ 
for Q'^ = 20 GeV^. The second - and major - source of discrepancy, for both VMs, is in the 
treatment of the to, and mostly p' backgrounds discussed in section I4.2.4[ Because of the 
non-detection of the decay photons, these backgrounds exhibit effective \t\ distributions which 
are much flatter than their genuine distributions and than the signal. Neglecting completely the 
presence of the p' background would lead in the present analysis to a decrease of the measure- 
ment of the elastic b slope by 0.4 GeV'^ for = 3 GeV^ and 0.2 GeV"^ for = 20 GeVl 



5.6.2 Universality of t slopes and hard diffraction 



In an optical model inspired approach, the t slopes for DVCS and VM production result from 
the sum of terms describing the form factors due to the transverse sizes of the scattered system 
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Figure 21: t dependence of the 7*p elastic (a)-(c) and proton dissociative (b)-(d) production 
cross sections for several values of Q^: (a)-(b) p production; (c)-(d) cj) production. Some distri- 
butions are multiplied by constant factors to improve the readability of the figures. The overall 
normalisation errors, not included in the error bars, are the same as in Fig. [121 The superim- 
posed curves correspond to exponential fits to the data (solide lines), to predictions from the 
MPS model [62J (dashed lines), and to fits of Eq. (l20l) parameterising the two-gluon form factor 
in the FS model [|54]| (dotted lines). The measurements are given in Tables [26M291 
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Figure 22: Evolution with the scale = {Q'^ + My)/A of the slope parameters b of the expo- 
nentially falling It I distributions of p and electroproduction: (a) elastic scattering; (b) proton 
dissociation. HI data for DVCS (H, p photoproduction ifTOl and J/ip production |[l2l[l4l| and 
ZEUS data for p ifTTlfTgll . Il23ll24l and J/^ If25ll26ll production are also presented. For DVCS, 
the scale is taken as p? = Q"^. The present measurements are given in Table [30l 



Y (by), of the qq dipole pair (bqq) and of the exchange (bp). An additional form factor reflecting 
the VM transverse size may also give a contribution, by, to the t slope for light VM production 
in models where the wave function plays an important role in the process, while being negligible 
for DVCS and for J/i^. The value of the slope can thus be decomposed as: 

b = bY + bgg + bip + by. (19) 

In elastic scattering, the slope by = bp reflects the colour distribution in the proton. For 
baryonic excited states with size larger than that of the proton, larger slopes (i.e. steeper t dis- 
tributions) than for elastic scattering may be expected. In contrast, when the proton is disrupted 
in the diffractive scattering, no form factor arises from the Y system and by is expected to be 
~ 0. The b]p contribution of the exchange is generally believed to be small and independent of 
Q^. There is indeed a priori no relation between and the transverse size of the exchange, at 
least for \t\ <^ and for as taken to be constant (LL BFKL). 

It is visible in Fig.[22]that, already for (Q^ + M^) /4 > 0.5 GeV^, the elastic b slopes for light 
VM electroproduction are significantly lower than in photoproduction, showing a departure 
from purely soft diffraction and a decrease of the relevant qq dipole transverse size. Until 
the scale (Q^ + My)/4 reaches values > 5 GeV^, light VM slopes are however significantly 
larger than for J/ip. This indicates the presence of dipoles with relatively large transverse sizes 
for light VMs in this domain. This is expected in the transverse amplitudes and also in 
longitudinal amplitudes until the fully hard regime is reached ("finite size" effects). Light VM 
and DVCS slopes are compatible when plotted as a function of the scales {Q'^+My)/A and Q'^, 
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respectively. For large scale values, they are consistent with the J /if) data, although they may 
be slightly higher. All these features confirm that the present Q'^ domain covers the transition 
from the regime where soft diffraction dominates light VM production to the regime where hard 
diffraction dominates. The comparable values of the slopes for p, (p and J /ip production in the 
harder regime suggests that light VM form factors are small. 

For proton dissociative diffraction, the t slopes shown in Fig.[22tb) have significantly smaller 
values than for elastic scattering. This is expected for Y systems above the nucleon resonance 
region, with vanishing values of 6y. The proton dissociative slopes for p and </> mesons are 
similar at the same {Q^ + My) / 4 value, but remain larger than for J / if), confirming the presence 
of large dipoles for {Q'^ + My)/A < 5 GeV^ or, alternatively, leaving room for a light VM form 
factor. 

5.6.3 Comparison with models 

In Figure [21] predictions of the MPS saturation model ll62l for the t dependence of the cross 
sections are shown, superimposed on the elastic measurements. The data fall faster with \t\ than 
predicted by the model, especially at small Q^. The discrepancy is particularly significant for 
production. 
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Figure 23: Dependence on the scale {Q'^ + My)/A of the parameter of the two-gluon form 
factor of the FS model [|54l . extracted from fits of Eq. (|20|) to the t distributions of p and 
elastic production cross sections. The J /if) measurements by HI in photoproduction lfT2]| and 
by ZEUS in electroproduction [|26ll are also presented. The present measurements are given in 
Table [HJ 

A dipole function with a t dependent two-gluon form factor has been proposed by Frankfurt 
and Strikman (FS) tMll, with 

da /At oc (1 + \t\/mlg)-\ (20) 

which tends to e^^'*' for t — )• 0, with h = A/ni^g. Fits of this parameterisation to the data for 
p and (j) elastic production in several bins in Q'^ are shown in Fig. [2T1 superimposed on the 
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measurements. The fit quality is good, similar to the exponential fits. Figure [23] presents the 
extracted values of the parameter as a function of (Q^ + My-)/4 for the p and (p elastic 
channels. The parameter increases with (Q'^ + My) / A, from about 0.6 GeV at 5 GeV^ to about 
0.8 GeV at 35 GeV^. A measurement in J/tp photoproduction is also shown. The (Q'^ + My) 
dependence of the form factor reflects the dependence of the t distributions, as summarised 
inFig.|22l 



5.6.4 Slope of the effective pomeron trajectory 
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Figure 24: W dependence of the 7*p cross sections for p meson production in four bins in \t\, 
for (a) g2 = 3.3 GeV^ and (b) = 8.6 GeV^. The lines are the resuks of power law fits. The 
notag (|t| < 0.5 GeV^) and tag (\t\ < 3 GeV^) samples are combined. The measurements are 
given in Table [32l 



The W dependences in four bins in |t| of the 7*p cross sections for p meson production are 
presented in Fig.[2lfor two values of Q^. The notag (\t\ < 0.5 GeV^) and tag (\t\ < 3 GeV^) 
samples are combined in order to extend the measurement lever arm in \t\. It was checked that, 
using only the notag events, compatible values of a' are obtained, although with much larger 
errors. The combination is also supported by the fact that the values of Q;jp(0) for the elastic 
and proton dissociative processes are compatible (see Fig.[T9l). 

The W dependences, which are observed to depend on \t\, are parameterised following 
the power law of Eq. (fT6l) . The extracted values of aip{t) = 5{t)/A + 1 are presented in 
Fig. [25] Linear fits to the t dependence of a]p{t), following Eq. ([T8l) . give the measurements 
of the slope a' of the effective pomeron trajectory reported in Table [T0[ Values slightly smaller 
than 0.25 GeV^^ and higher than are obtained. 
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Figure 25: t dependence of the values of the a]p{t) = S{t)/A + 1 parameters obtained from 
the linear fits to the W dependences shown in Fig. [24] (p production), for (a) = 3.3 GeV^; 
(b) Q'^ = 8.6 GeV^. The lines are the results of linear fits of the form of Eq. (fTSi) . The 
measurements are given in Table [33l 



Q' (GeV^) 


a' (GeV-^) 


3.3 
8.6 


0.19 ±0.07 
0.15 ± 0.09 tofe 



Table 10: Measurement of the slope of the effective pomeron trajectory a' for p production, 
from the \t\ evolution of the W dependence of the p cross section presented in Fig. [25l using 
Eqs. ([mill, for = 3.3 and 8.6 GeVl 



In soft diffraction, the non-zero value of the slope a' of the pomeron trajectory (a' ~ 
0.25 GeV^^) explains the shrinkage of the forward diffractive peak with increasing W: 

b = bo + 4a't \n{W/Wo) . (21) 

The parameter a' can thus in principle also be obtained from the evolution with W of the 
exponential \t\ slopes for elastic p production, but this measurement is affected by the large 
errors on b (not shown). 

Figure [26| summarises a' measurements by HI and ZEUS for DVCS and in photo- and 
electroproduction of p, (p and J/^ mesons. The a' measurement for p photoproduction [[TtI . 
which combines the ZEUS data at high energy with OMEGA results [J 04 J at low energy, is 
a' = 0.12 ± 0.04 GeV^^, which is lower than the value 0.25 typical for soft hadronic diffraction 
and is similar, within errors, to values of a' in electroproduction. Measurements of a' at large 
\t\ are consistent with 0, with small errors on the J/ip measurements [[T2l[25l[26l . 

In the BFKL description of hard scattering, the value of a', which reflects the average trans- 
verse momentum kt of partons along the diffractive ladder, is expected to be small. In Regge 
theory, the reggeon trajectories are fixed by the resonance positions, and slopes do not depend 
on Q"^. Evolutions of the effective pomeron trajectories with or |t| are thus an indication of 
additional effects, e.g. multiple exchanges and rescattering processes. 
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Figure 26: Slope of the effective pomeron trajectory a', presented as a function of the scale 
/i2 = (Q2 + M^)/4, together with measurements by HI and ZEUS lfT7l[T9ll24ti26ll for 

DVCS (upper limit 95% C.L., with the scale ^ = Q^) and p, and J/ip in photo- and elec- 
troproduction with |t| < 1.5 GeV^. The line a' = 0.25 GeV~^ represents a typical value in 
hadron-hadron interactions. 



5.7 Comparison of proton dissociative and elastic cross sections 

This section presents comparisons of the proton dissociative and elastic channels, for both p 
and (j) meson production. Measurements of the t integrated cross section ratios are first pre- 
sented, providing empirical information useful for experimental studies. The factorisation of 
VM production amplitudes into photon vertex and proton vertex contributions, which can be 
disentangled by comparing elastic and proton dissociative scatterings, is then discussed: the 
photon vertex contributions govern the dependence and the relative strength of the various 
helicity amplitudes, whereas proton vertex form factors govern the t dependence. Proton vertex 
factorisation ("Regge factorisation") has been observed to hold, within experimental uncertain- 
ties, for inclusive diffraction [11051 . Factorisation is tested here through the study of the 
independence of the VM production cross section ratios at t = and through the measurement 
of the difference bei. — bp^ ^j^^ between the elastic and the proton dissociative exponential t 
slopes. 



5.7.1 dependence of the cross section ratios 

Figure [27] presents, as a function of Q^, the ratio of the proton dissociative to elastic cross 
sections, for p and mesons. In the ratio, several systematic uncertainties cancel, in partic- 
ular those related to meson reconstruction. No significant dependence of the ratios on Q"^ is 
observed. 
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Figure 27: Q'^ dependence of the ratio of proton dissociative (My < 5 GeV) to elastic 7*p cross 
sections for W = 75 GeV: (a) p meson production; (b) (p production. The overall normalisation 
error on the ratios, which is not included in the error bars, is 2.4%. The measurements are given 
in Tables [Hand [351 



The average ratios of proton dissociative (with My < 5 GeV) to elastic cross sections, 
integrated over t, are: 

^My<5 GeV 

(p) = 0.56 ± 0.02 (stat.) tom (syst.) ± 0.01 (norm.) , 

^My<5 GeV 

tot,p. diss. ^ ^ +0.06 ^ Q (norm.) . (22) 

Ctoi,el. 

Within uncertainties, the values for the two VMs are compatible. Using the DIFFVM model to 
estimate the contributions of proton dissociative scattering with My > 5 GeV, the ratio of the 
proton dissociative cross section for the full My mass range to the elastic cross section is found 
to be close to 1. This value is used e.g. in 010511 . 



5.7.2 Cross section ratios for t = 

If the same object (e.g. a gluon ladder) is exchanged in proton dissociative and elastic scattering, 
proton vertex factorisation should be manifest through the independence of the cross section 
ratio for t = 0. 

For exponentially falling t distributions, the cross section ratio at t 
total cross sections and the b slopes as 

dcTp diss./*^*^ ^.j. g-^ (^tot,p. diss. ^p. diss. 

dcTei./di: <ytot,e\. K\. 

Figure [28]presents, as a function of Q"^, the cross section ratios at t = for p and produc- 
tion, as obtained from the total cross section ratios presented in Fig. [27] and the h slopes given 
in Fig.[22| 
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Figure 28: dependence of the ratio of the proton dissociative (with My < 5 GeV) to the 
elastic 7*p cross sections at t = and W = 75 GeV, -^-y^^5j-(t = 0), for (a) p meson 
production; (b) (p production. The overall normalisation errors, not included in the error bars, 
are the same as in Fig. [271 The measurements are given in Tables [34ll35l 



The average ratios for both VMs are measured as: 

j^My<5 GeV /^^ 

—2i^^^E—J—(t = 0)(p) = 0.159 ± 0.009 (stat.) tooll (syst-) ± 0.004 (norm.) , 

d(Jel./dt 

^^My<5 GeV 

P-"^"'- . / (t = O)(0) = 0.149 ± 0.021 (stat.) tooU (syst.) ± 0.003 (norm.) . 

d(Tel./dt 

(24) 



The ratios are observed to be independent of and consistent for the two VMs, which 
supports proton vertex factorisation. 

The ratios of the proton dissociative to elastic b slopes are also independent of Q^, with 
average values of 

V diss. I bei.ip) = 0.28 ± 0.01 (stat.) 1°;°^ (syst.) , 

bp. d^ss. / KiX(t>) = 0.27 ±0.05 (stat.) (25) 

This empirical observation is consistent with the independence of the total cross section 
ratios (Fig. [271) and of the cross section ratios at t = (Fig. [281). 



5.7.3 Difference in t slope between elastic and proton dissociative scattering 

In the optical model approach of Eq. ([T9l ). assuming pomeron universality, the difference be- 
tween the elastic and proton dissociative b slopes, bei. — bp, diss., is related only to the proton size 
and independent of the interaction scale at the photon vertex and of the VM species. 
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Figure 29: Slope differences bei. — bp, diss, between elastic and proton dissociative scattering for 
p and (j) meson production, as a function of (Q'^+My) / 4. Results of HI for DCVS JH and J/ ip 
photoproduction 012111411 and of ZEUS for p IfTTII and J/ [25','26] photo- and electroproduction 
are also shown. The present measurements are given in Table [361 



Figure [29] presents the slope difference b(,i — bp ^iss. for p and (p meson production, as a 
function of ((5^ + My)/4. Within errors, independent values for the slope differences are 
found, with consistent average values of 

bei. - bp. dissXp) = 5.31 ± 0.28 (stat.) (syst.) , 

bei. - bp. diss.i(l>) = 5.81 ± 1.14 (stat.) (syst.) . (26) 

These observations support proton vertex factorisation, with a proton form factor contribution 
of about 5.5 GeV~^ 

Measurements of J/4' photo- and electroproduction are also presented in Fig.[29l They are 
consistent with independence, with b^i. — bp. diss. = 3.50 ±0.07 GeV^^, a value significantly 
smaller than for p and (p production; for DVCS [[9l, the measurement is 3.88 ± 0.61 GeV^^. 
The difference observed between light and heavy VMs is difficult to understand in the optical 
model, since the contributions to the slopes of the qq dipole form factors and of possible VM 
form factors should cancel in the difference. It may indicate that the hard regime is not reached 
for p and (p mesons in the present kinematic domain. 



6 Polarisation Measurements 



Information on the spin and parity properties of the exchange and on the contribution of the 
various polarisation amplitudes are accessed in diffractive VM production through the distribu- 
tions of the angles 9, Lp and defined in Fig.[3l The present section presents, successively, the 
measurements of the spin density matrix elements, a discussion of the nature of the exchange, 
measurements of the longitudinal over transverse cross section ratio R, and measurements of 
the ratios and relative phases of the helicity amplitudes. The results are compared with QCD 
models. 
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6.1 Spin density matrix elements 
6.1.1 Measurement procedure 

In the formalism of Schilling and Wolf fT06'|, summarised in the Appendix, the angular distri- 
butions allow the measurement of spin density matrix elements given in the form r^f,, which are 
normalised bilinear combinations of the complex helicity amplitudes Ta^Aj^/.a^-Ajv' -^7 ^i^d Xy 
being the helicities of the virtual photon and of the VM, respectively, and Xn and Xn' those of 
the incoming proton and of the outgoing baryonic system Y. 

At HERA, the proton beam is not polarised and the helicity of the outgoing baryonic system 
Y is not measured; the helicities A ^ and A n' are thus integrated over. For the electron beam, 
transverse polarisation builds up progressively over the running period through the Sokolov- 
Temov effect but the related matrix elements are measurable only for ^ ml, where mg is 
the electron mass, and are not accessible in electroproduction. The electron beam is thus treated 
here as unpolarised. 

In these conditions, a total of 15 independent components of the spin density matrix re- 
main accessible to measurement. Under natural parity exchange (NPE) in the t channej^ five 
TxyX^i amplitudes are independent: two helicity conserving amplitudes (Too and Tn), two single 
helicity flip amplitudes (Toi and Tio) and one double flip amplitude (T_ii). 

The 15 matrix elements enter the normalised angular distribution W{9, Lp, 0) which is given 
in Eq. (|38l) of the Appendix. They are measured as projections of the W{9, ip, 0) distribution 
onto 15 orthogonal functions of the 0, Lp and angles, listed in Appendix C of [106]. In 
practice, each matrix element is given by the average value of the corresponding {6, Lp, cj)) 
function, calculated over the relevant data sample. For p production, the uj, cj) and p' background 
contributions to the angular distributions are subtracted following the results of the Monte Carlo 
simulations; no correction is performed for the interfering non-resonant txti channel but this 
is expected to have a small effect since the interference contribution is small, decreases with 
Q"^ and changes sign at the resonance mass value, so that it largely cancels when integrated 
over the selected mass range (see Fig. (9]). For production, the u, p' and dipion backgrounds 
are subtracted. Kinematic and angular distributions are corrected for detector acceptance and 
migration effects. The systematic errors on the measurements are estimated by varying the MC 
simulations according to the list given in Table |7] In addition, a systematic error related to the 
binning is assigned to the acceptance correction used for determining the average value of the 
projection functions; it is quantified by varying the number of bins in the 6', and angular 
variables. 

For both p and mesons, the matrix element measurements for the elastic and proton dis- 
sociative channels are found to be compatible within experimental errors. In order to improve 
the statistical significance of the measurements and to reach higher \t\ values, the notag and tag 
samples with \t\ < 0.5 GeV^ and \t\ < 3 GeV^, respectively, are combined. The large \t\ notag 
sample is not used because of the large p' background, as shown in Figs.|4tc)-(d). 
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Figure 30: Spin density matrix elements for the diffractive electroproduction of p mesons, as 
a function of Q^. The notag {\t\ < 0.5 GeV^) and tag {\t\ < 3 GeV^) samples are combined. 
ZEUS results fT9'| are also shown. Where appropriate, the dotted lines show the expected 
vanishing values of the matrix elements if only the SCHC amplitudes are non-zero. The shaded 
bands are predictions of the GK GPD model [|6T1l for the elements which are non-zero in the 
SCHC approximation; the curves are predictions of the INS /ct-unintegrated model [[651 for 
the compact (solid lines) and large (dashed lines) wave functions, respectively. The present 
measurements are given in Table [37l 
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Figure 31: Same as Fig.[30l for mesons. ZEUS results [[24l| for the r^^ matrix element are also 
shown. The present measurements are given in Table [38l 
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Figure 32: Spin density matrix elements for the diffractive electroproduction of p mesons, as a 
function of |t|, for two intervals in Q^: 2.5 < < 5 GeV^ and 5 < < QO GeV^ Where 
appropriate, the dotted lines show the expected vanishing values of the matrix elements if only 
the SCHC amplitudes are non-zero. The measurements are given in Table l42l 
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Figure 33: Same as Fig. [321 as a function of the mass mmr- The measurements are given in 
Tables m and SSI 
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6.1.2 Matrix element measurements 



The matrix element measurements are presented as a function of for p and (p production in 
Figs. [30] and [3T] and as a function of \t\ and the mass m^Tr for p production in two intervals of 
gMnFigs.[32]and[3l 
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Figure 34: Measurements, as a function of and of the p matrix element combinations 
Tqq + 2r\^ and Tqq + 2r\^, obtained from fits of Eq. (|4T]) to the </> angle distributions. The notag 
{\t\ < 0.5 GeV^) and tag < 3 GeV^) samples are combined. The dotted lines show the 
expected vanishing values of the matrix elements if only the SCHC amplitudes are non-zero. 
The measurements are given in Table l46l 



The present measurements as a function of and \t\ confirm with increased precision 
the previous HI results |l4l|5l[7]| and they are globally compatible with ZEUS measurements 
as a function of [[T9ll24l| . No significant dependence of the matrix elements with W is 
observed within the present data. Measurements (not shown) of the matrix elements Tqq and 
r°ip obtained from fits to the cos 6 and Lp distributions as given by Eqs. (I39H40I) of the Appendix, 
are in agreement with those presented in Figs. [30]to[33l For the combinations Tqq + 2r\^ and 
Tqq + 2r\^ for p mesons, measurements from fits of Eq. (|4T)) to the (j) distribution, which give 
smaller errors than the projection method, are presented in Fig. [34l 



6.1.3 Comparison with models 

Figures [301 and [311 present, superimposed on the p and (j) measurements, predictions of the GK 
GPD model [61] and of the INS /c^-unintegrated model ^65\ for two different wave functions; 
for the GK model, the SCHC approximation is used and only non-zero elements are shown. 

For p production (Fig. [30]), taking into account the experimental and theoretical uncertainties 
and the use of the SCHC approximation, the GK model [[6T]| gives a description of the data 
which is reasonable in shape but does not describe the normalisation well. The INS model |[65l 
reproduces the gross features of the evolution but there are problems in the details. The 
model with the compact wave function describes the Tqq matrix element evolution, but it fails 
for the other elements which are non-zero under SCHC {r\_i, Im r^_^. Re rfg, Im r^g); on the 

^NPE trajectories are defined as containing for t > Q poles with P = {—iy\P and J being the particle parity 
and spin, respectively. 
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other hand, the model with the large wave function gives a rather good description of these four 
elements, but fails badly for Tqq. In addi 
also in the regime with > 10 GeV^. 



elements, but fails badly for Tqq. In addition, both wave functions predict too low values for rf 



00' 



For mesons (Fig. [3TI) with less statistics, the picture is slightly different for the INS 
model ll65l . where the use of a large wave function gives a better description of all matrix 
elements, including Tqq, than the compact wave function. 



6.2 Nature of the exchange 
6.2.1 Natural parity exchange 



NPE,T 




[GeV^] 



NPE,T 




2 
1.5 
1 

0.5 




p 


2 


p 


; A H1 (|) 


-A H1 (|) 




1.5 

1 


T 


i 


1 


1 1 1 1 1 1 1 


0.5 



1 1 1 1 1 



10 



20 



[GeV] 



[GeV^] 



[GeV^] 



Figure 35: Asymmetry Pnpe.t between natural and unnatural parity exchange for transverse 
photons: (a)-(b) p mesons, as a function of and \t\; (c)-(d) mesons. The dotted lines 
indicate the value 1 expected for NPE. The measurements are given in Table l47l 



The observation at low energy [|29l[30l[34l of dominant natural parity exchange (NPE) sup- 
ports the attribution of the vacuum quantum numbers ( J^*" = 0^^) to the pomeron; the recent 
observation by the HERMES collaboration ll33l of the presence at low energy of a small con- 
tribution (about 6%) of unnatural parity exchange is attributed to quark exchange (tt, ai or bi 
exchange). At high energy, the modeling of diffraction as two gluon exchange implies a NPE 
character, in particular in the GK GPD model ||611 . 

With unpolarised beams and for a single value of the beam energies, the only accessible 
information about the parity of the exchange is the asymmetry Pnpe.t = (o"t ~ ) / + ) 
between natural (cr^) and unnatural (a^p) parity exchange for transverse photons, using Eq. (l43l) 
of the Appendix. Measurements of Pnpe,t as a function of and |t| for p and (p mesons are 
presented in Fig. |35l They are globally compatible with 1, which supports NPE for transverse 
photons. Natural parity exchange is assumed in the following. 



6.2.2 Helicity conserving amplitudes; SCHC approximation 



Inspection of Figs. |30] and |3T] shows that, for both p and meson electroproduction, the five 
matrix elements listed in Eq. (l44l) of the Appendix (r^Q, rj_^, Im rl_^. Re rf^, Im rfo), which 



60 



contain products of the two helicity conserving amplitudes, Tqo and Tn, are significantly differ- 
ent from zero, with the SCHC relations of Eq. (|45] ) being approximately satisfied. In addition, 
with the significant exception of Tqq, the other matrix elements are small or consistent with 0. 

In the present kinematic domain, SCHC is thus a reasonable approximation, which can be 
used to obtain information on the transition amplitudes. In order to decrease the sensitivity to 
the SCHC violating amplitudes, which increase with \t\ (see sections [6.2.3l and [64l) , only events 
with \t\ < 0.5 GeV^ are used in the rest of this section. 
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Figure 36: Cosine of the phase S between the Tqo and Tn helicity conserving amplitudes for p 
and production with |t| < 0.5 GeV^, measured as a function of from two-dimensional fits 
of Eq. (l47l) . in the SCHC approximation. The HERMES lf33l measurement on protons is also 
shown. The dotted line indicates the value 1 which corresponds to amplitudes in phase. The 
present measurements are given in Table |48| 

if) distributions; phase 6 between the SCHC amplitudes Under SCHC, the angular distri- 
bution W{9, if, (p) reduces to a function of the angles 9 and ip = (j) — ip, Eq. (l47l) . which allows 
the extraction in this approximation of the cross section ratio R = aL/cT and of the phase 5 
between the Too and Tn amplitudes. 

Measurements of cos 5 obtained from two-dimensional fits of Eq. (l47l) with R left free are 
presented in Fig. |36]as a function of for p and production (\t\ < 0.5 GeV^). They are in 
agreement with the measurements obtained with R fixed to the values measured in the SCHC 
approximation using the Tqq matrix element and Eq. (|28l) . 

The measurements of cos S are close to 1, indicating that the transverse and longitudinal 
amplitudes are nearly in phase. For p production with < 10 GeV^, cos 5 differs however 
significantly from 1, as is also observed for Q"^ around 2 GeV^ in the low energy measurement 
by HERMES [|33l . An indication of an increase of cos 5 toward 1 at high may be present in 
the data. An interpretation of a value of cos 5 different from 1 at high energy in terms of a 
dependence of aL far will be given in section [6!4.4[ 

6.2.3 Helicity flip amplitudes 

A significant violation of SCHC is observed in Figs. [30l and [311 through the non-zero value of 
the matrix element, for p and for mesons (see also Fig. [34lfor the Tqq + 2rf^ combination 
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measurement for p mesons). The Tqq matrix element is proportional to the product Re {TqqTI^) 
of Too, the leading SCHC amplitude, and Toi, the helicity flip amplitude describing the transition 
from a transverse photon to a longitudinal VM. In Figs. [32] and [33j non-zero values, with 
dependent strengths, are also observed in p production for the matrix elements Re rJo, Re 
and Im rl^, which contain the product Re {TiiTq-^) of Toi and the second SCHC amplitude Tn. 
The data tend to support the relation Im rfo = —Re r{o of Eq. (|49l ). Other matrix elements are, 
within errors, consistent with when integrated over t. 

These findings confirm the previous HI observation [|4l|71 in p production that the Toi he- 
licity flip amplitude is significantly different from in the present Q"^ domain and is dominant 
among the SCHC violating amplitudes, supporting the hierarchy (see for instance [|48jl) 

|Too| > |Tn| > |Toi| > |Tio| , |T_n|. (27) 

Note that helicity violation as such is not a signature for hard processes. When integrated over 
\t\, the Toi amplitude in the present kinematic domain is larger for low Q"^ than for large Q^, as 
shown by the Tqo matrix element measurement in Fig. [33l At low energy and for (Q^) around 
0.5 GeV^ the Tqi amplitude is non-zero, with |Toi| / ^jTooP + jTnP = 15 to 20% for W 
about 2.5 GeV and 11 to 14% for 10 < < 16 GeV CT. 

The Too matrix element increases with \t\, as observed in Fig. [32] (see also Fig. [34]) . This 
is expected on quite general grounds for helicity flip amplitudes, as will be discussed in sec- 
tion [63 



6.3 Cross section ratio R = ctl/ ctt 

The cross section ratio R = oi/aT is one of the most important observables in the study of 
light VM production since it is sensitive to the interaction dynamics, including effects related 
to the interacting dipole size or depending on the VM wave function. 

In the SCHC approximation, R can be calculated from the Tqq matrix element: 

p _ ^00 '00 /08^ 

iisCHC - 7fW - - 1 ^04 • 
-^11 ^ ^ '00 

In view of the observed violation of SCHC, a better approximation takes into account the 
dominant helicity flip amplitude Toi and uses in addition the measurement of Tqq: 

o ^ ^0^0 ^ 1 r',t-e{rl,)'+^{r'oty-Hr'ooy .... 

-^ll"r-^01 ^ ^'OO'T^VOOJ 

where NPE is assumed and the amplitudes are taken to be in phase. As expected, the effect 
of this improved approximation is mostly significant at large \t\ values, in view of the increase 
with \t\ of the helicity flip amplitudes: the corresponding measurement of R is lower than that 
obtained in the SCHC approximation by about 0.05 for \t\ =0.1 GeV^ and about 0.30 for |t| = 
1 GeV^, independently of Q^. Integrated over t, this makes a 7% difference. Measurements of 
R are presented in the following using the improved approximation of Eq. ([29]) . The general 
features of the kinematic variable dependences discussed below are similar when the SCHC 
approximation of Eq. ([28]) is used. 
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Figure 37: dependence of the ratio R = a^/aT of the longitudinal to transverse cross sec- 
tions measured using Eq. (|29l) for (a) p meson production; (b) (p production. Measurements of 
R in the SCHC approximation, for p photoproduction by HI [3] and ZEUS [17J and for p and 
(j) electroproduction by ZEUS [|17[|19[l24l are also shown. The superimposed curves are predic- 
tions of the models of GK [[611 (shaded bands), INS Il65l with the compact (solid lines) and the 
large (dashed lines) wave functions, MRT with the CTEQ6.5M PDF parameterisation [45] (dot- 
ted lines) and KMW [59| (dash-dotted lines). The present measurements are given in Table |49l 



6.3.1 dependence 

The measurements of R presented in Fig. [37] show a strong increase with Q^, which is tamed at 
large Q^, a feature already noted in previous HI [j4| and ZEUS HI 91 publications. 

For p production, the GK GPD model the MRT model [45] and the INS model [[651 
with the compact wave function give a good description of the measurements, whereas the 
KMW [59] predictions are too high and the INS model with the large wave function is ruled 
out. The predictions of the MPS model [62J (not shown) are very similar to those of KMW 
up to 10 GeV^, and then slightly lower. The dependence of the IK [|48l model (not shown) 
is similar to that of the MRT model, since it is derived in a similar way. For production, the 
KMW model gives a good description while the MRT predictions are too low; within the quoted 
uncertainty, the GK model describes the data; for the INS model, the large wave function gives 
a slightly better description than the compact wave function; the predictions of the MPS model 
(not shown) are again similar to those of KMW, although slightly higher at low Q^. 

R measurements for p, (p and J/'tp mesons are presented as a function of the scaling variable 
Q"^ /My in Fig.|38l The improved approximation, Eq. (|29|) . is used for the present data whereas 
the SCHC approximation is used for the other data, which makes little difference for the t 
integrated measurements. A smooth and common behaviour is observed for the three VMs over 
the full Q"^ /My range and the full energy range, from the fixed target experiments to the HERA 
collider measurements. 
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Figure 38: Ratio i? = (Tl/ctt as a function of the variable /My. Electroproduction mea- 
surements of p mesons by fixed target experiments (NMC [|3T1l . E665 [32] and HERMES L33J), 
of p and mesons by ZEUS [US El and of J/^ mesons by HI and ZEUS [[nilH are also 
shown. The dotted line represents the scaling behaviour R = Q'^ /My. 



The data are close to a law R = /My, represented by the dotted line, but they lie sys- 
tematically below the line, with a slower increase of R with increasing Q^. These features are 
easily understood in the MRT P31 and IK BSl models where the formal Q"^ /My evolution is 
damped by a factor 7^/(1 + 7)^ and the taming of the R evolution results from the decrease of 
7 with increasing Q'^. 



6.3.2 W dependence 

The W dependence of R is presented for p meson production in Fig. \39\a) for three intervals 
in Q^. Because of the strong correlation in detector acceptance between W and Q^, the lever 
arm in W for each domain in Q"^ is rather limited. As discussed in section [2l the onset of hard 
diffraction, characterised by a strong W dependence, is expected to be delayed for transverse 
amplitudes compared to longitudinal amplitudes. A harder W dependence is thus expected for 
(Tl than for ax, resulting in an increase of R with W. In view of the limited precision, no 
significant conclusion can be drawn from the present measurements. 
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Figure 39: Dependence, for p meson production, of the ratio R = ai/ of the longitudinal to 
the transverse cross sections, determined using Eq. ([291 ), on (a) W\ (b) \t\; (c) m^^, separately 
for 2.5 < g2 < 5 GeV^ and for 5 < < 59 QeV^ for W, the latter bin is divided into 
5 < g2 < 15.5 and 15.5 < < 60 GeV^. The curves in (c) are from the MRT model [fTOTll . 
The measurements are given in Tables [50ll52l 
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6.3.3 t dependence; — slope difference 

Figure \39\h) presents the measurement of i? as a function of \t\ for p mesons, in two bins 
in Q^. For exponentially falling t distributions, this can be translated into a measurement of 
the difference between the longitudinal and transverse t slopes, through the relation R(t) = 
(Jiit) /arit) oc e^*^^^^^^)!*!. Measurements of the slope difference — &t extracted from a fit 
of the t dependence of R are given in Table [TT] (for completeness, the result for production in 
one bin in t is also given in spite of the large errors). The errors are dominated by the systematic 
uncertainty on the p' background subtraction. A slight indication (1.5cr) is found for a negative 
value of 6l — 6t in the higher bin in Q"^. The use of the SCHC approximation of Eq. (l28l) 
instead of the improved approximation of Eq. (|29l) for the measurement of R does not affect the 
measurements of bi ~ bx- 



(GeV^) 


bL - br (GeV-2) 


p production 


3.3 
8.6 


-0.03 ± 0.27 

-0.65 ±0.U toil 


production 


5.3 


-0.16±0.56^?:f[; 



Table 11: Difference between the longitudinal and transverse slopes, bi — br, of the t distribu- 
tions for p (two bins in Q^) and meson production, calculated from the t dependence of the 
cross section ratio R = ai^jar obtained using Eq. (|29l ). 

A difference between the b slopes is expected to indicate a difference between the transverse 
size of the dominant dipoles for longitudinal and transverse amplitudes (see e.g. [65]). The 
indication for a negative value of 6l — 6t in the higher bin in Q'^ is consistent with the expectation 
that ol reaches a harder QCD regime than Of. Conversely, the absence of a \t\ dependence 
of R in the lower Q'^ range is consistent with the interpretation of b slope measurements in 
section 15.6.21 suggesting that large dipoles may be present in longitudinal amplitudes ("finite 
size" effects) for moderate values of the scale (C5^ + My)/4. 

6.3.4 m„„ dependences 

A striking decrease of the cross section ratio R with the increase of the m^^r niass, which was 
also reported by ZEUS [fT9ll , is observed in Fig.[39lc). This strong effect is not expected in cal- 
culations where the p meson is treated as a particle with well defined mass and wave function. 
A simple interpretation of the m^T^ dependence follows from the formal Q^/M^ dependence of 
the cross section ratio, if the mass M is understood as the dipion mass rather than the nominal 
resonance mass. Such an interpretation is in line with the open quark approach of the MRT 
parton-hadron duality model [45 J . and is qualitatively supported by the calculations superim- 
posed to the data in Fig.[39l;c) [[1071 . The mass dependence of R expected from the interference 
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Figure 40: Dependence of the exponential t slope for p elastic production as a function of the 
mass m^TT, for = 3.3 and 8.6 GeV^. The measurements are given in Table[ 



of resonant p and non-resonant nn production, discussed in [|99ll . is small compared to that 
observed here and should decrease with Q^. 

The b slopes of the \t \ distributions do not show any significant dependence on the mass (see 
Fig.|40l). which indicates that the m^j, dependence of R can not be explained by an hypothetic 
kinematic selection of dipoles with specific size, related either to transverse or longitudinal 
amplitudes. All this suggests that the VM wave function plays a limited role in the description 
of VM diffractive production. 



6.4 Helicity amplitude ratios and relative phases 

The measurements of the spin density matrix elements presented in Figs. [30] to [33] give access 
to the ratios and relative phases of the helicity amplitudes. Following the IK analysis B^ . 
four amplitude ratios, taken relative to the dominant Too amplitude, are measured from global 
fits to the 15 matrix element measurements, assuming NPE and taking all amplitudes as purely 
imaginary; negative values correspond to opposite phases. The measurements are presented in 
the following sections for p and mesons as a function of and \t\, and additionally for p 
mesons as a function of the m^^r invariant mass. The relative phases are then discussed. 



6.4.1 dependences 

The dependence of the four amplitude ratios for p and </> meson production are presented in 
Fig. [4TJ The strong decrease with Q'^ of the amplitude ratio Tu/Tqq for both VMs, which is 
consistent with a linear increase with 1/Q, is related to the increase of the cross section ratio 
R through the dominance of the SCHC amplitudes. For the first time, a dependence of 
the amplitude ratio Tqi/Tqo is also observed, for p meson production. This dependence is also 
visible in the comparison of the two Q'^ ranges in Figs. [42] and l43] No significant Q"^ dependence 
is observed for the amplitude ratios Tio/Tooand T^u/Tqo. 



67 



"''i 1 1 "''oo 





• H1 p 




1 






0.75 


• 

- • 




0.5 


• • 




0.25 


- 1 1 





20 

[GeV^] 



"^Ol I Tgo 



0.4 


-_• H1 p 




0.2 









b) 



20 

[GeV^] 



"^lo / "^oo 




[GeV^] 



'''-1 1 I "''oo 




[GeV^] 



"''i 1 I "''oo 



1 

0.75 
0.5 
0.25 



"''01 I "''00 



A H1 (]) 


0.4 


A H1 (j) 


▲ 






A 


0.2 




▲ 







e) 





20 
[GeV^] 



20 
[GeV^] 



0.2 
0.1 


-0.1 
-0.2 



"''10 ^ "''00 



: A HI (|) 












1 1 


9) 



20 
[GeV^] 



0.2 
0.1 

-0.1 
-0.2 



'''-1 1 I "''00 



E A HI (j) 




r I 




i 











20 

[GeV^] 



Figure 41: Ratios of the helicity amplitudes, calculated from global fits to the measurements of 
the 15 spin density matrix elements, as a function of Q'^: (a)-(d) p meson production; (e)-(h) cp 
production. NPE is assumed and all amplitudes are taken as purely imaginary. Where appropri- 
ate, the dotted lines show the expected null value of the ratio if the non-SCHC amplitudes are 
vanishing. The measurements are given in Table [54l 



In the IK ["481 model, the amplitude ratio Tu/Tqq is given by 

M 1 + 7 

Tn/Too = -r (30) 

Q 7 

where the decrease with Q"^ of the anomalous dimension 7 slows down the Q'^ evolution, and 
the amplitude ratio Tqi/Tqq is given by 

The model describes the Tn/Too evolution well for values of M = 0.6 GeV < nip and 7 = 0.7, 
or M = nip and 7 = 1.1 (not shown). The latter is preferred for the description of Tqi/Tqq, 
though the physical interpretation of this high value for the parameter 7 is unclear. 



6.4.2 \t\ dependences 



The t dependence of the amplitudes, empirically parameterised as exponentially falling, is 
mainly determined by the proton and VM form factors. It is a reasonable assumption that 
these form factors affect in a similar way all amplitudes, and that their effects cancel in matrix 
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elements and in amplitude ratios ||48l . The study of the t dependence of the amplitude ratios 
thus gives access, in the reaction dynamics, to features specific to the different amplitudes. 
Note, however, that this line of reasoning neglects the different t dependences for transverse 
and longitudinal amplitudes, related to different dipole sizes. 

Figure l42l shows, for both VMs, the \t\ dependences of the four amplitude ratios. For the 
first time, a decrease with \t\ of the ratio of amplitudes Tu/Tqq is observed, both for p and 
for (j) production (Figs. I42la) and (e)). The increase with \t\ of the normalised Tqi helicity flip 
amplitudes, which could be deduced from the behaviour of the Tqq matrix element, is confirmed 
in Figs.|42tb) and (f). For the second single flip amplitude, Tio, negative values with increased 
strength relative to Too are observed in Fig.|42tc) at large Q^. Finally, non-zero values are found 
in Fig.|42td) for the ratio of the double flip T_ii to the Too amplitude, with negative values of 
the ratio and intensity increasing with \t\ for both bins in Q'^. 

The \t\ dependence of the Tn to Too amplitude ratio, which is not predicted in the IK model, 
Eq. [30l may be understood as an indication of different transverse dipole sizes in transverse 
and longitudinal photon scattering, as discussed in section 16.3.31 for the t dependence of the 
cross section ratio R. This is substantiated by the calculation of the cross section ratio using the 
helicity amplitude ratios, the cross section ratio R = ai/ ctt being given by: 

o ^ 1 + 2 (Tio/Too)^ 

(Tn/Too)2 + (Toi/Too)2 + (T_n/Too)2- ^ ^ 

Following the procedure of section [6! 3. 3 [ the difference between the longitudinal and transverse 
slopes are extracted from the t dependence of R. The results are given in Table [l2l For p 
production, the same effect is observed as in Table [TTl where the value of R was obtained only 
from the measurements of the Tqq and Tqq matrix elements using Eq. (|29l ): a value of — bx 
consistent with for < 5 GeV^, and a negative value for > 5 GeV^. Errors are reduced 
due to the use of all amplitude ratios in the global fits, and the value of 6^ — 6t in the range 
with > 5 GeV^ is 3a away from 0. For production, the limited statistics do not allow to 
measure separately the slope difference in two bins in Q^. 



(g2) (GeV^) 


bi - br (GeV-^) 


p production 


3.3 
8.6 


-0.06 ± 0.22 toit 
-0.53 ±0.10 1^:57 


(p production 


5.3 


-0.70 ±0.23 toil 



Table 12: Difference between the longitudinal and transverse slopes of the t distributions for 
p (two bins in Q'^) and for (p meson production, calculated from the t dependence of the cross 
section ratio R = ai/cyT obtained using fits to the amplitude ratios, Eq. (|32l) . 

The t dependence of the helicity flip amplitudes for light quarks can be explained as follows. 
In the case of the Toi amplitude, the virtual photon with transverse polarisation fluctuates into a 
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quark and an antiquark which, given their opposite helicities, must be in an orbital momentum 
state with projection 1 onto the photon direction. During the hard interaction, the dipole size 
and the quark and antiquark helicities are unchanged, but a transverse momentum kt — ^/\t\ is 
transferred to the dipole, which modifies its line of flight and thus allows a change of the orbital 
momentum projection. The Tqi amplitude, which describes the production of a longitudinal 
meson from a transverse photon, is thus proportional to \/\t\. Similar reasons explain the 
t dependence of the Tio amplitude. Note that, at variance with the case of light VMs, for 
heavy VMs with a non-relativistic wave function — 2;~l/2), the exchange of orbital 

momentum cannot take place, thus implying SCHC. 

In the IK model the \t\ dependence of the single-flip to no-flip amplitude ratio Tqi /Too is 
given by Eq. (|3TI) . and that of Tiq/Too by 

M ^/\t\ V2 

Tio/Too = — , (33) 

Q 7 

respectively, where the negative value of the ratio is consistent with the p data in the higher Q'^ 
domain, Fig.l42lc'). 



In the two-gluon exchange picture of diffraction for the double flip Tn amplitude, the 
change by two units from the photon to the VM helicities requires in addition spin transfer by 
the exchanged gluons. The observation of a non-zero value for this amplitude may thus provide 
important information concerning gluon polarisation in the proton I1108II . The prediction of the 
IK model for T_ 11 /Too is 

T_ii/Too = v' + v\ (34) 
,0 = _ ^s\t\M 1 

with a dependence proportional to The model describes the \t\ dependence of the data, but 
the negative sign of T_ii/Too, both for < 5 GeV^ and > 5 GeV^, is at variance with the 
model expectation; this is attributed to the strong approximations involved in the parameterisa- 
tions mSL 



6.4.3 W and m„„ dependences 

No significant W dependence of the amplitude ratios is observed (not shown), which follows 
from the absence of a dependence of the matrix elements. The strong m^^r dependence of 
the ai/cT cross section ratio observed in Fig. [39lc) is confirmed in the ratio Tn/Too of the 
dominant SCHC amplitudes, as seen in Fig.|43l with a similar hint for Tqi/Tqq. As suggested 
in section [F.3. 41 these features may be related to the M/Q dependences in Eqs. (|30l) and (|3T1) . 
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Figure 42: Same as Fig. |4ll as a function of \t\: (a)-(d) p meson production, for two bins in 
Q2: 2.5 < g2 < 5 GeV^ (open circles) and 5 < < 59 GeV^ (closed circles); (e)-(h) 
production. The measurements are given in Table [55l 
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measurements are given in Table [56l 
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6.4.4 Amplitude relative phases 

In an extension of the fits performed in the previous sections, the phases between the amplitudes 
can be left free. To ensure proper convergence, the number of fitted quantities has to be reduced. 
In view of their small values, the approximation is made to put to the amplitudes Tiq and T_ii. 
When the phase difference cos(0oi — 0oo) is left free, it is pushed to the bound 1; it is therefore 
fixed to this value in the fii!l. 

The fitted phase difference cos(0ii — 0oo) is found to be systematically lower than 1, with 
the amplitude ratios Tu/Tqq and Toi/Tqo being compatible with those presented in the previous 
section. The average value of the phase difference for p mesons is 

cos(0ii - 0oo) = 0.936 ± 0.016 (stat.) (syst.), (37) 

which confirms the result of section 16.2.21 under the SCHC approximation, that the dominant 
longitudinal and transverse amplitudes are nearly but not completely in phase. 



7 Summary and Conclusions 

This paper reports on the measurement of diffractive p and (p meson electroproduction at high 
energy, both in the elastic and proton dissociative channels. The data were taken in the years 
1996 to 2000 with the HI detector at the ep collider HERA, in the kinematic domain 2.5 < 

g2 < 60 GeV^ 35<W < 180 GeV, \t\ < 3 GeV^ and My < 5 GeV. 

The total, longitudinal and transverse YP cross sections are measured as a function of the 
scaling variable Q'^ + My. They roughly follow power laws, and are well described by em- 
pirical parameterisations allowing the power to linearly depend on + My. The to p total 
cross section ratios are found to be independent of + My and consistent for elastic and 
proton dissociative scattering, with a value close to but slightly lower than the ratio expected 
from quark charge counting, : p = 2 : 9. The measurements significantly differ from the 
formal predictions n = 3 and n = A for the 1/(Q^ + My)" dependence of the longitudinal and 
transverse cross sections, respectively, which is attributed mainly to the increase with Q'^ of the 
gluon density at small x. 

The 7*j9 cross sections increase with the photon-proton centre of mass energy W , which is 
parameterised in the Regge inspired form oc W^, where 5 increases significantly with . This 
"hardening" of the W distribution is described in terms of the intercept aip(O) of the effective 
Regge trajectory. For values of the scale p^ = (Q^ + My)/4 up to about 3 GeV^, the W depen- 
dence of p and (j) production is slightly harder than the soft behaviour characteristic of hadron 
interactions and photoproduction, a;jp(0) = 1.08 to 1.11. For the higher {Q'^ + My)/A range, 
values of a]p{Q) of the order of 1.2 to 1.3 are reached, compatible with J /ij) measurements. 
DVCS measurements show a similar behaviour as a function of the scale p? = Q"^ . 

**The observation that cos((/)oi ~ 0oo) is close to 1 is at variance with calculations in pi09l|, where an attempt was 
made in a GPD approach to estimate the size of the Tqi amplitude within the handbag approach. The prediction 
in 01 091 that the amplitudes should be out of phase depends in fact strongly on a number of assumptions 111 101 . 
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The t dependences of the cross sections are well described as exponentially falling distri- 
butions oc e~''l*', up to \t\ values of 0.5 GcV^ for elastic production and 3 GcV^ for proton 
dissociation. The t slopes are measured for all four channels, providing the first precise de- 
termination at HERA of the proton dissociative slopes for light VM electroproduction. The 
values of the t slopes are lower than those in photoproduction and they decrease with increasing 
scale, in a way which is common to light VMs and DVCS. Values of the t slopes comparable 
to those for J/^' production, or slightly larger, are reached for a scale {Q'^ + My)/A > 5 GeV^, 
which suggests that light VM form factors are small and confirms that the dominant longitudi- 
nal amplitudes approach a perturbative behaviour for (Q'^ + My)/A around 3 to 5 GeV^. The 
correlation between the W and t dependences of the cross sections is parameterised in the form 
of the slope a' of the effective pomeron trajectory. For p meson production, this slope is smaller 
than that in soft hadron-hadron interactions, albeit with large errors. 

The ratio of the proton dissociative to elastic cross sections for |i| = and the difference 
between the elastic and proton dissociative slopes are measured to be independent of Q'^. These 
observations support the relevance of the factorisation of the process into a hard scattering 
contribution at the photon vertex and a soft diffractive scattering at the proton vertex ("Regge 
factorisation"). The value measured for p and production for the slope difference, bei. — 
bp. diss. — 5.5 GeV~^, however, is larger than for J/t/j production. 

Polarisation effects are studied through the measurement of 15 spin density matrix elements, 
which are normalised bilinear combinations of the complex helicity amplitudes T^yXj- The 
dependence on the kinematic variables and, for p mesons, on the dipion mass is measured. 
The main feature in the present domain is the dominance of the s-channel helicity conserving 
(SCHC) amplitudes. Too ^nd Tn, with Too > ^ii- In addition, a significant breaking of SCHC 
is manifest through the non-zero value of the Tqq matrix element, especially at large |t| values. 

The ratio R = ai/ ax oi the longitudinal to transverse cross sections increases strongly with 
Q"^, as predicted in pQCD, with a scaling behaviour as a function of Q"^ /My for the different 
VMs. The linear dependence it! = Q'^ / My predicted at LO, however, is damped for large 
values of Q^. No 1^ dependence of R is observed within errors. For t, an indication of the 
dependence of R is found for p meson production with > 5 GcV^. This can be interpreted 
as a difference between the longitudinal and transverse t slopes, hi — br, which differs from 
zero by 1.5a, with dominant systematic errors. A strong mTTTr dependence of R is observed 
for p meson production, both for smaller and larger than 5 GeV^. This behaviour may be 
interpreted as following from the general Q^/M^ dependence of VM production, if the mass M 
is understood as the dipion mass rather than the nominal resonance mass. 

The ratio of the helicity amplitudes is measured from global fits to the 15 matrix elements. 
Several features expected in pQCD are observed for the first time. A decrease with increasing 

is found for the amplitude ratio Toi/Too, which supports the higher twist nature of the he- 
licity flip amplitudes. The amplitude ratio Tn/Too is observed to decrease with increasing \t\, 
which may be related to different transverse sizes of transverse and longitudinal dipoles. This is 
substantiated by the non-zero value of the slope difference bi — &t obtained from the measure- 
ment of R from global fits of the helicity amplitudes, with a 3a significance. At large Q'^, the 
amplitude ratio Tiq/Tqq which involves the second single flip amplitude is found to exhibit a |t| 
dependence. Finally, a non-zero value at large \t\ is found for the ratio T_ii/Too which involves 
the double flip amplitude, an observation which may provide information on gluon polarisation 
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in the proton. The phase between the Too and Tu ampUtudes is measured to be non-zero, which 
may suggest different W dependences of the transverse and longitudinal amplitudes. 

The general features of the kinematic dependences of the cross sections and of the spin 
density matrix elements are understood qualitatively in QCD. In particular, the W and t depen- 
dences indicate that "hard", perturbative QCD features become dominant in the longitudinal 
cross section in the present kinematic domain, for (Q'^ + My)/A > 3 — 5 GeV^. The measure- 
ments are globally described by models using GPDs or a dipole approach, which differ in detail 
but agree on the gross features. 

The study of VM production at HERA thus provides new insights for the understanding of 
QCD and the interplay of soft and hard diffraction. 
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Appendix 



Matrix elements In the formalism of [II06II . the spin density matrix elements are normalised 
sums of products of two helicity amplitudes TxpX^,,x-,XN- They are given in the form r*^, where 
the notation of the upper index denotes the combination of unpolarised transverse and 
longitudinal photonfl the notations and are used for VM production by transverse photons 
with orthogonal linear polarisations, and and for the interference between VM production 
by transverse and longitudinal photons. The lower indices (j fc) refer to the VM helicities Ay of 
the pair of amplitudes. 



Angular distribution In the absence of longitudinal beam polarisation, 15 independent com- 
ponents of the spin density matrix can be measured (8 additional matrix elements are accessible 
with a longitudinally polarised lepton beam). They enter in the normalised angular distribution 

— V^Ker^Q sm26 cosif — r^ti sin^^ cos 2(^9 

— e cos 20 (^rl^ sin^ 6 + Tqq cos^ 9 — V2 Re r^Q sin 29 cos ip 

— rl_-^ sin^ 9 cos 2(p j 

— e sin 20 ^-\/2 Ini r^Q sin 2^ sin + Im r sin^ ^ sin 2(y9^ 
+ ^/2e {1 + e) cos rf-^ sin^ 9 + r^g cos^ 9 

— \p2 Re r^Q sin 29 cos — r\_Y sin^ 9 cos 2^9 ^ 
+ ^2e (1 + e) sin (72 Im r^Q sin 26' sin 99 

+Im r\_-^ sin^ ^ sin 2(^ j | . (38) 



Measurement of the matrix elements The matrix elements are measured as projections of 
the normalised angular distribution, Eq. (|38l) . onto orthogonal functions of the 9, ip and angles, 
with one specific function corresponding to each matrix element (see Appendix C of [ 106il ). In 
practice, each matrix element is measured as the average value of the corresponding function, 
taken over all events in the data sample. 

Alternatively, fits to the projections of the angular distribution W{9, y?, 0) onto each of the 
three angles provide measurements of the matrix elements Tqq and r\^^ and of the combinations 

^The separation of the and components is only possible through measurements with different polarisation 
parameters e, i.e. with different beam energies in the same detector configuration. In this case, 18 matrix elements 
in total can be measured. 
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(r^o + 2rfi) and (rl, + 2rl,y. 

W{e) oc 1 - r{]^ + (3 rg^ - 1) cos^ 6 (39) 

W{ip) oc 1 - 2r°ii cos (40) 

oc 1 + ^25(1 + e) cos (r^o + 2rfi) - e cos 20 (r^o + 2ri\). (41) 



Natural parity exchange Natural parity exchange (NPE) in the t channel implies the follow- 
ing relations between amplitude^ 

Av'AjY/ A^Aiv ( 1) ^ TxyXj^, ,X^Xff ■ (42) 

For unnatural parity exchange, an additional factor (—1) appears in the right hand side of 
Eq.dH. 

Under NPE and integrating over the nucleon polarisations, the number of independent TxyX^ 
amplitudes is reduced from 9 to 5: two helicity conserving amplitudes (Tqo and Tn = T_i_i), 
two single helicity flip amplitudes (Tqi = — To_i and Tio = — T_io) and one double flip ampli- 
tude (T_n = Ti_i). 

In general, longitudinally polarised lepton beams are required to separate natural and un- 
natural parity exchange process. However, unpolarised beams allow the measurement of the 
asymmetry Pnpe,t between natural (a^) and unnatural (a^) parity exchange for transverse 
photons: 

(jN _ U 

= Jn^Jj = 2 - + 2r?S - 2rl, - 2rU. (43) 

U rp -\- (J rp 

The measurement of the corresponding asymmetry for longitudinal photons requires different 
values of e, i.e. different beam energies. 



s-channel helicity conservation In the approximation of s-channel helicity conservation 
(SCHC) [1 1 IJ. the helicity of the virtual photon (measured in the helicity frame defined in sec- 
tion [3]4l) is retained by the final state VM (with the nucleon helicity also remaining unchanged). 
Single and double helicity flip amplitudes thus vanish (Tqi = Tiq = T_u = 0) and only five 
matrix elements are non-zero: 



0^, r|_i, Imr2_^, Re rfg, Imr^^; (44) 



Under SCHC and NPE, the following relations hold between these elements: 

r{_i = -Im r?_i = ^ ~ '^oo)' ^lo = "I™ ^lo- (45) 

In the case of SCHC, only two independent parameters are left, conveniently chosen as the 
cross section ratio R = ul/ot and the phase 5 between the Tqo and Tn amplitudes, with 

Toor;i = |roo||rn|e-^'^. (46) 



'"More precisely, Eq. ( |42| ) implies that, for |<| — |i|min> the trajectory exchanged in the t channel has natural 
parity. 
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The angular distribution W{9,ip, 0) then reduces to a function of 9 and ifj = (p — ip, the angle 
between the electron scattering plane and the p meson decay plane, in the 7*p frame: 

iy(cos^,V') = ^ ( sin^^e (1 + e cos2^/^) 

QTi \ + e R I 

+ 2eR cos^9- ^/2e{l + e)R cos 5 sin 2^ cos ^} . (47) 

In the SCHC approximation, the cross section ratio R is obtained from the measurement of 
the matrix element r^Q, as given by Eq. (|28l) . 

Dominant helicity flip amplitude Toi The precision of measurements performed in the 
SCHC approximation, especially at large \t\, can be improved by retaining the dominant helicity 
flip amplitude Tqi. Five additional matrix elements are then non-zero, supplementing the five 
elements given in Eq. (l44l) : 

Rer?^, r^o, Re rJo, Imr^g, r^Q. (48) 
Under NPE, the following relations hold in addition to the SCHC relations (l45l) : 

Re r% = -Re r^o = Im r^Q. (49) 

Assuming that the amplitudes are in phase, an improved approximation of the cross section 
ratio R is given by Eq. (|29l) , which uses the matrix elements Tqq and Tqq. 
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(GeV^) 


// 


n 


3.3 
6.5 
11.9 
23.0 


-0.470 ± 0.049 tofu 
-0.270 ± 0.059 to^olt 
-0.351 ± 0.095 +°:°^^ 
-0.100 ± 0.090 +°:°°^ 


1.73 ± 0.21 1°:°^ 
0.91 ± 0.24 tofs 
1.27 ± 0.39 +°:°? 
0.20 ± 0.37 



Table 13: dependence, for elastic p meson production, of the Soding skewing parameter // 
defined in Eq. (fT3l) and of the Ross-Stodolsky parameter n defined in Eq. 



(GeV^) 


(7(7*p — )■ pp) (nb) 


2.65 


563 


± 32 


+59 
-45 


2.95 


423 


± 24 


+49 
-34 


3.30 


383 


± 20 


+43 
-30 


3.70 


295 


± 17 


+33 
-25 


4.15 


232 


± 16 


+27 
-17 


4.65 


198 


± 13 


+20 
-19 


5.20 


154 


±9 


+ 17 

-8 


5.85 


131 


±9 


+ 12 

-7 


6.55 


102 


±8 


+8 
-7 


7.35 


79.1 


± 6.4 


+6.6 
-3.6 


8.20 


56.5 


±4.3 


+5.5 
-2.2 


9.20 


53.0 


±3.9 


+4.2 
-3.7 


10.3 


39.3 


±3.9 


+4.2 
-2.8 


11.5 


30.8 


± 2.9 


+2.8 
-2.5 


12.9 


25.8 


± 2.2 


+2.7 
-1.7 


14.5 


15.7 


± 1.3 


+1.7 
-1.2 


16.5 


12.5 


± 1.0 


+1.3 
-1.1 


18.8 


9.22 


±0.82 


+1.11 
-0.80 


21.7 


5.99 


±0.58 


+0.67 
-0.55 


25.0 


3.54 


± 0.42 


+0.39 
-0.31 


29.3 


2.24 


±0.31 


+0.31 
-0.28 


35.0 


1.68 


±0.27 


+0.24 
-0.21 


46.0 


0.742 


±0.105 


+0.101 
-0.088 



Table 14: dependence of the 7*p cross section for elastic p meson production for W = 
75 GeV. The overall normalisation uncertainty of 3.9% is not included in the systematic errors. 
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a(7*p —?■ pY) (nb) 


2 71 
91 


■^1 4 
217 


-1- '^9 
± 20 


+45 
-35 
+29 
-24 


3.82 


1 89 


-1- ifi 


+26 
-17 




lOl 




+17 
-16 


O.OKJ 


57 7 


zn u.z 


+6.8 
-6.0 




50.5 


zb 5.7 


+6.5 
-4.9 


7.60 


40.4 


±4.6 


+5.1 
-3.7 


9.30 
12.00 


25.9 
17.0 


± 2.3 
± 1.8 


+3.0 
-2.3 
+2.3 
-1.7 


14.85 


10.8 


± 1.4 


+1.2 
-1.2 


19.20 


3.39 


±0.45 


+0.53 
-0.56 


32.15 


1.01 


±0.12 


+0.16 
-0.16 



Table 15: Q'^ dependence of the 7*p cross section for proton dissociative p meson produc- 
tion for W — 75 GeV. The overall normalisation uncertainty of 4.6% is not included in the 
systematic errors. 



(GeV") 


(t(7*p — 7- (pp) (nb) 


2.71 


72.6 


±8.5 


+ 10.6 
-7.6 


3.21 


64.5 


± 6.2 


+8.3 
-6.0 


3.82 
4.52 


46.4 
35.0 


±4.5 
±4.0 


+5.2 
-3.6 
+4.1 
-3.4 


5.36 


25.1 


±2.7 


+2.6 
-1.4 


6.35 


18.2 


±2.0 


+1.9 
-1.5 


7.60 
9.30 


12.6 
7.04 


± 1.5 
±0.90 


+ 1.2 
-0.7 
+0.78 
-0.52 


12.00 


5.34 


±0.63 


+0.56 
-0.37 


14.85 


2.25 


±0.38 


+0.30 
-0.21 


19.20 


1.28 


±0.24 


+0.17 
-0.13 


32.15 


0.371 


± 0.076 


+0.019 
-I). 012 



Table 16: Q'^ dependence of the 7*p cross section for elastic meson production for W — 
75 GeV. The overall normalisation uncertainty of 4.7% is not included in the systematic errors. 
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(GeV^) 


(7(7*p ^ 0y) (nb) 


3.3 


33.4 ± 5.1 Hi 


6.6 


8.04 ± 0.89 tiM 


11.9 


2.66 ± 0.51 toil 


18.6 


0.779 ± 0.216 toill 


31.3 


0.273 ± 0.090 tomt 



Table 17: dependence of the 7*p cross section for proton dissociative p meson produc- 
tion for 1^ = 75 GeV. The overall normalisation uncertainty of 5.3% is not included in the 
systematic errors. 



(GeV^) 


(7(7*p (t>p)lcr{YP PP) 


Q2 + M2 (GeV^) 


(7(7*p (f>p)/cr{l*P PP) 


2.9 


0.148 


± 0.012 


+0.011 
-0.009 


3.94 


0.186 


± 0.015 


+0.014 
-0.012 


4.1 


0.171 


± 0.014 


+0.009 
-0.007 


5.14 


0.205 


± 0.017 


+0.011 
-0.008 


6.6 


0.169 


± 0.011 


+0.007 
-0.005 


7.64 


0.192 


± 0.013 


+0.008 
-0.005 


11.9 


0.171 


± 0.018 


+0.007 
-0.005 


12.94 


0.185 


±0.019 


+0.007 
-0.005 


18.6 


0.146 


± 0.025 


+0.009 
-0.006 


19.64 


0.154 


±0.027 


+0.010 
-0.007 


31.3 


0.187 


± 0.047 


+0.010 
-0.009 


32.34 


0.195 


± 0.049 


+0.010 
-0.009 



Table 18: Ratio of the to p elastic production cross sections for 1^ = 75 GeV, as a function 
of and of (g^±My). The overall normalisation uncertainty of 4.0% is not included in the 
systematic errors. 



(GcV") 


o'riY'P ~^ PP) (lib) 


o"l(7*P PP) (nb) 


2.9 


180 


± 11 


+23 

-17 


288 


± 26 tfo 


4.1 


78.1 


± 5.5 


+10.2 
-8.2 


165 


± 17 ^1,^ 


6.6 


23.7 


± 1.6 


+3.1 
-2.1 


74.7 


± 6.9 til 


11.9 


5.0 


±0.6 


+0.7 

-0.5 


24.0 


± 3.4 


18.6 


1.49 


±0.25 


+0.30 
-0.27 


7.7 


± 1.8 1°:? 


31.3 


0.27 


±0.06 


+0.05 
-0.05 


1.76 


±0.49 



Table 19: dependence of the transverse and longitudinal 7*p cross sections for elastic p 
meson production with W = 75 GeV. The overall normalisation uncertainty of 3.9% is not 
included in the systematic errors. 
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(GeV") 


(7r(7*P (t>P) (nb) 


(Jl(7*P 0P) (nb) 


3.3 
6.6 
15.8 


19.6 ± 2.4 til 
3.8 ± 0.5 t^l 
0.34 ± 0.10 1^6 


38.1 ± 7.3 
12.8 ± 2.6 til 
2.2 ± 1.0 til 



Table 20: dependence of the transverse and longitudinal 7*p cross sections for elastic 
meson production with 1^ = 75 GeV. The overall normalisation uncertainty of 4.7% is not 
included in the systematic errors. 



\Kje\ ) 


VV \\jiG V ) 


(t(7*P — )■ pp) (nb) 


3.3 
3.3 


41 
54 


308 
294 


± 14 
± 16 


+30 

-16 

+ ZD 

-17 


3.3 


67 


346 


± 19 


+30 
-21 


3.3 


80 


416 


±23 


+40 
-27 


3.3 


93 


397 


±27 


+40 
-25 


6.6 


48 


72.6 


±4.3 


+7.8 
-4.2 


6.6 
6.6 


64 
80 


97.2 
99.3 


± 5.9 
±6.2 


1 n A 
+y.4 

-5.5 

+7.2 

-5.0 


6.6 


96 


120 


±9 


+10 
-6 


6.6 


114 


115 


± 10 


+10 
-6 


11.9 


59 


25.3 


±2.6 


+2.4 
-1.3 


11.9 


77 


33.1 


± 2.7 


-1.5 


11.9 


95 


32.1 


±3.5 


+3.0 
-1.8 


11.9 


113 


27.5 


±4.1 


+3.5 
-1.9 


11.9 


131 


34.9 


±3.2 


+3.1 
-2.3 


19.5 


61 


6.5 


± 0.6 


+0.6 
-0.4 


19.5 


83 


9.6 


± 0.9 


+0.7 
-0.5 


19.5 


105 


9.3 


± 1.0 


+0.8 
-0.5 


19.5 


127 


9.8 


± 1.1 


+0.8 
-0.6 


19.5 


149 


16.9 


± 1.8 


+1.4 

-1.0 


35.6 


71 


1.2 


± 0.2 


+0.1 
-0.1 


35.6 


97 


2.0 


± 0.4 


+0.2 
-0.1 


35.6 


116 


2.3 


± 0.5 


+0.2 
-0.1 


35.6 


139 


3.4 


± 0.6 


+0.3 
-0.2 


35.6 


165 


2.9 


± 0.7 


+0.3 
-0.2 



Table 21: W dependence of the 7*p cross section for elastic p meson production, for several 
values of Q"^. The overall normalisation uncertainty of 3.9% is not included in the systematic 
errors. 
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(GeV") 


W (GeV) 


(7(7*p pY) (nb) 


3.3 
3.3 
3.3 
3.3 
3.3 


41 
54 
67 
80 
93 


177 ± 17 tfi 
209 ± 20 lif 
213 ± 21 tfo 
228 ± 26 
226 ± 33 tf^ 


7.5 
7.5 
7.5 

7.5 
7.5 


48 
64 
80 

96 
114 


36.6 ± 2.8 ttl 

38.8 ± 4.2 til 

34.9 ± 3.7 ttl 
40.2 ± 4.3 Its 
46.2 ± 5.2 til 


22.5 
22.5 
22.5 
22.5 
22.5 


71 
97 
116 
139 
165 


3.0 ± 0.4 

2.5 ± 0.5 tol 

3.6 ± 0.7 1°:^ 
4.6 ± 0.8 tol 
4.5 ±1.0 



Table 22: W dependence of the 7*p cross section for proton dissociative p meson production, 
for several values of Q^. The overall normalisation uncertainty of 4.6% is not included in the 
systematic errors. 
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(GeV^) 


W (GeV) 


(7(7*p (f)p) (nb) 


3.3 


41 


41.2 


±4.3 


+5.3 
-2.9 


3.3 


54 


55.1 


±5.7 


+5.4 
-3.7 


3.3 


67 


49.2 


±6.9 


+6.2 
-4.4 


3.3 


80 


57.5 


± 7.4 


+6.8 
-4.6 


3.3 


93 


69.6 


± 8.4 


+7.5 
-5.2 


6.6 


48 


13.2 


± 1.5 


+1.5 
-0.9 


6.6 


64 


13.0 


± 1.7 


+1.5 
-1.0 


6.6 


80 


20.5 


± 2.3 


+1.9 
-1.1 


6.6 


96 


14.7 


± 2.4 


+1.7 
-1.0 


6.6 


114 


23.3 


±4.0 


+2.2 
-1.6 


15.8 


71 


2.3 


±0.3 


+0.2 
-0.1 


15.8 


97 


2.5 


± 0.4 


+0.2 
-0.2 


15.8 


116 


3.9 


±0.6 


+0.3 
-0.2 


15.8 


139 


4.4 


±0.8 


+0.5 
-0.3 


15.8 


165 


7.6 


± 3.2 


+0.6 
-0.4 



Table 23: W dependence of the 7*p cross section for elastic meson production for several 
values of Q'^. The overall normalisation uncertainty of 4.7% is not included in the systematic 
errors. 



(GeV") 


W (GeV) 


(7(7*p ^ (f)Y) (nb) 


5.0 


50 


9.3 ± 1.6 t\i 


5.0 


70 


17 A ± 2.3 tli 


5.0 


90 


15.2 ± 2.5 tli 


5.0 


110 


11.6 ±2.9 tl'l 


5.0 


130 


32.7 ± 11.7 tli 


5.0 


150 


22.4 ± 13.7 



Table 24: W dependence of the 7*p cross section for proton dissociative meson produc- 
tion for Q"^ — 5 GeV^. The overall normalisation uncertainty of 5.3% is not included in the 
systematic errors. 
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(GeV^) 


6 


ap{0) 


7*p — 7- pp 


3.3 
6.6 
11.9 
19.5 
35.6 


0.40 ± 0.08 tofe 
0.57 ± 0.10 +°:°^ 
0.28 ±0.15 torn 
0.77 ±0.15 torn 
1.17 ±0.26 


1.10 ± 0.02 torn 
1.14 ±0.02 t'o'ol 
1.07 ± 0.04 tofi 
1.19 ±0.04 
1.29 ± 0.07 tofi 


'~f*p — i- (pp 


3.3 
6.6 
15.8 


0.53 ±0.17 t'^om 
0.52 ± 0.21 tofs 
1.09 ±0.34 1°:°^ 


1.13 ±0.04 f^of^ 

1.13 ±0.05 

1.27 ± 0.08 to^^ol 


yp py 


3.3 
7.5 
22.5 


0.32 ± 0.17 tofg 
0.17 ±0.14 t'Z 
0.58 ± 0.29 1°;}° 


1.17 ± 0.04 
1.12 ±0.04 
1.23 ± 0.07 tofr 


7*p — > 


5.0 


0.50 ± 0.24 toit 


1.20 ± 0.06 torn 



Table 25: dependence of the parameters 6 and a;jp(0), for elastic and proton dissociative p 
and meson production, computed from the W dependence of the cross section using Eqs. (fT6l - 
[TSl) . The values of a;jp(0) are obtained using the measured values of {t) and the measurements 
of a' for p production given in Table [TOl 
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Q- (GoV^) 


\f\ (Gov-) 


(la/dt p f>i 


) (iil)/GoV-) 


3.3 


0.025 


2156 


±82 


+202 
-141 


3.3 


0.075 


1379 


±73 


+132 
—90 


3.3 


0.125 


858 


±54 


+96 
-65 


3.3 


0.175 


665 


±48 


+72 
-60 


3.3 


0.250 


346 


±27 


+46 
-38 


3.3 


0.350 


234 


±25 


+35 
-28 


3.3 


0.450 


61.2 


± 13.9 


+18.8 
-16.4 


6.6 


0.025 


604 


±31 


+44 
-21 


6.6 


0.075 


392 


±24 


+35 
-16 


6.6 


0.125 


214 


± 18 


+19 
-12 


6.6 


0.175 


198 


± 16 


+15 
-9 


6.6 


0.250 


99.1 


±9.0 


+10.0 

-8.4 


6.6 


0.350 


50.0 


± 7.2 


+8.7 
-5.6 


6.6 


0.450 


31.5 


±5.5 


+5.8 
-4.7 


11.5 


0.025 


181 


± 18 


+14 

-10 


11.5 


0.075 


123 


± 12 


+11 
-8 


11.5 


0.125 


89.8 


± 10.6 


+9.5 
-5.4 


11.5 


0.175 


61.4 


±9.8 


+7.3 
-6.0 


11.5 


0.250 


44.0 


±8.0 


+5.0 
-4.5 


11.5 


0.350 


22.0 


±5.0 


+2.3 
-2.3 


11.5 


0.450 


8.58 


±3.02 


+2.06 
-1.35 


17.4 


0.025 


51.1 


±3.8 


+4.8 
-3.9 


17.4 


0.075 


35.6 


±3.5 


+3.6 
-2.8 


17.4 


0.125 


24.3 


± 2.9 


+2.6 
-2.2 


17.4 


0.175 


26.4 


± 3.4 


+2.7 

-2.4 


17.4 


0.250 


12.3 


± 1.4 


+2.0 
-1.2 


17.4 


0.350 


8.48 


± 1.40 


+1.16 

-1.04 


17.4 


0.450 


3.87 


±0.91 


+0.95 
-0.64 


33.0 


0.025 


6.52 


±0.93 


+0.87 
-0.77 


33.0 


0.075 


4.90 


±0.78 


+0.73 
-0.62 


33.0 


0.125 


4.43 


±0.71 


+0.56 
-0.50 


33.0 


0.175 


2.59 


±0.56 


+0.52 
-0.32 


33.0 


0.250 


2.28 


±0.39 


+0.28 
-0.28 


33.0 


0.350 


1.85 


±0.52 


+0.29 
-0.24 


33.0 


0.450 


0.660 


± 0.272 


+0.164 
-0.127 



Table 26: t dependences of the Yp cross section for elastic p meson production for several 
values of Q^. The overall normalisation uncertainty of 3.9% is not included in the systematic 
errors. 
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(GeV') 


\t\ (GeV') 

1 1 V ' 


dff/dt {Yp pY) (nb/GeV^) 




n inn 

\j . ±\j\j 


o 


± oz 


+47 

-48 




n snn 


01 /I 


± iy 


+28 
-19 




n 5nn 




± ib 


+21 

-15 




n 7nn 


7/1 1 


1 IIP 
± ii.o 


+11.1 
-8.0 


(J .(J 


n Qnn 


/ i.o 


-I- 1 Q 7 


+ 11.0 
-8.0 


3.3 


1.100 


9Q 1 


-1- a zL 


+6.0 
-4.3 


3.3 


1.300 


oU.o 


Zt iZ.4 


+6.6 
-4.0 


3.3 


1.500 


14. 


Zt 4.0 


+2.7 
-2.9 


3.3 


1.800 


o.DO 


-I- 9 9Q 

zt z.zy 


+1.78 
-1.78 


3.3 


2.250 


81^ 


-I- 1 Q9 
Zt ±.oZ 


+0.94 
-1.02 




2.750 


U.oU / 


± U.OOo 


+0,258 
-0. 10() 




100 


7R c: 


-I- 7 R 
ztz / .0 


+7.5 
-7.2 




300 


KQ 7 


zh /.D 


+8.8 
-5.2 




500 


OK n 
zo.U 


zh o.O 


+3.3 
-2.6 




700 

yj . 1 yjyj 


60 A 


zh 4.y 


+4.0 
-3.0 




QOO 

yj . tjyjyj 


1 Q 1 

io.i 


4- 1 O 

zh i.y 


+2.1 
-1.7 




1 100 


7 R"? 
/ .Do 


-1- 1 48 
zt 1.4(5 


+1.23 
-1.19 


6.6 


1.300 


u.yo 


zc i .o ( 


+ 1.30 

-0.86 


6.6 


1.500 


0. iZ 


4- 1 1 Q 
zh l.io 


+0.69 
-0.54 


6.6 


1.800 


o.Ui 


4- n R/i 

zh U.u4 


+0.59 
-0.52 


6.6 


2.250 


1 71 


4- n /I9 
zt yjAZ 


+0.39 

-0.36 


6 6 


2 75 n 


u.ozu 


4- n 978 
zt U.Z ( o 


+0.182 
-0.222 


i(J .(J 


n inn 


y.oo 


4-1 1 'J 
zt i.ij 


+ 1.32 
-1.33 


i(J .(J 


n '^nn 


4.00 


4- n R4 

zt 


+0.57 
-0.70 


15.8 


n.5nn 


/I fi V 

4.0 / 


±0.(0 


+0.71 
-0.56 


15 8 


n 7nn 


9 39 


+ n 35 


+0.33 
-0.33 


15.8 


n.9nn 


1 /I 


1 n on 
zt u.z / 


+0.28 
-0.24 


15 8 


1 inn 

± . ±\j\j 


1 8Q 


4- n 4^ 

zn yj.^j 


+0.23 
-0.35 


15.8 


i.snn 


0.882 


± n.2i6 


+0.225 
-0.093 


15.8 


1.500 


n.6i3 


± n.i93 


+0.093 
-0.611 


15.8 


1.800 


0.426 


± 0.108 


+0.096 
-0.095 


15.8 


2.250 


0.370 


± 0.089 


+0.095 
-0.089 


15.8 


2.750 


0.417 


± 0.245 


+0.037 
-1.058 



Table 27: t dependences of the 7*^ cross section for proton dissociative p meson production 
for several values of Q^. The overall normalisation uncertainty of 4.6% is not included in the 
systematic errors. 
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(GeV^) 




t (GeV^) 


da/dt (7*p ^ #) (nb/GeV^) 


3.3 


0.025 


431 


± 


34 


+40 
-25 


3.3 


0.075 


209 


± 


24 


+24 
-22 


3.3 


0.125 


120 


± 


21 


+24 

-16 


3.3 


0.175 


85.4 


± 


17.1 


+15.2 
-9.6 


3.3 


0.250 


64.6 


± 


10.5 


+9.0 
-7.4 


3.3 


0.350 


27.5 


± 


7.4 


+5.8 
-3.6 


3.3 


0.450 


27.4 


± 


7.1 


+5.4 
-4.5 


6.6 


0.025 


93.1 


± 


10.4 


+7.8 
-3.9 


6.6 


0.075 


77.7 


± 


9.1 


+6.8 
-3.4 


6.6 


0.125 


34.1 


± 


6.1 


+3.1 
-2.6 


6.6 


0.175 


24.9 


± 


5.0 


+2.8 
-1.7 


6.6 


0.250 


21.2 


± 


3.8 


+2.7 
-2.7 


6.6 


0.350 


8.77 


± 


2.37 


+1.54 
-1.00 


6.6 


0.450 


6.41 


± 


2.12 


+1.29 
-1.26 


15.8 


0.025 


8.24 


± 


1.14 


+1.09 
-0.89 


15.8 


0.075 


10.7 


± 


1.4 


+1.3 
-1.1 


15.8 


0.125 


3.89 


± 


0.85 


+0.59 
-0.42 


15.8 


0.175 


3.96 


± 


0.80 


+0.48 
-0.56 


15.8 


0.250 


2.32 


± 


0.47 


+0.39 
-0.29 


15.8 


0.350 


0.702 


± 


0.296 


+0.188 
-0.140 


15.8 


0.450 


0.349 


± 


0.278 


+0.160 
-0.150 



Table 28: t dependences of the 7*p cross section for elastic meson production for several 
values of Q'^. The overall normalisation uncertainty of 4.7% is not included in the systematic 
errors. 



(GeV") 




t (GeV^) 


dcr/dt (7*p ^ (f)Y) (nb/GeV^) 


5.0 


0.150 


58.2 ± 11.8 


5.0 


0.500 


23.1 ± 5.5 ttl 


5.0 


1.100 


6.17 ± 2.76 tlit 


5.0 


2.250 


0.681 ± 0.418 toill 



Table 29: t dependence of the 7*p cross section for proton dissociative meson production for 
Q"^ — 5 GeV^. The overall normalisation uncertainty of 5.3% is not included in the systematic 
errors. 
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(GeV^) 


b 


7*p — 7- pp 


3.3 
6.6 
11.5 
17.4 

33.0 


7.82 ± 0.33 tin 
7.57 ± 0.35 toll 
6.72 ± 0.53 toll 

5.86 ± 0.40 toll 

4.87 ± 0.66 toll 


7*]5 — i- (f)p 


3.3 
6.6 

15.8 


8.28 ± 0.80 tlil 
7.17 ± 0.73 toll 
7.08 ± 0.71 toll 


Yp — ^ pY 


3.3 
6.6 

15.8 


2.29 ± 0.12 tin 

1.91 ± 0.26 

1.70 ±0.15 tlfo 




fp (f)Y 


5.0 


2.21 ± 0.37 



Table 30: dependence of the b slope parameters of the exponentially falling \t\ distributions 
of p and (j) elastic and proton dissociative production. 



(GeV^) 


m2g (GeV) 


7*p — )■ pp 


3.3 
6.6 
11.5 
17.4 
33.0 


0.59 ± 0.01 toll 
0.60 ± 0.02 
0.65 ± 0.03 toll 
0.71 ± 0.03 toll 
0.80 ± 0.06 toll 


7*p — )■ (pp 


3.3 
6.6 
15.8 


0.57 ± 0.03 toll 
0.63 ± 0.04 +°:°^ 
0.64 ± 0.04 toll 



Table 3 1 : Parameter of the two-gluon form factor of the FS model [|54ll . extracted from fits 
of Eq. (|20l ) to the t distributions of p and elastic production cross sections. 
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t (GeV^) 


W (GeV) 


a{Yp pY) (nb) 


= 3.3 GeV^ 


0.08 


45 


1451 


±51 


+124 
-56 


U.Uo 


oo 


1 fi77 


-1- fifi 
zn oo 


+123 
-75 


0.08 


87 


2030 


± 82 


+173 
— 102 


0.32 


45 


393 


±24 


+45 
-29 


0.32 


65 


433 


±30 


+48 
-37 


0.32 


87 


m 


± 33 


— o5 
-ii 


0.69 


45 


95.1 


± 9.4 


+16.1 
-13.5 


0.69 


65 


110 


± 13 


+20 
-17 


0.69 


87 


124 


± 19 


+22 
-17 


1.45 
1.45 


45 
65 


14.5 
9.0 


±2.5 
± 1.8 


+3.2 
-2.7 
+2.7 
-2.4 


1.45 


87 


8.3 


±2.1 


+2.3 
-2.2 


g2 = 8.6 GeV^ 


0.08 


65 


261 


± 14 


+23 
-11 


0.08 


79 


261 


± 11 


+17 

-10 


0.08 


104 


320 


± 14 


+22 
-13 


0.32 


65 


64.8 


±5.9 


+6.4 
-3.8 


0.32 


79 


72.8 


±5.3 


+6.0 
-5.3 


0.32 


104 


82.1 


± 5.6 


+8.7 
-5.9 


0.69 


65 


16.5 


± 1.7 


+3.1 
-2.6 


0.69 


79 


17.5 


± 1.6 


+3.0 
-2.5 


0.69 


104 


17.3 


± 1.7 


+3.6 
-3.2 


1.47 


65 


2.8 


±0.5 


+0.8 
-0.8 


1.47 


79 


1.5 


±0.3 


+0.9 
-0.8 


1.47 


104 


2.3 


± 0.4 


+1.0 
-0.9 



Table 32: W dependence of the 7*p cross sections for p meson production in four bins in 

for g2 = 3.3 GeV^ and = 8.6 GeVl The notag {\t\ < 0.5 GeV=^) and tag (\t\ < 3 GeV^) 
samples are combined. The overall normalisation uncertainty of 4% is not included in the 
systematic errors. 
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t (GeV^) 


ajp(t) 


Q2= 3.3 GeV^ 


-0.08 
-0.32 
-0.69 
-1.45 


1.12 ± 0.02 tofi 
1.05 ±0.04 torn 
1.10 ±0.06 
0.76 ± 0.12 tofr 


Q2= 8.6 GeV^ 


-0.08 
-0.32 
-0.69 
-1.48 


1.12 ± 0.04 tofi 
1.12 ± 0.06 tofi 
1.02 ± 0.07 toZ 
0.93 ± 0.16 toll 



Table 33: t dependence of a]p{t) for p meson production, for two values of Q^. 



(GeV") 


ftot.p. diss. / \ 
0-tot,el. 


do-p. d 

dcTel 


ss. /dt f^. 

/At - 


0)(P) 


2.9 


0.58 


±0.04 


+0.12 
-0.08 


0.169 


± 0.017 


+0.011 
-0.015 


4.1 


0.65 


±0.05 


+0.07 
-0.06 


0.191 


± 0.019 


+0.014 
-0.016 


6.6 


0.53 


±0.03 


+0.02 
-0.06 


0.133 


± 0.021 


+0.009 
-0.011 


11.9 


0.58 


±0.05 


+0.09 
-0.07 


0.147 


± 0.022 


+0.041 
-0.016 


18.6 


0.45 


±0.05 


+0.04 
-0.05 


0.131 


± 0.021 


+0.042 
-0.019 


31.3 


0.57 


±0.09 


+0.04 

-0.06 


0.198 


± 0.044 


+0.062 
-0.030 



Table 34: dependences, for W = 75 GeV, of the ratios of proton dissociative (My < 
5 GeV) to elastic p meson production cross sections integrated over t and for t = 0. The overall 
normalisation uncertainty of 2.4% is not included in the systematic errors. 



(GeV^) 


f^toi.p. diss. /±\ 




3.3 
6.6 
15.8 


0.58 ±0.09 toZ 
0.48 ± 0.06 to% 
0.47 ± 0.08 torn 


0.155 ± 0.039 tomQ 
0.148 ± 0.034 tomt 
0.146 ± 0.038 to^^^ 



Table 35: dependences, for 1^ = 75 GeV, of the ratios of proton dissociative (My < 
5 GeV) to elastic meson production cross sections integrated over t and for t — Q. The overall 
normalisation uncertainty of 2.4% is not included in the systematic errors. 
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(GeV^) 


hei. - hp. diss. (GeV ^) 


P 


3.3 

6.6 
15.8 


5.52 ± 0.40 +^:26 

5.74 ± 0.62 

4.76 ± 0.48 to^^ 




5.0 


5.81 ± 1.15 to^l 



Table 36: Slope differences hei. — hp. diss, between elastic and proton dissociative scattering for 
p and (p meson production as a function of Q"^. 
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(g2) (GeV^) 


3.1 


4.8 


7.5 


12.2 


23.1 


00 





597 


± 





013 


+0.014 





680 


± 





017 


+0.017 

09S 





789 


± 





013 


+0.018 

n 09Q 

— u.uzy 





793 


± 





016 


+0.018 

09*^ 





877 


± 





019 


+0.019 

090 
— U.UZU 







049 


± 





009 


+0.010 

— U-UIU 





019 


± 





012 


+0.012 

01 
— U.UiO 





034 


± 





010 


+0.011 

010 


-0 


Oil 


± 





012 


+0.009 

OOQ 

— u-uuy 





046 


± 





014 


+0.008 

OOQ 

— u.uuy 


^04 
1 — 1 





000 


± 





Oil 


+0.004 
—0.004 





Oil 


± 





013 


+0.002 
—0.003 


-0 


003 


± 





010 


+0.002 
—0.002 


-0 


031 


± 





013 


+0.003 
—0.004 





012 


± 





014 


+0.002 
—0.004 


^00 


-0 


001 


± 





029 


+0.035 

019 





021 


± 





024 


+0.015 
n 09S 





001 


± 





039 


+0.020 

01 


-0 


081 


± 





095 


+0.016 

01 


-0 


015 


± 





061 


+0.026 

01 1 
— U.Ui 1 


^11 


-0 


019 


± 





031 


+0.006 

091 


-0 


034 


± 





036 


+0.010 
— u.uuo 


-0 


028 


± 





044 


+0.011 

01 1 
— U.Ui ± 





027 


± 





051 


+0.015 

01 
— U • U 1 





058 


± 





077 


+0.014 

01 fi 
— U.UiO 


Re tIq 


-0 


029 


± 





013 


+0.014 

010 


-0 


043 


± 





017 


+0.011 

019 


-0 


007 


± 





014 


+0.013 

010 


-0 


019 


± 





017 


+0.010 

OOQ 

— u-uuy 





023 


± 





021 


+0.013 

01 1 
— U.Ui 1 







157 


± 





015 


+0.008 
—0.008 





088 


± 





018 


+0.005 
—0.008 





117 


± 





014 


+0.006 
—0.007 





068 


± 





017 


+0.006 
—0.006 


-0 


019 


± 





021 


+0.008 
—0.009 







031 


± 





013 


+0.011 
-0.016 





033 


± 





016 


+0.010 
-0.009 


-0 


040 


± 





014 


+0.012 
-0.011 


-0 


024 


± 





016 


+0.007 
-0.008 


-0 


012 


± 





020 


+0.012 
-0.009 


Im rf_^ 


-0 


176 


± 





015 


+0.010 
-0.006 


-0 


133 


± 





018 


+0.005 
-0.007 


-0 


083 


± 





014 


+0.003 
-0.004 


-0 


045 


± 





016 


+0.003 
-0.005 


-0 


041 


± 





020 


+0.009 
-0.006 


'00 





156 


± 





019 


+0.040 
-0.065 





171 


± 





025 


+0.0.38 
-0.035 





080 


± 





022 


+0.040 
-0.041 





130 


± 





026 


+0.039 
-0.039 





135 


± 





033 


+0.032 
-0.034 




-0 


008 


± 





014 


+0.028 
-0.014 





Oil 


± 





017 


+0.008 
-0.014 





010 


± 





015 


+0.010 
-0.009 





001 


± 





018 


+0.009 
-0.010 





006 


± 





022 


+0.010 
-0.009 


Re rfo 





168 


± 





006 


+0.004 
-0.005 





141 


± 





008 


+0.004 
-0.005 





158 


± 





007 


+0.006 
-0.007 





128 


± 





008 


+0.005 
-0.005 





085 


± 





010 


+0.003 
-0.004 


^1-1 





001 


± 





008 


+0.002 
-0.002 





017 


± 





010 


+0.002 
-0.004 


-0 


009 


± 





007 


+0.003 
-0.003 


-0 


006 


± 





009 


+0.003 
-0.005 


-0 


003 


± 





010 


+0.007 
-0.004 


^6 
'10 


-0 


156 


± 





006 


+0.006 
-0.005 


-0 


141 


± 





007 


+0.005 
-0.004 


-0 


134 


± 





007 


+0.006 
-0.005 


-0 


117 


± 





008 


+0.004 
-0.004 


-0 


095 


± 





010 


+0.003 
-0.003 


^6 
' 1-1 


-0 


003 


± 





007 


+0.001 
-0.003 





003 


± 





009 


+0.004 
-0.002 


-0 


Oil 


± 





007 


+0.003 
-0.002 


-0 


035 


± 





008 


+0.004 
-0.004 





022 


± 





010 


+0.002 
-0.011 



Table 37: Spin density matrix elements for the diffractive electroproduction of p mesons, as a function of Q^. The notag (\t\ < 0.5 GeV^) 
and tag (\t\ < 3 GeV^) samples are combined. 



(g^) (GeV^) 


3.3 


6.6 


'00 


0.581 ± 


0.023 iHi^ 


0.746 ± 0.024 


Re r?o^ 


-0.004 ± 


0.017 


0.011 ± 0.018 


'1-1 


-0.059 ± 


020 +° °°^ 


0.047 ± 0.020 


'^00 


-0.060 ± 


n 1 74 +0-011 

^■^1^ -0.013 


-0.049 ± 0.070 


ril 


-0.059 ± 


n 1 c:q +0.008 
U.lOO _o,oiO 


0.006 ± 0.044 




-0.044 ± 


023 +° °2i 

U.UZO _o,oi5 


-0.073 ± 0.025 


rl-i 


0.220 ± 


027 +° °^^ 

U.UZ/ -0.013 


0.104 ± 0.029 


Mo 


-0.038 ± 


0.023 1°;°?^ 


0.075 ± 0.027 


Im r'l_i 


-0.152 ± 


0.028 


-0.111 ± 0.029 


'00 


0.053 ± 


0.034 


0.080 ± 0.040 




0.004 ± 


0.025 t',^ 


0.015 ± 0.028 


Rerfo 


0.220 ± 


r. n-, -, +0.009 
U.Ull _o,o08 


0.139 ±0.012 


rl-i 


-0.010 ± 


01 5 +° °°^ 

U.UiO _o.O05 


-0.002 ± 0.015 


^6 

Mo 


-0.147 ± 


01 +° °°^ 

U.UIU _o.O06 


-0.174 ± 0.012 


' 1-1 


-0.039 ± 


A Q-l o +0.006 
U.UIO _0.009 


0.006 ± 0.014 



+0.017 
-0.016 
+0.012 
-0.013 
+0.004 
-0.003 

+0.018 
-0.018 
+0.012 
-0.009 

+0.017 
-0.013 

+0.009 
-0.007 
+0.011 
-0.014 

+0.008 
-0.016 
+0.030 
-0.036 

+0.011 
-0.012 
+0.005 
-0.004 

+0.003 
-0.004 

+0.005 
-0.006 
+0.004 
-0.005 



15.8 



0.864 ± 
0.007 ± 
-0.020 ± 
-0.008 ± 
-0.004 ± 
-0.028 ± 
0.058 ± 
-0.017 ± 
-0.094 ± 
0.112 ± 
-0.010 ± 
0.091 ± 
0.035 ± 
-0.121 ± 
-0.003 ± 



0.031 
0.024 
0.025 
0.018 
0.014 
0.034 
0.036 
0.034 
0.034 
0.055 
0.037 
0.017 
0.019 
0.016 
0.017 



+0.019 
-0.016 
+0.014 
-0.011 

+0.006 
-0.002 

+0.020 
-0.026 

+0.015 
-0.011 

+0.013 
-0.018 

+0.010 
-0.017 
+0.019 
-0.012 

+0.005 
-0.016 
+0.041 
-0.034 

+0.009 
-0.013 
+0.006 
-0.003 
+0.082 
-0.012 
+0.004 
-0.005 
+0.003 
-0.003 



Table 38: Spin 
function of Q^. 



density matrix elements for 
Thenotag(|t| < 0.5 GeV^) 



the diffractive electroproduction of mesons, as a 
and tag < 3 GeV^) samples are combined. 
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= 3.3 GeV^ 



{W) (GeV) 


43 


59 


76 


92 


'00 





619 


± 





021 


+0 
-0 


012 
048 





600 


± 





021 


+0.012 
-0.020 





612 


± 





023 


+0.014 
-0.029 





580 


± 





028 


+0 
-0 


014 
027 


Re rlt 





580 


± 





028 


+0 
-0 


014 
027 





729 


± 





019 


+0.016 
-0.020 





755 


± 





019 


+0.020 
-0.024 





795 


± 





021 


+0 
-0 


023 
045 


».04 
' 1-1 





795 


± 





021 


+0 
-0 


023 
045 





731 


± 





021 


+0.021 
-0.021 





798 


± 





037 


+0.024 
-0.018 





878 


± 





036 


+0 
-0 


018 
027 


^1 

'oo 





878 


± 





036 


+0 
-0 


018 
027 





872 


± 





040 


+0.023 
-0.032 





856 


± 





045 


+0.041 
-0.059 





047 


± 





014 


+0 
-0 


012 
013 


r\. 





047 


± 





014 


+0 
-0 


012 
013 





056 


± 





015 


+0.014 
-0.010 





029 


± 





016 


+0.010 
-0.010 





028 


± 





019 


+0 
-0 


010 

oil 


Rerlo 





028 


± 





019 


+0 
-0 


010 

oil 





043 


± 





014 


+0.011 
-0.011 





028 


± 





014 


+0.012 
-0.012 


-0 


001 


± 





015 


+0 
-0 


012 
010 


ri-i 


-0 


001 


± 





015 


+0 
-0 


012 
010 





006 


± 





015 


+0.008 
-0.007 





047 


± 





027 


+0.009 
-0.010 





047 


± 





029 


+0 
-0 


009 
Oil 


^?o 





047 


± 





029 


+0 
-0 


009 
Oil 





033 


± 





028 


+0.008 
-0.015 





004 


± 





032 


+0.013 
-0.012 





006 


± 





016 


+0 
-0 


003 
009 


Im ri_i 





006 


± 





016 


+0 
-0 


003 
009 


-0 


Oil 


± 





017 


+0.005 
-0.002 





024 


± 





019 


+0.004 
-0.002 


-0 


012 


± 





022 


+0 
-0 


003 
006 


'oo 


-0 


012 


± 





022 


+0 
-0 


003 
006 


-0 


024 


± 





015 


+0.002 
-0.004 


-0 


005 


± 





015 


+0.004 
-0.003 


-0 


032 


± 





016 


+0 
-0 


004 
003 


'11 


-0 


032 


± 





016 


+0 
-0 


004 
003 


-0 


015 


± 





018 


+0.004 
-0.005 





021 


± 





027 


+0.020 
-0.017 


-0 


003 


± 





029 


+0 
-0 


002 
024 


Rerfo 


-0 


003 


± 





029 


+0 
-0 


002 
024 





029 


± 





030 


+0.006 
-0.020 





037 


± 





033 


+0.006 
-0.012 


-0 


015 


± 





100 


+0 
-0 


018 
064 


'1-1 


-0 


015 


± 





100 


+0 
-0 


018 
064 





134 


± 





147 


+0.022 
-0.030 





025 


± 





041 


+0.035 
-0.018 





018 


± 





040 


+0 
-0 


077 
012 


'^lO 





018 


± 





040 


+0 
-0 


077 

012 


-0 


082 


± 





118 


+0.019 
-0.011 


-0 


130 


± 





266 


+0.020 
-0.013 





163 


± 





143 


+0 
-0 


023 
027 


'1-1 





163 


± 





143 


+0 
-0 


023 
027 





154 


± 





263 


+0.032 
-0.012 


-0 


071 


± 





080 


+0.021 
-0.014 





120 


± 





291 


+0 
-0 


081 
036 



Table 39: Spin density matrix elements for the diffractive electroproduction of p mesons as a function of W, for 2.5 < < 5, 5 < 
15.5 and 15.5 < < 60 GeVl The notag {\t\ < 0.5 GeV^) and tag {\t\ < 3 GeV^) samples are combined. 



(g2)=7.5 GeV' 



{W) (GeV) 


58 


76 


93 


111 


'00 





729 


± 





019 


+0.016 
-0.020 





755 


± 





019 


+0.020 
-0.024 





795 


± 





021 


+0.023 
-0.045 





731 


± 





021 


+0.021 
-0.021 







731 


± 





021 


+0.021 
-0.021 





798 


± 





037 


+0.024 
-0.018 





878 


± 





036 


+0.018 
-0.027 





872 


± 





040 


+0.023 
-0.032 


^04 
'1-1 





872 


± 





040 


+0.023 
-0.032 





856 


± 





045 


+0.041 
-0.059 





047 


± 





014 


+0.012 
-0.013 





056 


± 





015 


+0.014 
-0.010 


'^OO 





056 


± 





015 


+0.014 
-0.010 





029 


± 





016 


+0.010 
-0.010 





028 


± 





019 


+0.010 
-0.011 





043 


± 





014 


+0.011 
-0.011 







043 


± 





014 


+0.011 
-0.011 





028 


± 





014 


+0.012 
-0.012 


-0 


001 


± 





015 


+0.012 
-0.010 





006 


± 





015 


+0.008 
-0.007 


Re rJo 





006 


± 





015 


+0.008 
-0.007 





047 


± 





027 


+0.009 
-0.010 





047 


± 





029 


+0.009 
-0.011 





033 


± 





028 


+0.008 
-0.015 


^Li 





033 


± 





028 


+0.008 
-0.015 





004 


± 





032 


+0.013 
-0.012 





006 


± 





016 


+0.003 
-0.009 


-0 


Oil 


± 





017 


+0.005 
-0.002 




-0 


Oil 


± 





017 


+0.005 
-0.002 





024 


± 





019 


+0.004 
-0.002 


-0 


012 


± 





022 


+0.003 
-0.006 


-0 


024 


± 





015 


+0.002 
-0.004 


Im r'l_^ 


-0 


024 


± 





015 


+0.002 
-0.004 


-0 


005 


± 





015 


+0.004 
-0.003 


-0 


032 


± 





016 


+0.004 
-0.003 


-0 


015 


± 





018 


+0.004 
-0.005 


'oo 


-0 


015 


± 





018 


+0.004 
-0.005 





021 


± 





027 


+0.020 
-0.017 


-0 


003 


± 





029 


+0.002 
-0.024 





029 


± 





030 


+0.006 
-0.020 







029 


± 





030 


+0.006 
-0.020 





037 


± 





033 


+0.006 
-0.012 


-0 


015 


± 





100 


+0.018 
-0.064 





134 


± 





147 


+0.022 
-0.030 


Re rfo 





134 


± 





147 


+0.022 
-0.030 





025 


± 





041 


+0.035 
-0.018 





018 


± 





040 


+0.077 
-0.012 


-0 


082 


± 





118 


+0.019 
-0.011 


rti 


-0 


082 


± 





118 


+0.019 
-0.011 


-0 


130 


± 





266 


+0.020 
-0.013 





163 


± 





143 


+0.023 
-0.027 





154 


± 





263 


+0.032 
-0.012 


' 10 





154 


± 





263 


+0.032 
-0.012 


-0 


071 


± 





080 


+0.021 
-0.014 





120 


± 





291 


+0.081 
-0.036 


-0 


032 


± 





033 


+0.028 
-0.026 


».6 

' 1-1 


-0 


032 


± 





033 


+0.028 
-0.026 


-0 


474 


± 





458 


+0.300 
-0.051 


-0 


055 


± 





101 


+0.031 
-0.009 


-0 


054 


± 





076 


+0.012 
-0.013 



Table 40: Spin density matrix elements for the diffractive electroproduction of p mesons as a function of W , for 2.5 < < h,h < 
15.5 and 15.5 < < go GeV^ continued from TableHH 



(g2) = 22.5 GeV^ 



{W) (GeV) 


72 


97 


122 


147 


'00 





798 


± 





037 


+0 
-0 


024 
018 





878 


± 





036 


+0.018 
-0.027 





872 


± 





040 


+0.023 
-0.032 





856 


± 





045 


+0.041 
-0.059 







856 


± 





045 


+0 
-0 


041 
059 





047 


± 





014 


+0.012 
-0.013 





056 


± 





015 


+0.014 
-0.010 





029 


± 





016 


+0.010 
-0.010 


^04 
'1-1 





029 


± 





016 


+0 
-0 


010 
010 





028 


± 





019 


+0.010 
-0.011 





043 


± 





014 


+0.011 
-0.011 





028 


± 





014 


+0.012 
-0.012 


'^OO 





028 


± 





014 


+0 
-0 


012 
012 


-0 


001 


± 





015 


+0.012 
-0.010 





006 


± 





015 


+0.008 
-0.007 





047 


± 





027 


+0.009 
-0.010 







047 


± 





027 


+0 
-0 


009 
010 





047 


± 





029 


+0.009 
-0.011 





033 


± 





028 


+0.008 
-0.015 





004 


± 





032 


+0.013 
-0.012 


Re rJo 





004 


± 





032 


+0 
-0 


013 
012 





006 


± 





016 


+0.003 
-0.009 


-0 


Oil 


± 





017 


+0.005 
-0.002 





024 


± 





019 


+0.004 
-0.002 







024 


± 





019 


+0 

-0 


004 
002 


-0 


012 


± 





022 


+0.003 
-0.006 


-0 


024 


± 





015 


+0.002 
-0.004 


-0 


005 


± 





015 


+0.004 
-0.003 




-0 


005 


± 





015 


+0 

-0 


004 
003 


-0 


032 


± 





016 


+0.004 
-0.003 


-0 


015 


± 





018 


+0.004 
-0.005 





021 


± 





027 


+0.020 
-0.017 


Im rj_]^ 





021 


± 





027 


+0 
-0 


020 
017 


-0 


003 


± 





029 


+0.002 
-0.024 





029 


± 





030 


+0.006 
-0.020 





037 


± 





033 


+0.006 
-0.012 


'oo 





037 


± 





033 


+0 
-0 


006 
012 


-0 


015 


± 





100 


+0.018 
-0.064 





134 


± 





147 


+0.022 
-0.030 





025 


± 





041 


+0.035 
-0.018 







025 


± 





041 


+0 
-0 


035 
018 





018 


± 





040 


+0.077 
-0.012 


-0 


082 


± 





118 


+0.019 
-0.011 


-0 


130 


± 





266 


+0.020 
-0.013 


Re rfo 


-0 


130 


± 





266 


+0 
-0 


020 
013 





163 


± 





143 


+0.023 
-0.027 





154 


± 





263 


+0.032 
-0.012 


-0 


071 


± 





080 


+0.021 
-0.014 


rti 


-0 


071 


± 





080 


+0 
-0 


021 
014 





120 


± 





291 


+0.081 
-0.036 


-0 


032 


± 





033 


+0.028 
-0.026 


-0 


474 


± 





458 


+0.300 
-0.051 


' 10 


-0 


474 


± 





458 


+0 
-0 


300 
051 


-0 


055 


± 





101 


+0.031 
-0.009 


-0 


054 


± 





076 


+0.012 
-0.013 





001 


± 





028 


+0.014 
-0.021 


».6 

' 1-1 





001 


± 





028 


+0 

-0 


014 
021 


-0 


045 


± 





055 


+0.012 
-0.035 





009 


± 





074 


+0.011 
-0.012 


-0 


047 


± 





206 


+0.010 
-0.014 



Table 41: Spin density matrix elements for the diffractive electroproduction of p mesons as a function of W , for 2.5 < < h,h < 
15.5 and 15.5 < < go GeV^ continued from TableHH 



\t\) (GeV^) 



0.08 \ 0^ 

(g2) = 3.3 GeV^ 



1.05 



r04 
' 00 





602 


± 





014 


+0.011 
-0.031 





593 


± 





020 


+0.019 
-0.029 





660 


± 





040 


+0.033 
-0.045 


Re r% 





030 


± 





009 


+0.011 
-0.012 





026 


± 





014 


+0.013 
-0.009 





150 


± 





032 


+0.017 
-0.028 


' 1-1 


-0 


004 


± 





oil 


+0.002 
-0.004 





037 


± 





017 


+0.006 
-0.002 


-0 


080 


± 





036 


+0.010 
-0.015 


'"oo 





070 


± 





077 


+0.020 
-0.014 


-0 


198 


± 





277 


+0.259 
-0.025 





163 


± 





119 


+0.067 
-0.076 


ril 


-0 


028 


± 





040 


+0.006 
-0.011 





029 


± 





171 


+0.025 
-0.132 


-0 


144 


± 





106 


+0.015 
-0.015 


Rerlo 


-0 


010 


± 





014 


+0.008 
-0.006 


-0 


086 


± 





020 


+0.015 
-0.017 


-0 


126 


± 





045 


+0.028 
-0.016 


^1-1 





143 


± 





016 


+0.006 
-0.009 





149 


± 





024 


+0.011 
-0.013 





075 


± 





048 


+0.010 
-0.022 


' 10 





018 


± 





013 


+0.012 
-0.012 





072 


± 





021 


+0.016 
-0.016 





104 


± 





046 


+0.017 
-0.025 


Im r^_i 


-0 


192 


± 





015 


+0.007 
-0.006 


-0 


108 


± 





023 


+0.016 
-0.008 


-0 


207 


± 





056 


+0.018 
-0.022 


'00 





125 


± 





020 


+0.032 
-0.060 





199 


± 





030 


+0.038 
-0.047 





197 


± 





060 


+0.051 

-0.033 


' 11 


-0 


014 


± 





014 


+0.024 
-0.007 





004 


± 





022 


+0.025 
-0.028 





100 


± 





043 


+0.017 
-0.022 


Rer^o 





160 


± 





006 


+0.003 
-0.004 





154 


± 





010 


+0.008 
-0.009 





138 


± 





024 


+0.015 
-0.207 







005 


± 





008 


+0.001 
-0.003 





024 


± 





013 


+0.003 
-0.005 


-0 


065 


± 





030 


+0.014 
-0.006 


Mo 


-0 


162 


± 





006 


+0.005 
-0.004 


-0 


130 


± 





009 


+0.008 
-0.009 


-0 


131 


± 





020 


+0.012 
-0.010 


' 1-1 





003 


± 





008 


+0.002 
-0.002 


-0 


012 


± 





Oil 


+0.004 
-0.004 


-0 


006 


± 





025 


+0.011 
-0.004 














iQ') = ^ 


^.6 GeV^ 


















^04 
'00 





734 


± 





Oil 


+0.014 
-0.022 





817 


± 





014 


+0.031 
-0.040 





841 


± 





016 


+0.069 
-0.067 


Rer?^ 





030 


d= 





008 


+0.009 
-0.008 





004 


=fc 





Oil 


+0.016 
-0.014 


-0 


049 


± 





019 


+0.014 
-0.013 


^04 
' 1-1 


-0 


006 


± 





008 


+0.001 
-0.001 


-0 


008 


± 





012 


+0.004 
-0.006 


-0 


061 


± 





019 


+0.015 
-0.016 


^00 


-0 


009 


± 





032 


+0.008 
-0.007 


-0 


052 


± 





067 


+0.049 
-0.023 


-0 


046 


± 





063 


+0.070 
-0.042 




-0 


013 


± 





030 


+0.004 
-0.004 





012 


± 





039 


+0.026 
-0.026 





008 


± 





038 


+0.074 
-0.073 


Re tIq 


-0 


028 


± 





Oil 


+0.009 
-0.007 





045 


± 





016 


+0.014 
-0.020 


-0 


094 


± 





027 


+0.024 
-0.020 


^i-i 





133 


± 





012 


+0.004 
-0.005 





102 


d= 





016 


+0.010 
-0.010 


-0 


007 


± 





027 


+0.018 
-0.019 


'^fo 


-0 


036 


± 





Oil 


+0.007 
-0.008 


-0 


Oil 


± 





015 


+0.013 
-0.013 


-0 


044 


± 





025 


+0.019 
-0.018 


Im r^_i 


-0 


081 


± 





Oil 


+0.002 
-0.003 


-0 


077 


± 





016 


+0.006 
-0.022 


-0 


129 


± 





026 


+0.022 
-0.017 


'00 





071 


± 





017 


+0.029 
-0.028 





169 


± 





024 


+0.051 
-0.055 





115 


± 





042 


+0.062 
-0.046 


^11 





023 


± 





012 


+0.005 
-0.006 


-0 


036 


± 





017 


+0.020 
-0.021 


-0 


010 


± 





029 


+0.042 
-0.052 


Re r^o 





146 


± 





005 


+0.003 
-0.003 





137 


± 





008 


+0.008 
-0.008 





152 


± 





013 


+0.023 
-0.021 




-0 


013 


± 





006 


+0.002 
-0.003 


-0 


007 


± 





009 


+0.004 
-0.004 





060 


± 





016 


+0.029 
-0.053 


^6 
'10 


-0 


145 


± 





005 


+0.003 
-0.003 


-0 


113 


± 





007 


+0.006 
-0.007 


-0 


103 


± 





Oil 


+0.014 
-0.015 


^6 
' 1-1 


-0 


020 


± 





006 


+0.001 
-0.002 





017 


± 





007 


+0.003 
-0.003 


-0 


029 


± 





012 


+0.002 
-0.005 



Table 42: Spin density matrix elements for the diffractive electroproduction of p mesons as a 
function of \t\, for 2.5 < < 5 and 5 < < gO GeVl The notag (\t\ < 0.5 GeV^) and tag 
(\t\ < 3 GeV^) samples are combined. 



103 



( 


t ) (GeV^) 


0.08 


0.34 


1.05 







667 


± 





019 


+0.012 
-0.011 





641 


± 





027 


+0.021 
-0.021 





830 


± 





039 


+0.087 
-0.088 


Re r?o^ 





002 


± 





014 


+0.010 
-0.012 





008 


± 





020 


+0.020 
-0.021 


-0 


010 


± 





041 


+0.033 
-0.026 


'1-1 


-0 


014 


± 





016 


+0.004 
-0.002 


-0 


003 


± 





023 


+0.010 
-0.006 


-0 


036 


± 





044 


+0.012 
-0.024 


'^00 


-0 


058 


± 





108 


+0.005 
-0.007 


-0 


069 


± 





166 


+0.037 
-0.020 





007 


± 





032 


+0.140 
-0.127 


ril 


-0 


012 


± 





079 


+0.005 
-0.002 


-0 


044 


± 





138 


+0.020 
-0.015 


-0 


029 


± 





040 


+0.068 
-0.065 


Rerlo 


-0 


028 


± 





019 


+0.016 
-0.009 


-0 


071 


± 





027 


+0.016 
-0.020 


-0 


139 


± 





059 


+0.026 
-0.025 


rl-i 





215 


± 





021 


+0.009 
-0.007 





064 


± 





032 


+0.010 
-0.013 





114 


± 





067 


+0.060 
-0.031 


Mo 





014 


± 





019 


+0.015 
-O.Ol,'? 


-0 


018 


± 





027 


+0.017 
-0.019 





094 


± 





053 


+0.037 
-0.021 


Im r^_j^ 


-0 


132 


± 





022 


+0.007 
-0.008 


-0 


129 


± 





033 


+0.008 
-0.013 


-0 


204 


± 





060 


+0.0.36 
-0.045 


'oo 





037 


± 





029 


+0.025 
-0.026 





099 


± 





042 


+0.046 
-0.048 





138 


± 





097 


+0.071 
-0.081 




-0 


002 


± 





020 


+0.009 
-0.010 





049 


± 





030 


+0.015 
-0.018 





006 


± 





066 


+0.044 
-0.053 


Rerfo 





176 


± 





009 


+0.00,3 
-0.004 





155 


± 





013 


+0.009 
-0.008 





094 


± 





032 


+0.021 
-0.030 




-0 


015 


± 





Oil 


+0.00,3 
-0.003 





043 


± 





018 


+0.006 
-0.009 


-0 


028 


± 





038 


+0.030 
-0.029 


^6 

Mo 


-0 


163 


± 





009 


+0.003 
-0.004 


-0 


132 


± 





012 


+0.007 
-0.009 


-0 


128 


± 





023 


+0.026 
-0.023 


' 1-1 


-0 


014 


± 





Oil 


+0.004 
-0.003 





005 


± 





015 


+0.005 
-0.003 


-0 


037 


± 





027 


+0.027 
-0.012 



Table 43: Spin density matrix elements for the diffractive electroproduction of mesons as a 
function of The notag {\t\ < 0.5 GeV^) and tag < 3 GeV^) samples are combined. 



104 



= 3.3 GeV^ 



(m^+^-) (GeV) 



0.65 



0.75 



0.85 



0.95 



1.05 



T.04 
'00 


0.677 ± 0.021 


Rer?^ 


0.067 ± 0.016 


' 1-1 


-0.002 ± 0.019 


'^00 


-0.033 ± 0.091 


ril 


-0.029 ± 0.079 


Rerlo 


0.020 ± 0.023 


ri-i 


0.093 ± 0.028 


r?o 


0.025 ± 0.022 


Im ri_i 


-0.122 ± 0.027 


'go 


0.139 ± 0.036 


' 11 


0.005 ± 0.025 


Re rfo 


0.147 ± 0.011 




0.018 ± 0.016 


'^lO 


-0.129 ±0.010 


' 1-1 


-0.005 ± 0.013 



+0.029 
-0.037 
+0.013 
-0.012 

+0.009 
-0.003 
+0.026 
-0.033 
+0.027 
-0.027 
+0.009 
-0.013 

+0.011 
-0.010 

+0.009 
-0.010 
+0.021 
-0.010 
+0.046 
-0.059 

+0.025 
-0.021 

+0.009 
-0.011 
+0.004 
-0.009 
+0.009 
-0.012 
+0.002 
-0.005 



0.635 ± 
0.028 ± 
0.002 ± 
0.068 ± 

-0.033 ± 

-0.069 ± 
0.135 ± 
0.034 ± 

-0.156 ± 
0.142 ± 
0.010 ± 
0.151 ± 
0.000 ± 

-0.149 ± 
0.006 ± 



0.017 
0.012 
0.013 
0.070 
0.035 
0.017 
0.019 
0.016 
0.017 
0.024 
0.017 
0.007 
0.009 
0.007 
0.009 



+0.012 
-0.021 
+0.011 
-0.013 
+0.002 
-0.002 
+0.084 
-0.017 
+0.011 
-0.041 

+0.011 
-0.011 

+0.005 
-0.007 

+0.013 
-0.010 
+0.009 
-0.006 
+0.049 
-0.034 

+0.006 
-0.013 

+0.003 
-0.005 
+0.002 
-0.002 

+0.005 
-0.004 
+0.002 
-0.004 



0.559 ± 
0.033 ± 

-0.016 ± 
0.084 ± 

-0.058 ± 

-0.043 ± 
0.188 ± 
0.058 ± 

-0.171 ± 
0.302 ± 

-0.084 ± 
0.181 ± 
0.003 ± 

-0.155 ± 

-0.009 ± 



0.023 
0.016 
0.019 
0.070 
0.042 
0.022 
0.026 
0.023 
0.027 
0.036 
0.025 
0.011 
0.014 
0.010 
0.013 



+0.013 
-0.046 
+0.011 
-0.012 
+0.00.5 
-0.007 
+0.045 
-0.054 
+0.025 
-0.019 
+0.009 
-0.009 

+0.011 
-0.008 

+0.013 
-0.020 
+0.011 
-0.007 
+0.032 
-0.284 

+0.145 
-0.028 

+0.006 
-0.004 

+0.003 
-0.008 
+0.006 
-0.007 
+0.002 
-0.004 



0.456 ± 
0.046 ± 
0.038 ± 
-0.089 ± 
-0.062 ± 
-0.035 ± 
0.175 ± 
0.022 ± 
-0.286 ± 
0.202 ± 
-0.005 ± 
0.167 ± 
0.013 ± 
-0.149 ± 
0.005 ± 



0.037 
0.023 
0.032 
0.223 
0.187 
0.036 
0.045 
0.036 
0.042 
0.049 
0.037 
0.017 
0.025 
0.016 
0.020 



+0.024 
-0.012 
+0.013 
-0.010 
+0.008 
-0.010 
+0.061 
-0.029 
+0.015 
-0.034 
+0.014 
-0.014 

+0.033 
-0.021 

+0.015 
-0.023 
+0.017 
-0.026 
+0.044 
-0.032 

+0.028 
-0.034 

+0.008 
-0.009 
+0.011 
-0.007 
+0.013 
-0.018 
+0.003 
-0.005 



0.458 ± 
0.027 ± 
0.003 ± 
0.212 ± 
-0.135 ± 
-0.003 ± 
0.113 ± 
-0.041 ± 
-0.297 ± 
0.263 ± 
-0.027 ± 
0.084 ± 
0.005 ± 
-0.163 ± 
0.002 ± 



0.055 
0.034 
0.042 
0.186 
0.101 
0.051 
0.058 
0.047 
0.060 
0.071 
0.051 
0.025 
0.032 
0.021 
0.031 



+0.027 
-0.051 
+0.014 
-0.016 
+0.014 
-0.014 
+0.060 
-0.077 
+0.030 
-0.037 
+0.038 
-0.018 

+0.080 
-0.057 
+0.012 
-0.016 
+0.056 
-0.037 
+0.079 
-0.130 
+0.041 
-0.031 
+0.022 
-0.034 
+0.010 
-0.010 
+0.020 
-0.015 
+0.031 
-0.007 



Table 44: Spin density matrix elements for the diffractive electroproduction of p mesons as a function of mTrjr, for 2.5 < < b and 
5 < < 60 GeVl The notag {\t\ < 0.5 GeV^) and tag < 3 GeV^) samples are combined. 



(g2) = 8.6 GeV^ 



(m^+n-) (GeV) 


0.65 


0.75 









.8 


5 




0.95 


1.05 


'00 





857 


± 





012 


+0.048 
-0.049 





793 


± 





012 


+0.019 
-0.030 





747 


± 





016 


+0.032 
-0.032 





586 


± 





027 


+0.031 
-0.028 





473 


± 





035 


+0.021 
-0.026 


Rer% 


-0 


002 


± 





013 


+0.015 
-0.014 





016 


± 





009 


+0.011 
-0.010 





028 


± 





012 


+0.009 
-0.012 


-0 


039 


± 





018 


+0.012 
-0.009 





125 


± 





024 


+0.021 
-0.015 


' 1-1 


-0 


051 


± 





014 


+0.007 
-0.008 





005 


± 





009 


+0.002 
-0.001 





033 


± 





013 


+0.002 
-0.004 


-0 


063 


± 





022 


+0.012 
-0.010 


-0 


116 


± 





031 


+0.024 
-0.025 


^00 


-0 


053 


± 





109 


+0.022 
-0.036 


-0 


008 


± 





025 


+0.026 
-0.010 


-0 


020 


± 





037 


+0.031 
-0.022 





047 


± 





137 


+0.105 
-0.033 





156 


± 





263 


+0.074 
-0.052 




-0 


019 


± 





084 


+0.042 
-0.032 


-0 


009 


± 





022 


+0.016 
-0.018 


-0 


004 


± 





028 


+0.024 
-0.026 


-0 


141 


± 





180 


+0.025 
-0.054 





013 


± 





184 


+0.041 
-0.040 


RerJo 





038 


± 





020 


+0.017 
-0.014 


-0 


028 


± 





012 


+0.012 
-0.010 


-0 


027 


± 





016 


+0.012 
-0.010 


-0 


061 


± 





025 


+0.012 
-0.014 





066 


± 





035 


+0.042 
-0.009 







117 


± 





020 


+0.017 
-0.016 





087 


± 





013 


+0.004 
-0.008 





121 


± 





018 


+0.009 
-0.008 





251 


± 





032 


+0.128 
-0.039 


-0 


024 


± 





039 


+0.040 
-0.016 




-0 


044 


± 





017 


+0.014 
-0.010 


-0 


046 


± 





012 


+0.012 
-0.010 





004 


± 





017 


+0.009 
-0.011 





025 


± 





025 


+0.004 
-0.010 


-0 


034 


± 





032 


+0.021 
-0.009 


Im r'l_^ 


-0 


046 


± 





018 


+0.013 
-0.006 


-0 


084 


± 





012 


+0.004 
-0.005 


-0 


087 


± 





018 


+0.010 
-0.008 


-0 


104 


± 





030 


+0.012 
-0.020 


-0 


146 


± 





039 


+0.039 
-0.019 


'oo 





106 


± 





029 


+0.042 
-0.044 





098 


± 





019 


+0.038 
-0.041 





155 


± 





026 


+0.037 
-0.038 





075 


± 





039 


+0.031 
-0.038 





035 


± 





047 


+0.034 
-0.039 




-0 


002 


± 





020 


+0.025 
-0.028 





010 


± 





013 


+0.011 
-0.012 


-0 


016 


± 





018 


+0.020 
-0.019 





019 


± 





029 


+0.029 
-0.029 





060 


± 





035 


+0.020 
-0.022 


Re rfo 





173 


± 





010 


+0.017 
-0.016 





128 


± 





006 


+0.004 
-0.005 





142 


± 





008 


+0.008 
-0.008 





113 


± 





014 


+0.014 
-0.011 





161 


± 





015 


+0.018 
-0.015 


rti 





005 


± 





Oil 


+0.005 
-0.005 


-0 


017 


± 





007 


+0.003 
-0.003 


-0 


001 


± 





009 


+0.003 
-0.003 





007 


± 





018 


+0.010 
-0.007 


-0 


068 


± 





023 


+0.013 
-0.015 


».6 
' 10 


-0 


115 


± 





008 


+0.009 
-0.010 


-0 


133 


± 





006 


+0.004 
-0.004 


-0 


147 


± 





008 


+0.008 
-0.007 


-0 


124 


± 





Oil 


+0.010 
-0.012 


-0 


107 


± 





014 


+0.009 
-0.010 


»,6 
' 1-1 





013 


± 





009 


+0.003 
-0.003 


-0 


016 


± 





006 


+0.002 
-0.001 


-0 


023 


± 





008 


+0.002 
-0.004 





014 


± 





013 


+0.005 
-0.010 





012 


± 





020 


+0.019 
-0.009 



Table 45: Spin density matrix elements for the diffractive electroproduction of p mesons as a function of m^7r, for 2.5 < < 5 and 
5 < < 60 GeV^ continued from Table 113 



(g2) (GeV^) 


+ 2rf^ 


r^o + 2rl, 


2.9 


0.149 


± 0.018 


+0.062 
-0.060 


-0.041 


± 0.034 


+0.029 
-0.034 


4.1 


0.144 


± 0.019 


+0.046 
-0.055 


-0.036 


± 0.041 


+0.023 
-0.025 


6.6 


0.122 


± 0.016 


+0.049 
-0.049 


-0.078 


±0.032 


+0.023 
-0.014 


11.9 


0.088 


± 0.024 


+0.049 
-0.047 


-0.093 


± 0.045 


+0.024 
-0.018 


18.5 


0.178 


± 0.031 


+0.037 
-0.049 


-0.010 


± 0.059 


+0.034 
-0.034 


31.3 


0.149 


± 0.028 


+0.041 
-0.037 


0.123 


± 0.059 


+0.030 
-0.032 


{\t\) (GeV^) 




r^o + 2rl, 


0.03 


0.074 


± 0.022 


+0.032 
-0.036 


0.003 


± 0.045 


+0.004 
-0.013 


0.13 


0.109 


± 0.022 


+0.062 
-0.052 


-0.024 


± 0.044 


+0.017 
-0.024 


0.31 


0.200 


± 0.028 


+0.064 
-0.086 


-0.131 


± 0.054 


+0.041 
-0.028 


0.70 


0.249 


± 0.051 


+0.063 
-0.038 


-0.015 


± 0.099 


+0.040 
-0.067 


1.44 


0.308 


± 0.077 


+0.086 
-0.096 


-0.162 


± 0.129 


+0.086 
-0.080 



Table 46: and \t\ dependences of the matrix element combinations Tqq ± 2r^]^ and Tqq ± 2rl^, 
obtained from fits of Eq. (|4T]) to the (p distribution, for p meson electroproduction. The notag 
{\t\ < 0.5 GeV^) and tag < 3 GeV^) samples are combined. 



(g2) (GeV^) 


-Pnpe.t 


{\t\){GeV') 


-Pnpe.t 


p production 


3.1 
4.8 
7.5 
12.2 
23.1 


1.13 ±0.07 t'om 
1.24 ± 0.09 1°:°^ 
1.03 ±0.10 to^ol 
0.95 ±0.11 1°:°^ 
1.07 ±0.16 to^ot 


0.03 
0.13 
0.31 
0.70 
1.44 


1.10 ±0.06 tof2 

1.02 ± 0.12 tofs 

1.03 ±0.12 t'ofs 
0.96 ± 0.39 toll 
1.08 ±0.12 tofi 


(p production 


3.3 
6.6 
15.8 


0.98 ± 0.31 to^ot 
1.13 ±0.11 t'om 
0.99 ± 0.10 torn 


0.08 
0.34 
1.05 


0.90 ±0.17 tof2 
1.31 ±0.29 torn 
0.93 ±0.18 toll 



Table 47: Asymmetry Pnpe.t between natural and unnatural parity exchange for transverse 
photons, as a function of and \t\, for p and meson production. 
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cos 6 


p production 


3.3 
6.6 
15.8 


0.914 ± 0.014 to'^oll 
0.915 ± 0.026 to^oll 
0.978 ± 0.030 to^oll 


(p production 


5.3 


0.966 ± 0.027 tofil 



Table 48: Cosine of the phase 6 between the Too and Tn helicity conserving amplitudes for p 
and (p meson production, measured as a function of from two-dimensional fits to Eq. (1471) . 
in the SCHC approximation < 0.5 GeV^) 



(g2) (GeV^) 


R = ctl/o-t 


p production 


3.1 
4.8 
7.5 
12.2 
23.1 


1 qc +0.08 +0.11 
-0.08 -0.13 

1 Q9 +0.16 +0.15 
^■•^■^ -0.15 -0.12 
q ecr +0.31 +0.41 
O.UO _o.27 -0.38 

q an +0.38 +0.39 
'^■^^ -0.34 -0.37 

c-9 +1.27 +1.23 
"•'-'^ -0.98 -0.96 


production 


3.3 
6.6 
15.8 


1 q7 +0.14 +0.09 
-0.13 -0.09 

9 87 +0A0 +0.28 
^•O' -0.34 -0.26 

c m +2.00 +1.00 

-1.33 -0.78 



Table 49: dependence of the ratio R = ai/ax of the longitudinal to transverse cross 
sections, for p and (p meson production. 
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{W) (GeV) 


R = ul/ot 


(g2) = 3.3 GeV^ 


43 
59 
76 
92 


-"-■"J" -0.14 -0.08 

1 AR +0-14 +0.09 
^■^^ -0.13 -0.08 
1 +0.16 +0.10 
^■^'^ -0.15 -0.08 
-1 97 +0.17 +0.10 
-0.15 -0.11 


= 7.5 GeV' 


58 
76 
93 
111 


9 f^7 +0.28 +0.22 
_Q24 _o.l8 

q nn +0-34 +0.29 
O.UU _Q29 _o.30 

q qi +0.47 +0.46 
-0.41 -0.45 
9 79 +0.32 +0.26 
^•'^ -0.28 -0.22 


(g2) = 22.5 GeV' 


72 
97 
122 
147 


q CTQ +0.98 +0.59 

O.OJ _Q74 _Q3y 

r ofi +2.48 +1.19 
-1.60 -1.05 
97 +2.98 +1.41 
-1.76 -1.02 

Ai +3.13 +1.08 
U.Ul _2.i4 



Table 50: W dependence of the ratio i? = (Ti/cTT of the longitudinal to transverse cross sections 
for p meson production, for 2.5 < Q' < 5 GeV', 5 < Q' < 15.5 GeV' and 15.5 < Q' < 

60 GeVl 



(10 (GeV) 


R = gl/ot 


= 3.3 GeV' 


0.08 
0.34 

1.05 


-1 AO +0.09 +0.08 
-'-■^'J -0.09 -0.08 

-1 9R +0.13 +0.16 
^■^^ -0.12 -0.16 

-0.32 -0.23 


(g2) = 8.6 GeV' 


0.08 
0.34 
1.05 


9 7q +0.16 +0.15 
^■'"J -0.15 -0.15 
q AO +0.40 +0.92 
-0.35 -0.79 

A QQ +0-69 +3.63 
^•^^ -0.59 -1.98 



Table 51: |t | dependence of the ratio R = ai/cT of the longitudinal to transverse cross sections 
for p meson production, for 2.5 < Q' < 5 GeV' and 5 < Q' < 60 GeVl 
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(^TT+TT-) (GeV) 


R = (7l/ (7t 


(g2) = 3.3 GeV^ 


0.65 
0.75 
0.85 
0.95 
1.05 


1 Q7 +0.21 +0.32 
-0.19 -0.30 

1 +0-13 +0.11 
-■-■"^ -0.12 -0.13 

n QC +0.14 +0.12 
-0.15 -0.09 

n RQ +0-14 +0.05 
'-'■"^ -0.15 -0.10 
A rn +0.24 +0.14 
-0.34 -0.19 


(Q2) = 8.6 GeV^ 


0.65 
0.75 
0.85 
0.95 
1.05 


f- 7n +0.59 +2.85 
-0.52 -1.75 

q 7n +0.28 +0.46 
' -0.26 -0.43 
9 79 +0.25 +0.40 
^•'^ -0.23 -0.36 
1 qo +0.17 +0.20 
-'-•'-'^ -0.15 -0.17 

u.yu _Qj^2 -0.06 



Table 52: TOtttt dependence of the ratio R = ai/cTT of the longitudinal to transverse cross 
sections for p meson production, for 2.5 < < 5 GeV^ and 5 < < go GeVl 



(m^+n-) (GeV) 






iQ')- 


= 3.3 GeV^ 




0.67 




7.33 ± 0.51 


+0.38 
-0.36 


0.87 




8.25 ± 0.40 


+0.26 
-0.29 


1.07 




7.73 ± 0.51 


+0.46 
-0.42 




iQ')- 


= 9.0 GeV^ 




0.67 




7.15 ± 0.48 


+0.27 
-0.30 


0.87 




7.65 ± 0.43 


+0.30 
-0.32 


1.07 




6.83 ± 0.45 


+0.27 
-0.35 



Table 53: Dependence of the exponential t slope for p elastic production, as a function of the 
mass m^^, for two domains in Q"^: 2.5 < < 5 QeV^ and 5 < < 60 GeV^ 
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{Q') (GeV^) 


Til /Too 


Toi/Too 


Tio/Too 


T-u/Too 


cos(0ii - 0oo) 


p production 


3.1 

4.8 
7.5 
12.2 
23.1 


0.786 ± 0.019 t'^omo 
0.619 ± 0.021 toml 
0.511 ± 0.015 t^omi 
0.444 ± 0.018 1^08 
0.287 ± 0.020 tol^oot 


0.177 ± 0.019 

0.169 ± 0.022 toml 
0.056 ± 0.018 toml 
0.106 ± 0.022 toml 
0.090 ± 0.025 


-0.002 ± 0.011 

-0.025 ± 0.014 tofol 
0.018 ± 0.010 tomt 

-0.049 ± 0.013 tomr 
0.047 ± 0.016 


-0.019 ± 0.014 toml 

0.001 ±0.017 1[5:E5[J5 

-0.034 ± 0.014 tomt 
-0.038 ± 0.017 tomo 
0.029 ± 0.020 tom5 


0.949 ± 0.024 toml 
0.882 ± 0.035 f^oml 
0.997 ± 0.038 toml 
0.900 ± 0.049 toml 
0.879 ± 0.094 toml 


4> production 


3.4 
6.6 
15.7 


0.823 ± 0.033 tofit 
0.580 ± 0.029 torn! 
0.375 ± 0.035 to^ou 


0.049 ± 0.031 to!^o^ 

0.115 ± 0.035 

0.097 ± 0.041 tofrl 


-0.027 ± 0.020 toml 
-0.001 ± 0.020 1°:°09 
-0.028 ± 0.027 12:02! 


-0.151 ± 0.027 t^oml 
0.076 ± 0.026 t^oxm 
0.008 ± 0.027 tol^oll 


0.982 ± 0.048 tomo 
0.993 ± 0.057 toml 
0.867 ±0.115 1^:25? 



Table 54: Ratios of the helicity amplitudes (taken to be purely imaginary) and phase difference between the Tn and Too amplitudes (the 
amplitude ratios Tio/Too and T_ii/Too and the phase difference 0oi — 0oo are taken to be 0) for p and (j) meson production, computed from 
global fits to the measurements of the 15 spin density matrix elements, as a function of Q'^ (NPE is assumed). 



( t ) (GeV^) 








T-u/Tqq 


cos(0ii - 4>oo) 


p production 


(Q2) = 3.3 GcV^ 


0.08 
0.34 
1.05 


0.773 ± 0.020 to^of^ 
0.769 ± 0.031 to^oll 
0.673 ± 0.067 


0.128 ± 0.020 to^ofr 
0.281 ± 0.034 1°:°^° 
0.248 ± 0.070 toml 


-0.007 ± 0.012 1°:°^^ 
-0.062 ± 0.019 tal]^ 
0.138 ± 0.046 to^^ 


-0.004 ± 0.015 to^Q^ol 
-0.015 ± 0.026 IHra 
-0.095 ± 0.051 tomt 


0.939 ± 0.026 toml 
0.851 ± 0.040 to^oll 
0.929 ± 0.094 toill 


(g2) = 8.6 GeV^ 


0.08 
0.34 
1.05 


0.559 ± 0.013 tofil 
0.445 ± 0.016 
0.422 ± 0.021 to^r 


0.059 ± 0.014 to^ll 
0.104 ± 0.020 toM6 
0.113 ± 0.033 toml 


0.012 ± 0.009 t^^^^ 
-0.007 ± 0.012 
-0.103 ± 0.019 t^l^ 


-0.009 ± 0.011 
-0.037 ± 0.016 tomt 
-0.090 ± 0.023 tomb 


0.947 ± 0.026 
0.894 ± 0.044 tomt 
0.955 ± 0.076 toill 


4> production 


0.08 
0.34 
1.05 


0.713 ± 0.024 tofil 
0.650 ± 0.036 tofil 
0.478 ± 0.053 toilr 


0.046 ± 0.025 to^il 
0.117 ± 0.039 toml 
0.164 ± 0.065 toill 


-0.007 ± 0.016 
-0.033 ± 0.025 
-0.056 ± 0.044 


-0.029 ± 0.021 tofol 
-0.050 ± 0.031 tofog 
-0.015 ± 0.049 to^oll 


0.969 ± 0.042 
0.869 ± 0.057 tofi 
0.851 ± 0.146 



Table 55: Ratios of the helicity amplitudes (taken to be purely imaginary) and phase difference between the Tn and Too amplitudes (the 
amplitude ratios Tio/Too and T-u/Tqo and the phase difference 0oi — 0oo are taken to be 0), computed from global fits to the measurements 
of the 15 spin density matrix elements, as a function of separately for 2.5 < < 5 GeV^ and 5 < < 60 GeV^ for p meson 
production and for 2.5 < Q"^ < 60 GeV^ for production (NPE is assumed). 



(^TT+TT-) (GeV) 


Til/ Too 








COs(0ii - 0oo) 


= 3.3 GeV^ 


0.65 
0.75 
0.85 
0.95 
1.05 


0.618 ± 0.029 toml 
0.711 ± 0.023 to^ll 
0.861 ± 0.036 iHio 
1.066 ± 0.073 1^38 
0.981 ± 0.106 1^66 


0.138 ± 0.034 tolfb 
0.166 ± 0.024 
0.292 ± 0.035 tof^ 
0.250 ± 0.060 1^49 
0.275 ± 0.090 


0.020 ± 0.021 

-0.006 ± 0.014 1°:°;^ 

-0.058 ± 0.020 in^^ 
-0.018 ± 0.035 l^ig 
-0.077 ± 0.057 t^^i 


-0.031 ± 0.025 t^^fi 

-0.007 ± 0.018 

-0.062 ± 0.026 1^06 
0.009 ± 0.048 1^27 
0.050 ± 0.063 t^^i 


0.855 ± 0.047 t^^l 
0.902 ± 0.033 
0.971 ± 0.042 
0.898 ± 0.065 1^56 
0.722 ± 0.091 to^Q 


(Q2) = 8.6 GeV^ 


0.65 
0.75 
0.85 
0.95 
1.05 


0.429 ± 0.015 to'^ofg 
0.472 ± 0.013 t'^oml 
0.558 ± 0.019 f^omi 
0.746 ±0.041 t^oll 
0.717 ± 0.052 toMl 


0.050 ± 0.022 
0.069 ± 0.015 t^^^l 
0.145 ± 0.023 

0.110 ± 0.041 

-0.062 ± 0.045 t^MO 


-0.006 ± 0.013 t'^^m^ 

0.001 ± 0.009 
-0.017 ± 0.013 to^^ 
-0.066 ± 0.027 toml 

0.261 ± 0.038 t^^^^ 


-0.095 ± 0.017 t^Q^^^i 
0.006 ± 0.012 1°:°°^ 
0.028 ± 0.017 

-0.101 ± 0.033 

-0.183 ±0.0411°:°^^ 


0.997 ± 0.052 1°:°°^ 
0.923 ± 0.034 t^lf^ 
0.959 ± 0.041 
0.690 ± 0.049 
0.761 ± 0.059 to^f^^ 



Table 56: Ratios of the helicity amplitudes (taken to be purely imaginary) and phase difference between the Tn and Tqq amplitudes (the 
amplitude ratios Tio/Too and T_ii/Too and the phase difference 0oi — 0oo are taken to be 0), computed from global fits to the measurements 
of the 15 spin density matrix elements, as a function of m^Tr separately for 2.5 < < 5 GeV^ and 5 < < 60 GeV^ for p meson 
production and for 2.5 < < 60 GeV^ for production (NPE is assumed). 



